ANNALEN DER PHYSIK 


- HEFT 3+- JUNI 1939 


Uber die Raumverteilung der Trägerwolken 2 


A: aus Photoabsorption und Streuabsorption 
bei der Luftionisation durch Réntgenstrahlen -¥ 
(Aus dem Institut für Medizinische Physik der Universität Göttingen) A A 

4 

(Mit 13 Abbildungen) 

In einer Untersuchung über diese Fragestellung. die kürzlich !) d 

in dieser Zeitschrift erschien und im folgenden kurz mit (I) be- af 
zeichnet werden soll, konnte der Verf. sowohl über den Anteil beider x 
Arten von Trägerbildung an der Gesamtionisation als auch über ” 


deren räumliche Verteilung quantitative Angaben machen. Als 
experimentelle Unterlage dienten Untersuchungen des Verf.?2) über 
die Luftionisation in zylindrischen Faßkammern von gleicher wirk- 
samer Länge, aber von verschiedenem Radius, die mit extrem stark 
gefilterten und daher weitgehend homogenisierten Röntgenstrahlen 
durchgeführt worden waren. Drückt man die bei einer bestimmten 
Strahlung in Kammern von verschiedenem Radius o beobachteten 


b 


Ionisationen als Bruchteil derjenigen Ionisation aus, die eine Kammer +4 
von praktisch unendlich großem Radius liefert*), und trägt man diese . Br 
- 


relative Ionisation in Abhängigkeit vom Radius o auf, so erhält man 
eine „Integralkurve“, die angibt, welcher Bruchteil aller überhaupt 


durch die gegebene Strahlung erzeugbaren Träger innerhalb eines _ 
mit dem Strahlenbündel koaxialen Zylinders o entsteht. Führt . 
man dasselbe auch für andere Strahlungen durch, so erhält man 
eine Schar von Integralkurven; diese zeigen gegenseitige Über- 4 


schneidungen und Wendepunkte, die sich, wie in (I) dargetan wurde, 
ganz systematisch mit der Härte der Strahlung verschieben. Durch 
graphische Differentiation erhält man aus den Integralkurven eine 
Schar von „Differentialkurven“. Jede derselben gibt für die zu- 
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*) Der Radius einer Kammer ist dann „praktisch“ unendlich groß, wenn 7 
eine noch größere Kammer keine größere Ionisation mehr liefert, d. h., wenn eh 
keine Elektronenbahnen mehr die Kammerwandung erreichen und somit kein 
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gehörige Strahlung an, welcher Bruchteil aller Träger zwischen zwei 
koaxialen Nachbarzylindern der Radien o und 9+ 4e pro Ab- 
stand Jo entsteht. Dabei zeigte sich, daß zwei Wolken von Trägern 
vorhanden sind, deren eine, um das Röntgenstrahlenbündel gelagerte, 1 
durch Streuabsorption (also Rückstoß- oder Comptonelektronen), und 
deren andere, weiter außen liegende, durch Photoabsorption (also F 
Photoelektronen) entsteht. Beide Wolken waren entweder räumlich 
völlig getrennt oder sie überlagerten sich so, daß ihre analytische 
Trennung in vielen Fällen möglich war. Dabei zeigte sich hin- 
sichtlich Ausdehnung und Trägerinhalt der Wolken ein systematischer 
Zusammenhang mit der Härte der Röntgenstrahlen, der auch quan- 
titativ so gut mit unseren Kenntnissen übereinstimmt, als erwartet 
werden kann. 

Der kleinste Kammerradius o, für den Untersuchungen vor- 
liegen, beträgt 1 cm. Es wurde nun in (I) angenommen, daß es 
erlaubt sei, mit einem unendlich dünnen Röntgenstrahlenbündel zu 
rechnen und demgemäß die Integralkurven für 9 < 1 so zu extra- 
polieren, daß sie für o = 0 den Wert Null annehmen. Dabei ergab 
sich zwangsläufig die Notwendigkeit, die Integralkurve für die weichste 
Strahlung, bei der die Primärelektronen die kleinste Reichweite haben, 
vom Nullpunkt aus am steilsten, und die Integralkurve für die härteste 
Strahlung mit der größten Reichweite der Primärelektronen unter 
geringster Steilheit ansteigen zu lassen. Dies ließ sich auch recht- 
fertigen durch die mit unseren Erfahrungen im Einklange stehende 
Annahme, daß die Trägerbildung in Kolonnen erfolgt. Daß die 
Annahme eines unendlich dünnen Röntgenstrahlenbündels indessen 
nur eine für die bisherigen Betrachtungen zulässige Näherung dar- 
stellt, zeigt sich besonders deutlich dann, wenn wir einen Schritt 
weitergehen und nach der Dichte der Träger im Abstande o von 
der Achse der Zylinderkammer fragen. Die relative Dichte der 
Träger — relativ deshalb, weil wir uns auf die in einer unendlich 
großen Kammer entstandene Trägerzahl beziehen — erhalten wir 
offenbar dann, wenn wir jeden Wert einer Differentialkurve durch 
das zugehörige 220 dividieren. Für große Werte von o, gegen die 
der wirkliche Radius des Strahlenbündels klein ist, erhalten wir 
dann ein fast völlig richtiges Ergebnis. Je kleiner wir aber o wählen, 
desto größere Werte nimmt die Dichte an, und für o = 0 wird ihr 
Wert unendlich. Das kann natürlich nicht sein und beruht auf d 
unserer für die Dichtebetrachtung nicht mehr zulässigen Annahme, 
daß alle Primärelektronen innerhalb eines unendlich dünnen Röntgen- 
strahlenbündels, d.h. innerhalb eines Volumens Null, ihren Ursprung 
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Für die Sicherung der bisherigen Ergebnisse (I) sowie für An- 
gaben über die Dichte der Träger an den verschiedenen Punkten 
des Raumes scheint es daher notwendig zu sein, die Verteilung der 
Träger innerhalb eines Röntgenstrahlenbündels vom Radius r > 0 
sowie in seiner Umgebung einer näheren Überprüfung zu unterziehen; 
die Rechnung wird, im Rahmen der Zuverlässigkeit unserer ver- 
einfachenden Annahmen, auch Aussagen über die Dichteverteilung 
innerhalb des Bündels und in seiner Umgebung zulassen. 

Als Grundlage unserer Betrachtung dient wieder die Annahme 
der Ionisation in Kolonnen, die mit der Erfahrung im Einklang 
steht. Darüber hinaus müssen wir aber einige vereinfachende An- 
nahmen machen, ohne welche die Rechnung zu verwickelt wird, um 
sich durchführen zu lassen. Und zwar nehmen wir an, die Primär- 
elektronen sollten erstens gleiche Richtungsverteilung, zweitens die- 
selbe Reichweite R, drittens geradlinige Bahn besitzen und viertens 
in jedem Zentimeter ihrer Bahn dieselbe Zahl s von Trägern erzeugen. 
Nun hat zwar Kulenkampff*) überzeugend dargetan, daß einem 
Elektron bestimmter Anfangsgeschwindigkeit in Luft von Atmosphären- 
druck auch eine bestimmte Bahnlänge zukommt; aber die Bahnen 
verlaufen gekrümmt und geknickt und außerdem nimmt die Volt- 
geschwindigkeit des Elektrons bei Erzeugung jedes Trägerpaares um 
einen bestimmten Betrag ab [Kulenkampff°), *) sowie Eisl)], 
so daß gegen das Bahnende zu s immer höhere Werte erreicht. 
Weiterhin besitzen bei einer gegebenen Wellenlänge der Röntgen- 
strahlen die Photoelektronen zwar eine einheitliche Anfangsgeschwin- 
digkeit, aber eine Richtungsverteilung im Sinne einer ,,Voreilung“; 
die Comptonelektronen aber zeigen nicht nur eine ausgeprägte Ver- 
teilung in der Häufigkeit ihrer Emissionsrichtung, sondern auch eine 
ganz ausgeprägte Verteilung ihrer Anfangsgeschwindigkeit in dem 
Sinne, daß die Comptonelektronen mit der größeren Geschwindigkeit 
in Richtung der Röntgenstrahlen, die Elektronen rechtwinklig dazu mit 
der Geschwindigkeit Null emittiert werden. Indessen tritt, durch die 
schnell einsetzenden Bahnkrümmungen, alsbald ein Zustand ein, der 
den anfänglichen Zustand der Richtungsverteilung verwischt und dem 
der allseitig gleichen Richtungsverteilung nahe kommt. Mit dem 
Einfluß dieser Erscheinungen auf das Versuchsergebnis haben wir 
uns bereits in (I) eingehend beschäftigt. Wir müssen uns aber klar 
darüber sein, daß infolge unserer vereinfachenden Annahmen unsere 
Rechnung nur eine Näherung darstellen kann. Immerhin wird sich 
zeigen, daß ihre Ergebnisse in allen grundlegenden Dingen entweder 
mit den Versuchsergebnissen übereinstimmen, oder daß Abweichungen 
zwischen beiden auf unseren vereinfachenden Annahmen beruhen. 
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Die Rechnung wird daher auch in der Lage sein, uns — zum 
wenigsten qualitativ — Fragen zu beantworten, für die die experi- 
mentellen Grundlagen fehlen. Das gilt z. B. für die Verteilung 
der Trägerdichte innerhalb des Röntgenstrahlenbündels und in seiner 
Umgebung. 

Für unsere Rechnung gehen wir von folgendem aus: 

Durch Röntgenstrahlen gegebener Intensität und Härte mögen 
bei Bestrahlung von 1 cm? Luft von Normaldichte Z Primärelektronen- 
bahnen in der Sekunde entstehen, deren jede geradlinig verläuft und 
eine Reichweite von R cm besitzt. Sehen wir von der Luftschwächung 
der Röntgenstrahlen längs ihres Weges durch die Kammer ab [was 
für unsere harten Strahlungen der Untersuchung (1) zulässig ist], so 
werden also durch ein Röntgenstrahlenbündel von 1 cm? Querschnitt, 
dessen wirksame Länge in der Kammer L cm beträgt, insgesamt 
Z.L Primirelektronenbahnen pro Sekunde erzeugt. Wir nehmen 
weiterhin an, das Röntgenstrahlenbündel habe kreisrunden Quer- 
schnitt, sein Radius sei r cm und seine Achse falle mit der Achse 
der betrachteten zylindrischen Räume vom Radius o und der Achse 
der Ionisationskammern zusammen. 


I. Berechnung für den Bereich innerhalb des Röntgenstrahlenbündels 


Wir legen unseren Betrachtungen die Grenzbedingung 
(1) rss 


9 
zugrunde Dabei haben wir zwei Fälle getrennt zu behandeln, 
nämlich erstens Elektronenbahnen, die innerhalb o ihren Ursprung 
haben (Fall A), und zweitens solche, die zwischen o und r entstehen 
(Fall B). 


A. Ursprung der Elektronenbahnen zwischen 9 = 0 und g=r 


In Abb. 1a ist eine Ebene dargestellt, die unsere Zylinder senk- 
recht schneidet; die Achse der Zylinder geht durch M. Innerhalb o 
habe ein Punkt P den Abstand ö von M; für alle Punkte im Ab- 
stande ö von M gelten offenbar gleiche Bedingungen. Zwei zur 
Zylinderachse parallele Ebenen, die gegen MP um die Winkel y 
bzw. w + dy geneigt sind, werden vom Kreise o in den Abständen g, 
bzw. 9, geschnitten. Dann ist offenbar 


9, = + Vo? — sin?yw—d-cosy, 
9, = + Vo? — 0? - sin? w+ 0- cosy. 


Aus dem doppelt keilférmigen Raume, den die Ebenen wy baw. wy+dy 
bilden, ist in Abb. 1b eine Parallelebene zur Zylinderachse dargestellt, 
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die durch P, T und S geht. Die einander parallelen Geraden S’U 
und T’V stellen die Schnitte der zur Zylinderachse parallelen Ebene 
mit dem Zylinder o dar. Schließlich ist in Abb. 1b um P der Kreis 
mit der Reichweite R geschlagen. 

Wir betrachten nun Elektronen, deren Bahnen unter verschiedenen 
Winkeln g gegen die Parallele zu S’U oder T’V durch P in den 
Doppelkeil eintreten. Dabei zeigt sich: ist < 9, bzw. < g,, 80 
verläuft die ganze Bahn innerhalb 9; ist aber 9>g, baw. 9 >, 
so ist nur der Weg w, bzw. w, der Elektronenbahn innerhalb o an : 


Abb. 1. Zur Integration des Bereichs 0 = d= 


der Trägerbildung beteiligt. Wir haben also den Bereich von gy = 0 
bis = bzw. 9, und den Bereich von g = g, bzw. = g, 


bis = gesondert zu behandeln. Es gilt nun offenbar für ¢, 


und g,: 
(3) sing, = und sin = 
1 R 3 R 
Und für. die w, und w, folgt: 


In den Raumwinkel zwischen g und g + dq und zwischen wy und 
yw + dw tritt der Bruchteil 

sing-dg-dy 

4n 
aller Elektronenbahnen. Zwischen zwei koaxialen Zylindern mit 
den Radien d und 6+ d06 entstehen längs eines Strahlwegs von 
L cm in der Sekunde 
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Elektronenbahnen. Erzeugt jedes Primärelektron pro Zentimeter 
Bahnlänge s Träger (eines Vorzeichens), so entstehen längs seiner 
ganzen Bahn 6 


R-s 
Träger, und längs seines Weges w 
w+s 
Träger. 
Integrieren wir von w=0 bis „=, und in jedem Teile 


des Doppelkeiles von g=0 bis =, so erstreckt sich unser 


Ergebnis nur über ?/, der Vollkugel; um diese ganz zu erfassen, 
haben wir also mit dem Faktor 4 zu multiplizieren. 

Die Gesamtzahl aller Träger, die das Röntgenstrahlenbündel der 
Länge L und des Radius r in 1 Sek. erzeugen kann, ist 


rn-L-Z-R:s. 
Dividieren wir hierdurch unser von 0 = 0 bis ö= @ zu erstreckendes 
Integral, so erhalten wir den Bruchteil aller erzeugten Träger, der 
innerhalb des Zylinders o von solchen Elektronen gebildet wird, die 
innerhalb o ihren Ursprung haben. Dieser Bruchteil ist, entsprechend 
dem Vorstehenden, nach Kürzen gleicher Größen 


6=0 y=0 


B. Ursprung der Elektronenbahnen zwischen g und r 


Wir führen nun die Berechnung für solche Primärelektronen 
durch, die ihren Ursprung P zwischen o und r haben und inner- 
halb o Träger bilden. Wir haben hier, wie Abb. 2a zeigt, einen 


-[Vo?—d?-sin?y +0-cosy|-dgdwdo. 
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entsprechenden räumlichen Keil zwischen w und w+ dw, der vom 
Zylinder o bei K und L begrenzt wird. Dabei gilt offenbar: 


\ h, = + Vo? — ö?. sin? y. 


Wie Abb. 2b zeigt, ist hier immer nur ein Teil der ganzen Elektronen- 
bahn wirksam; und zwar zwischen g, und g, das Bahnende, zwischen ¢, 
und 2/2 ein mittlerer Teil der Bahn. Für diese Winkel gilt: 
und sing, = =, | 


j 
Z 


2 / 
h, 
hy 
6 
Abb. 2. Zur Integration des Bereichs =r 
Für die wirksamen Bahnlängen folgt: 3 ‘at * 
‘ i 
(8) w,= k— =R+ = und w,= 


Maximalwert läuft, der durch die Beziehung ute Pi 


0) e 


Wmax = are sin - 3 


bestimmt ist. 

Unter Beriicksichtigung der Ausfiihrungen zum Bereiche A er- 

halten wir dann für den Bruchteil aller erzeugten Träger, der inner- 

halb des Zylinders o von solchen ny EEE gebildet wird, die 

zwischen g und r ihren Ursprung haben: 
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5=r Wwax 9 
J J sing — d-cos y}-dgdwdo 
= y=0 gy; 
F 2 ö- Vo?—0?. sin? 
d= y= Pr 
d=r Wmax 7/2 


C. Die Integration 


Wir müssen zuerst nach integrieren, da die Integrations- 
grenzen wieder yw und ö enthalten. Setzen wir zur Abkürzung 


(11) a= b= ‘cosy, 


so erhalten wir 


y= 


—a.[arc sin (a—b) + arc sin (a+b)] 


to] 


an sin (a + b)] 
— Vl (a +b)}}.dwdö 


bee, 
rd 
2 


| y J 


+b.[aresin(b — a) — arcsin(b + a)] 
+ Y1-(b- a)? — y1— (b-+a)}-dwäd. 


Zur weiteren Integration schreiben wir für den ersten Summanden 
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Hierin ist O=S — 1. Das Integral über w ist also das elliptische 


Normalintegral zweiter Gattung, und der erste Integrand von (12) wird 


20° 0 6 3 


ö 
e 


Entwicklung und Integration nach Gliedern liefert, für das Integral if 


ohne den davorstehenden Faktor die Reihe 


ö\2 1\2 1 1-3\2 1 1 ö \4 
1 2 ö \6 
n \2 5 \4 \6 
deren Koeffizienten die Werte haben: ed 
ur A, = 0,125000, A, = 0,000661, 
A, = 0,015625, A, = 0,000422, 7 
(16) A, = 0,004883, A, = 0,000286, ‘po 
A, = 0,002136, A, = 0,000202, 
A, = 0,001122,  A,, = 0,000148. 


An der unteren Grenze wird das Integral also Null; fir = 1 hat 
die Reihe in der Klammer den Wert C = =, der sich bei 
10 Gliedern bereits bis auf einen Fehler von — 0,81°/,, ergibt. Der 
Beweis für den Grenzwert C darf hier übergangen werden, weil 
sich C, wie wir im folgenden sehen werden, weghebt. 
Somit wird der erste Summand von (12) Rei 


-C. 


n ) 
(17) 2° 
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)] it 
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ER 
Gehen wir beim ersten Summanden von (13) ebenso vor, so wird 
dieser unter Beachtung von | 
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Hier ist aber ie Durch die Substitution 


Die Integrale über x sind die beiden elliptischen Normalintegrale 
erster und zweiter Gattung. Damit wird der erste Summand von (13) 


Indem wir wieder entwickeln, summieren und gliedweise integrieren, 

erhalten wir für (20) nach Einsetzen der Integrationsgrenzen: 
nu r \? o \4 

(21) 6 

- 

wobei die Koeffizienten (16) A,, 4,,... dieselben sind wie in (15) 

und C den oben angegebenen Wert besitzt. Indem wir (17) und 


(21) zusammenfassen, erhalten wir für die ersten Summanden von 
(12) und (13) 


(24) 


: geht (18) über in 
/ — 4 
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Um alle Summanden von (12), die den arcsin enthalten, zu 
integrieren, entwickeln wir die arcsin einzeln als Reihen und fassen 
Summanden, die gleiche Koeffizienten der arc sin-Reihe enthalten, 


zusammen. Damit wird 
— a. fare sin (a — b) + arc sin(a + b)] 
+ b- [are sin (a — b) — arc sin (a + b)] ra 
2{T, - (a? +53) + T, - (at + 6a?b? + bt) 
+ T, +(a® + 15a*b? + 15a? b* + bP) 
+ T,-(a® + 28a°b? + b* + 28a b° + b9) 
+ T,-(a!’ + 45a®b? + 
wobei 
die bekannten Koeffizienten der arcsin-Reihe sind. Da a und b : oa 
in (23) nur mit geradzahligen Exponenten auftreten, so haben die my % 
Integrale über die Einzelsummanden vermöge (11) die Form ; 2 
f sin?” we cos!" yw dy _-* 


und lassen sich dank der günstigen Grenzen sehr einfach lösen; 
jedes Einzelergebnis ist das Produkt aus 2/2 und einem Zahlen- 
koeffizienten. Die Integration von (23) nach w liefert damit 


Gehen wir bei den arcsin-Gliedern in (13) ebenso vor, so wird der rg 
Integrand für diese 

(25) +T,-(4a?b + 28a5b? + 28a? b° + 4a b?) 

+ T,-(5a%b + 60a™b? + 1260555 + 60a°b? + 5ab%)+--.}. 
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Wenden wir die Substitution N 


zZ siny=z dz .cosw-dy =b-dy 


an, so erhalten wir, da hier a und b nur mit ungeraden Exponenten 
auftreten, vermöge (9) und (11) lauter Einzelintegrale der Form q 


die sich dank der günstigen Grenzen sehr einfach lösen lassen; 
die Integration nach ergibt wiederum die Formel (24. Da zu 


den Integralen in (12) und (13) derselbe Faktor = tritt, so 


2 


können wir nunmehr die arc sin-Glieder von (12) und (13) zusammen- 
fassen und nach 6 von 0 bis r durchintegrieren. Das ergibt 


Wir haben nun noch die Wurzeln in (12) und (13) zu integrieren. 
Wir entwickeln wieder und erhalten fiir den entsprechenden Teil 
des Integranden von (12): SE 


2 — yl — (a — by? yı + b)2 a 


+8,-(a + 15atb? + 1502 b* + 
| + S,- (a8 + 28a°b? + 70a*b* + 28a7b* + bs) 
Und für den Integranden mit den Wurzeln in (13) wird: 
yl—(b+a) 

-ab +S,-(2a?b+2ab?) + S,-(3a®b+ 10a?b? + 3a 
| + S,-(4a7b + 28a°b? + 28a°b + 4ab') 

+ S,-(5a%b + 60a?b? + 126a5b® + 600° + 
Dabei sind in (27) und (28) PIE 
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die bekannten Koeffizienten der Reihenentwicklung der Quadrat- 
wurzel. Vergleichen wir nun (27) mit (23) und (28) mit (25), so 
zeigt sich, daß Koeffizienten S und T mit demselben Index auch 
gleiche Klammern als Faktoren haben. Wir können also die 
Wurzelglieder mit den arc sin-Gliedern zusammenfassen. Setzen wir 


2 6 8 24 2 12 

T S => 3 A — 1 1 »i. 

35 7 7 

. . . . . . . . . . . . . - 5 


so erhalten wir mit Hilfe der Einzellösungen für den Bruchteil p, 
aller gebildeten Träger, der innerhalb eines Zylinders mit dem 


P= ?P, + Ps 


II. Bereich außerhalb des Röntgenstrahlenbündels 


Das Ergebnis erhalten wir ohne Schwierigkeit im Anschluß 
an I,A; nur haben wir hier von ö=0 bis d=r zu integrieren. 
Als Grenzbedingung gilt hier an Stelle von (1) offenbar 


(31) o+r=sR; 


dies folgt aus (2) für den größten Weg g, = R, den ein vom Rande 
des Strahlenbiindels (wy = 0 und 6 =r) ausgehendes und die Zylinder- 
achse senkrecht durchsetzendes Primärelektron zurücklegt. Es bleibt 


die Integration nach g unter den alten Bedingungen und, da 21 


ist, auch bei den arcsin und den Wurzeln die alte Integration 


= 


— 


* 


: 
as 
u 
0 
|- 
> 
4 
1 | 
| 
) 
. 
) 


206 Annalen der Physik. 5. Folge. Band 35. 1939 wag are 


bestehen. Wir erhalten fiir das elliptische aidan: mit 
Hilfe von (15): 


@ 
=0 
r)s 
Für den arcsin und die Wurzeln zusammen folgt wegen derselben 
Integrationsgrenzen dasselbe Ergebnis wie bei I,A und I,B zu- 
sammen. 
Damit erhalten wir für den Bruchteil p, aller gebildeten 
Träger, der innerhalb eines Zylinders mit dem Radius o im Be- 
reiche r=o0= R —r entsteht 
1 r 2 r\4 r \6 
r @ 1 2 2 @ r \4 


1 + (2) + (2) + (2) 


Am Rande des Röntgenstrahlenbündels ist o =r; hierfür wird 


offenbar p, = p,; ein Sprung besteht hier also nicht. x 


III. Differentiation von p; und p, nach o 


Die Differentiation der ersten Zeile von (30) nach o gibt unter 
Beriicksichtigung der den Beclänkenten A,, A,, ... nach (16) zu- 
kommenden Werte 


Die Differentiation der ersten Zeile von (32) nach o gibt Ser r 2 
(>) + (2) +54 (2) + 
(34) ‘ 
| 


was sich leicht durch die Reihenentwicklung der elliptischen Normal- 
integrale auf der rechten Seite beweisen läßt. Damit erhalten wir 


(33) a- 
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R'r =. | r 
-|1420(£ +60(£ )’ +40(2)" (a) 
und 
|| 
r » 
An der Grenze des Strahlenbündels, wird sowohl 
i als auch K (=) K (~) E(=)| gleich 1. Die 


Gl. (33) und (34) haben dann denselben Wert, so daß am Rande 
des Strahlenbündels auch kein Sprung des Differentialquotienten 
auftritt. 

IV. Der Bereich R-r=e=R+r ily 

Wie die geometrische Betrachtung lehrt, muß für diesen Be- 
reich eine noch weitgehendere Unterteilung der Winkelgebiete vor- 
genommen werden. Das hat das Auftreten von noch mehr Inte- 
gralen zur Folge. Einfacher gestaltet sich der Sonderfall o = R. 
Die Durchführung der Rechnung ist indessen nicht lohnend, wenn 
man bedenkt, daß außerhalb o = R + r keine Trägerbildung mehr 
stattfinden kann und p, hier den Wert 1 annehmen muß. Im 
allgemeinen kommt man, wie wir sehen werden, bei der rechnerischen 
Auswertung des Bereiches 0= 9 = R —r bereits so nahe an den 
Endwert 1 heran, daß eine graphische Interpolation zwischen R—r 
und R+r völlig ausreicht. Es wurde daher auf die Durchführung 
der Rechnung zwischen R— r und R-+r verzichtet. agents! 
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V. Die Tabulierung der Reihen in (30), (32) und (36) Ir 
1. Die Tabulierung der Reihe 
1— A,-2*— A,-2*— A, 2° — 
läßt sich leicht mit Hilfe der Koeffizienten (16) durchführen. Dies 
geschah in Tab. 1. Für (30) kann man hiernach die Werte von z = © 


und für (32) die Werte = z durch graphische Interpolation ent- 


nehmen. 


Tabelle 1 
Tabulierung der Reihe in (30) und (32) u 
x 1—A,a2*-— A,at*—.-- 


| 


1,00000 
0,99875 
0.99498 
098862 

0,97958 


0.96769 
0,95270 
0,93426 
0,91177 
0,88426 


SOS 
SO © 


.0,84883 = — 
Pr be Die Tabulierung von 
ove 2 , , ay 
in (36) führt man für sehr kleine Werte von r/o am besten mit der 
linken Seite, fiir gréBere Werte von r/o am besten nach der rechten 
Seite von (34) durch. Hierbei wurden die Tafeln der elliptischen 
Normalintegrale von Kiepert‘) benutzt, die 5 Stellen hinter dem 
Komma haben und damit um eine Dezimale weiter gehen als die 
von Jahnke-Emde’). Entsprechend diesen Tabellen ist noch das 


zur Berechnung verwandte @ mitgeteilt, wobei sina = -"- gesetzt 
wurde. Die berechnete Funktion 1) fällt in der Umgebung von 
= 1 schnell ab und konvergiert gut gegen ihren Grenzwert 2/4, 


den sie für — = - schon auf 1°/, erreicht hat. Für die weiteren 


@ 
Berechnungen erforderliche Zwischenwerte wurden durch graphische 

Interpolation gewonnen. Were 
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Tabulierung der Funktion f in (36) 

0 r r 
r | | 
1,000000 | 1,0000 1,000000 
1.000006 0,9997 0,999994 
1,000025 0.99988 0,999975 
1,00005 0,9976 0,99995 
1,00015 0,99940 0,99985 
1,00061 0,99802 0,99939 
1,0014 0,99610 0,9986 
1,0038 0,99108 0,9962 
1,015 0.97438 0,9848 
1,064 0,93476 0,9397 
1,155 0,89591 0,8660 
1305 | 0,86190 0,7660 
1,556 0,83409 0,6428 
2,000 0,81259 0,5000 
. 2,923 0,79739 0,3420 
3,862 0,79219 0,2588 
5,755 0,78856 0,1736 

8,200 0,78688 0,1219 

0,78657 0,0872 

20 0,78567 | 0,050000 

50 0,78544 0,020000 

100 0010000 

0,78540 = 0,000000 
VI. Die Darstellung der Integralkurven = 
In den Gl. (30) für p, und (32) fiir p, treten nur die u re 7 
hiltnisse r/R und o/r auf. 
Ist eine Strahlung bestimmter Härte und damit bestimmter u 

Reichweite R der Primärelektronen gegeben, so kann man R als 

Einheit wählen und r und o als Bruchteile derselben ausdrücken, 

wie das in Abb. 3 in zwei verschiedenen Maßstäben geschehen ist. Be 

Dabei zeigt sich nun, daß sich für jeden von Null verschiedenen my ee 

Radius r des Réntgenstrahlenbiindels die Integralkurven offenbar um di \ a 

so besser der Grenzkurve fiir r = 0 anschmiegen, je kleiner r ist. ee. 

Diese Grenzkurve erhält man leicht, wenn mau in (32) zunächst das r er N 

kürzt und es dann gleich Null setzt. Indem man bei den dann übrig Ta 

bleibenden Gliedern die Koeffizienten wieder gemäß 29) in TundS 

zerlegt, erhält man für die Grenzkurve 

(37) P,|r=0 = R are cos +1 (+) 


Das ist aber die in der Untersuchung (I) den Integralkurven zu- 
grunde gelegte Gl. (1). Am stärksten tritt der Unterschied zwischen 
Annalen der Physik. 5. Folge. 35. 14 
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der Grenzkurve fiir r = 0 und den wirklichen Integralkurven fiir r>0 
innerhalb des Réntgenstrahlenbiindels hervor: hier sinken die letzteren 
mit abnehmendem o mehr und mehr unter die Grenzkurve und be- 
rühren fiir o = 0 die Abszissenachse, während diese von der Grenz- 
kurve unter einem Winkel geschnitten wird. Je kleiner aber r gegen R 
ist, desto weniger fällt der Unterschied ins Gewicht. Dieser Umstand 
rechtfertigt unsere direkte Fortführung der Integralkurven nach dem 
Nullpunkt, wie sie in der Untersuchung (I) durchgeführt wurde. 


0 QR GER QUR Q5R Q6R OBR _Q9R 


S 


GR 
Abb. 3. Integralkurven für verschiedene r als Bruchteile 
derselben Reichweite R, sowie Grenzkurve für r = 0 


Vergleicht man Strahlungen verschiedener Härte und daher ver- 
schiedener Reichweite R miteinander, so kann man r als Einheit 
wählen und R und o als Vielfache (oder als Bruchteil) von r aus- 
drücken (Abb. 4). Diese Darstellung rechtfertigt noch deutlicher das 
in der Arbeit (I) angewandte Vorgehen: im Hinblick darauf, daß es 
sich um ein graphisches Verfahren handelt, ist die Abweichung der 
wahren Integralkurven von den durch den Nullpunkt extrapolierten 
Grenzkurven so gering, daß sie kaum ins Gewicht fallen dürfte. 
Demgegenüber lassen die Integralkurven der Abb. 4 aber sehr ein- 
drucksvoll den Tatbestand erkennen, daß die Integralkurven von der 
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Umgebung des Nullpunktes aus um so steiler ansteigen und um so 


4% Mr kr Mr br br 20r 

Abb. 4. Integralkurven für r = 1 und für verschiedene Reichweiten R & ; 


als Vielfache von r 
A 


\ 


Abb. 5. Integralkurven für eine weiche (I) und eine harte (JZ) Röntgen- 
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weicher also die Strahlung ist. Gerade dieser Umstand bildete eine 
der Grundlagen der Untersuchung (I). 

Dieser Tatbestand legt es nahe, zu prüfen, was sich ergibt, wenn 
man dem Umstande Rechnung trägt, daß in Wirklichkeit kurzbahnige 
Comptonelektronen und langbahnige Photoelektronen gleichzeitig 
vorhanden sind. Es sind demgemäß in Abb. 5 zwei Integralkurven 
dargestellt, die für zwei Röntgenstrahlungen verschiedener Härte 
gelten. Dabei wurde angenommen, bei den Strahlungen J und II 
verteile sich der Anteil an der gesamten Trägerbildung und die 
Reichweite entsprechend Tab. 3 auf Comptonelektronen und Photo- 
elektronen. Die hiermit berechneten Integralkurven zeigen nun 
folgendes: wiederum steigt die Kurve für die weichere Strahlung I 


Tabelle 3 


Beispiel für den Anteil der Compton- und Photoelektronen 
an der Gesamtionisation sowie für deren Reichweiten als Grundlage für Abb. 5 


Anteil an der Gesamtionisation Reichweiten R (als Vielfaches 


bei den von r) bei den 
Strahlung 
Comptonelektr. | Photoelektr. | Comptonelektr. | Photoelektr. 
in | in °/, in °/, in 
I weich 60 40 5 ts 
IT hart 90 10 10 50 


zunächst steiler an. Außerhalb von 6r sind aber keine Compton- 
elektronen mehr wirksam: hier neigt sich die Kurve und strebt nun 
ganz allmählich dem Werte 1 zu, den sie bei 21r erreicht, über 
welche Entfernung hinaus auch keine Photoelektronen mehr gelangen. 
Bei der härteren Strahlung II reichen die Comptonelektronen bis 11 r: 
hier biegt die Kurve JJ um und strebt dem Werte 1 zu, den sie 
(außerhalb der Zeichnung) bei 5lr erreichen würde*). 

Dieser Befund bestätigt die folgenden Punkte, die als Grund- 
lagen der Untersuchung I dienten: 

1. Die Integralkurve der weicheren Strahlung I steigt vom Null- 
punkte aus steiler an als die der härteren Strahlung II; 

2. Die Integralkurve der weicheren Strahlung J sinkt in einem 
Schnittpunkte A unter die der härteren Strahlung II; 

3. Die Integralkurve der weicheren Strahlung I steigt in einem 
zweiten Schnittpunkte B über die der härteren Strahlung II wieder 
empor; 


*) Die Entfernungen, in der eine Elektronenart aufhört, zu ionisieren, er- 


geben sich nach dem unter IV Gesagten zu +n 
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4. Die Integralkurve der weicheren Strahlung I verläuft von 
da ab über der der härteren Strahlung II und erreicht den Wert 1 
eher als letztere. 

Demgegenüber weichen die berechneten Integralkurven in einem 
wesentlichen Punkte von den experimentell gewonnenen ab: die be- 
rechneten Kurven zeigen keine Wendepunkte. Eine Integralkurve 
hat aber dort einen Wendepunkt, wo ihre Differentialkurve ein 
Extremum besitzt. Daß unsere berechneten Differentialkurven außer- 
halb des Röntgenstrahlenbündels aber kein Maximum aufweisen 
können, ist dem Umstande zuzuschreiben, daß wir, wie wir schon 
betonten, einerseits die Zahl s der Träger pro Zentimeter Bahnlänge 
des Primärelektrons als konstant ansahen, und andererseits mit bis 
ans Ende geradlinigen Elektronenbahnen rechneten. Die Träger- 
wolke muß daher mit wachsendem o monoton abnehmen. In Wirk- 
lichkeit nimmt demgegenüber sowohl s als auch die Bahnkrümmung 
vom Ursprung an mehr und mehr zu und bewirkt dadurch das 
Maximum der Trägerwolke, welcher Vorgang bereits sehr eingehend 
in der Untersuchung I besprochen wurde. Da unsere Rechnung, 
wegen unserer vereinfachenden Annahmen, also kein Maximum der 
Differentialkurve liefert, so kann sich auch kein Minimum ergeben, 
das nur aus der Überlagerung der abfallenden Differentialkurve für 
die Comptonträger und aus der zu einem Maximum ansteigenden 
Differentialkurve für die Phototräger entsteht. Das Fehlen der 
Wendepunkte in unseren berechneten Integralkurven ist also nur 
durch unsere vereinfachenden Annahmen verursacht. 


VII. Die Darstellung der Differentialkurven 

In den Gl. (35) für dp,/do und (36) für dp,/do tritt außer den 
Verhältnissen r/R und o/r auch noch R allein auf. 

Ist die Härte einer Strahlung und damit die Reichweite R ge- 
geben und wählen wir R als Einheit, so können wir in (35) und (36) 
das R im Nenner vor der geschweiften Klammer gleich 1 setzen. 
Wir können dann 9 und r wieder als Bruchteile von R ausdrücken 
und dp,/do bzw. dp, /do als Ordinaten auftragen. Wir erhalten so 
die Kurven der Abb. 6. Die Grenzkurve für r= (0 gewinnen wir 
entweder, indem wir (37) nach o differentiieren, oder, indem wir 
in (36) das r kürzen und dann r = 0 setzen. Dabei wird nach Tab. 2 
oder nach (34) der erste Summand in der geschweiften Klammer 2/2; 


die übrigen Summanden liefern die Reihe für — arc sin, so daß 
wir auf beiden Wegen fiir die Grenzkurve das Ergebnis erhalten 


r=0 R R 
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Es ist dies die Differentialkurve, die in der Untersuchung I der 
graphisch gewonnenen Differentialkurve entspricht. 

Wie man erkennt, schmiegen sich die für ein Strahlenbündel 
vom Radius r berechneten Differentialkurven in einigem Abstande 
vom Strahlenbündel der Grenzkurve befriedigend an. Indessen findet 
sich an seiner Grenze und in seinem Innern eine auffällige Er- 
scheinung: dicht außerhalb des Röntgenstrahlenbündels steigt jede 
der Differentialkurven über die Grenzkurve empor, erreicht im Innern 


DER DER DAR D5R 


Abb. 6. Differentialkurven für verschiedene r 
als Bruchteile derselben Reichweite R, sowie Grenzkurve für r = 0 


des Bündels kurz hinter seiner Begrenzung ein Maximum und fällt 
dann zum Nullpunkt ab, den sie unter einem um so größeren 
Neigungswinkel gegen die Abszissenachse schneidet, je kleiner r ist. 

Es soll nun die Lage dieser Maxima und ihr Ordinatenwert be- 
stimmt werden. Hierzu setzen wir den Differentialquotienten von (35) 
nach o gleich Null. Differenzieren des ersten Summanden von (35) gibt 


Entwickeln wir jeden der Summanden in der eckigen Klammer 
rechts, und ziehen wir den zweiten vom ersten ab, so folgt 


| 
] 
I 
| 
> 
| 


der Einsetzen in die vorhergehende Gleichung liefert 
‘ d | 2e@ e\|__2 @ \|. 
de de |rR rR |" \r) r 
sp Differenzieren wir auch die übrigen Glieder, und multiplizieren wir In 


mit rR/2, so erhalten wir 


Indem wir nun für jedes vorgegebene r/R einzeln durch sukzessive 
Approximation ein zugehöriges o/r finden, das (40) genügt, erhalten 
wir die Lage der Maxima und nach Einsetzen in (35) deren Höhe. 
Die Ergebnisse zeigt Tab. 4. 


Tabelle 4 2 
Lage und Höhe der Maxima im Innern des Bündels bei der Differentialkurve 
r | Re | BS | B10 | B20 | 20 0 
| | 0,8635 | 0,8940 | 0,9019 | 0,9055 | 0,9075 | 0,9093 
ap | 1,6422 | 1,9277 | 2,0187 | 2,0644 | 2,0914 | 2,1098 
dg max 


Eigenartig ist dabei der Befund, daß auch beim Grenzüber- 
gange r —> 0 die Grenzkurve, die für g = 0 nach (38) die Ordinaten- 
achse beim Werte 2/2 trifft, auf dieser selbst noch bis zum Werte 2,1098 
ansteigt und dann auf ihr zu Null abfällt. Es ist dies vielleicht ein 
Analogon zur Gibbsschen Erscheinung bei der Fourieranalyse. 

Sollen die Differentialkurven für Strahlungen verschiedener Härte 


. und daher verschiedener Reichweite R bei gleichem Radius r des 
‘ Röntgenstrahlenbündels miteinander verglichen werden, so kann r als 
@ Einheit dienen und R und o können als Vielfache (oder als Bruch- 
5) teil) von r dargestellt werden. Dazu miissen aber (35) und (36) noch 


mit der Einheit r multipliziert werden. In (35) und (36) steht dann 
der Bruch r/R vor g, und beide Gleichungen enthalten dann wiederum 
nur noch die Verhältnisse r/R und o/r. 

Diese Darstellung zeigt Abb. 7. Die ausgezogenen Differential- 
r kurven gelten für den Bündelradius r = 1. Die punktierten Kurven 
sind die Grenzkurven für r = 0 nach (38), die den Differentialkurven ; 
der Untersuchung I entsprechen; für jede punktierte Kurve wurde 


h der Ordinat t fü Od tellt. 
auc er raina engrenzwer arges e 


Ti > 


> = >= 


q H. Kiistner. Über die Raumverteilung der Trägerwolken usw. II 215 Be 
4 
4 
Ba 
Te 
2 ER 
5 
| 
ig 
x 


216 Annalen der Physik. 5. Folge. Band 35. 1939 — 


Fiir die Neigung der Grenzkurve im Nullpunkt erhalten wir ‘ 
durch Differenzieren von (38): 


1 
(+) lo =0 
Je kleiner R, je weicher also die Réntgenstrahlen, desto steiler ist 
die Tangente an die Grenzkurve in der Ordinatenachse nach unten 


tang 7,. 


d [ı 
(41) | are cos, | =— 


r=] 


| [ke 


r Ya br Ar a 


Er Abb. 7. Differentialkurven für r = 1 und für verschiedene Reichweiten R 
18 als Vielfache von r, sowie Grenzkurven für r = 0 


geneigt. Das zeigt auch Abb. 7. Sehr deutlich wird durch die vor- 
stehende Rechnung der Verlauf der Kurven in den Abb. 9 und 10 
der Untersuchung I für die Wolken der Comptonträger bestätigt, 
nämlich: 


1. Je härter die Strahlung und je größer daher die Reichweite R 
ist, desto niedrigere Werte hat die Ordinate in der Ordinatenachse 
und in deren Umgebung; 

2. desto geringere Neigung weist die Kurve in der Odinatenachse 
und in deren Umgebung auf; 

3. desto weiter erstreckt sich die Kurve, bis sie Null wird; und 
4. Die Differentialkurven verschiedener Reichweite R über- 
schneiden sich gegenseitig, und die Schnittpunkte verschieben sich 
systematisch im gleichen Sinne wie bei den Beobachtungen in I. 
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VIII. Die Dichtekurven 
1. Die relative Dichte Tae . Te 


Dividieren wir (35) und (36) durch 220, so erhalten wir die 
„relative Dichte“, nämlich die Zahl der Träger, die in 1 cm’ im Ab- | 
stande g von der Achse gebildet werden, dividiert durch die Zahl h 
aller Träger, die das zylindrische Röntgenstrahlenbündel vom Radius r ‘ 
und 1 cm Lange erzeugt. Dabei zeigt eine einfache Uberlegung, daB 1 
die Länge des Röntgenstrahlenbündels nicht in das Ergebnis ein- 
gehen kann, wenn es nur als hinreichend lang angesehen wird. Es 
werden nämlich zwar die Träger, die durch die Elektronenemission . 
aus jedem Zentimeter Länge des Röntgenstrahlenbündels im Ab- 7 
stande o vom Bündel innerhalb einer Röhre entstehen, die 1 cm? | 
Querschnitt hat und dem Strahlenbiindel parallel läuft, über ein Stück > aw 
dieser Röhre verteilt sein, wobei die Länge, über die sich diese Ver- 
teilung erstreckt, von der Reichweite R der Elektronen abhängt; 
aber jedes Zentimeter Länge des Röntgenstrahlenbündels liefert eine 
solche Verteilung innerhalb der Röhre, und die einzelnen Ver- 
teilungen überlagern sich gegenseitig. Bei hinreichender Länge des er. 
Strahlenbündels kommt es also auf dasselbe hinaus, als wären alle 
Träger, die ein 1 cm langes Stück des Röntgenstrahlenbündels iber- 
haupt in der Röhre erzeugt, in dem diesem Zentimeter zugehörigen, 
1 cm langen Stück der Röhre konzentriert. Es ist also ebenso, als 
erzeugte jedes Zentimeter Röntgenstrahllänge alle Träger im Volumen 
1 cm? und im Abstande von o cm selbst, während in Wirklichkeit 
ein äquivalenter Austausch stattfindet. BE 
Ist die Härte der Röntgenstrahlen und damit die ReichweteR 
gegeben, so können wir wieder R als Einheit wählen. Damit nur en 
die Verhältnisse r/R und o/r in den Gleichungen auftreten, multi 
plizieren wir mit R? und drücken r und o als Bruchteile von R 
aus. Die Darstellung in Abb. 8 zeigt nun sehr deutlich das, worauf 
schon in der Einleitung hingewiesen wurde: die Größe des Radius r 
des Röntgenstrahlenbündels beeinflußt die relative Dichte in der ae" 
Nähe der Achse des Biindels ganz grundlegend; eine willkürliche .; 
Annahme. über r ist hier nicht mehr zulässig. Je kleiner wir r 


annehmen, desto größer wird die Dichte in der Nähe der Achse; = u 
und wird r= 0, so erhalten wir die punktierte Grenzkurve. Sie 4 
ergibt sich mit Hilfe von (36) und (38) zu ech 
R* dp. 1 R 
(42) 2no do R=1 Q 
r=0 


In einigem Abstande vom Röntgenstrahlenbündel schmiegen sich 
alle Kurven der Grenzkurve gut an: wie schon in der Einleitung 
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betont wurde, ergibt sich dort die relative Dichte richtig und zwar d 


unabhängig von der Annahme über den Bündelradius r. Daß die v 


DIR 


Abb. 8. Relative Trägerdichte für verechiedener 
als Bruchteile derselben Reichweite R, sowie Grenzkurve für r = 0 


relative Dichte aber in der Nähe der Achse für r= 0 unendlich 
groß wird, folgt aus der Annahme einer „relativen Dichte“: da wir 
uns dabei auf gleiche Ionisation und damit auf gleichviel gebildete 
Primärelektronen beziehen, so muß mit abnehmendem r die Zahl 
ihrer Ursprungspunkte pro Kubikzentimeter umgekehrt proportional 
dem Quadrate von r zunehmen. 

Für Strahlen verschiedener Härte und damit verschiedener 
Reichweite R, aber bei gleichem Radius r des Strahlenbündels, 
können wir wieder r als Einheit ansehen und R und o als Vielfache 
(oder als Bruchteile) von r ausdrücken. Wir dürfen dann, damit 
nur die Verhältnisse r/R und o/r in den Gleichungen auftreten, 
mit r? multiplizieren und erhalten so die Abb. 9. Hier nimmt die 
relative Dichte, ebenso wie die Höhe der Differentialkurven in 
Abb. 7, in der Umgebung der Bündelachse mit steigendem R ab, 
und die auch bei den Differentialkurven festgestellte Überschneidung 
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der Kurven für verschiedene Reichweiten R wird besonders bei dem 


vergrößerten Maßstabe in Abb. 9b deutlich. 


r 
m 
4 
> 
phy 
? — 7 & 
Ir 5r Wr Ihr 2r ar 
Abb. 9. Relative Trigerdichte für r = 1 und für verschiedene Reichweiten R > 


als Vielfaches von r; bei Multiplikation der Ordinatenwerte mit a erhält man 


die absolute Trägerdichte 


Uri 2. Die absolute Dichte D = 


des Röntgenstrahlenbündels 


(43) D, 20 Z 8 de 2 . R F 
she: . == 


r 


@ 


R 
20 
Um die absolute Dichte aus den Gl. (35) und (36) zu erhalten, = 

haben wir diese nicht allein durch 220 zu dividieren, sondern =. 
wir müssen sie auch noch mit dem Faktor r?2-R-L-Z-s für = , 
L =1cm, d.h. mit der Gesamtzahl aller Träger multiplizieren, die __ 
das Röntgenstrahlenbündel vom Radius r und 1 cm Länge pro Sekunde 


erzeugt. Zu den Gl. (35) und (36) tritt also der Faktor 


Damit erhalten wir für die absolute Dichte innerhalb und außerhalb 


r 
Da Z die Zahl der Primärelektronen pro Kubikzentimeter und pro 
Sekunde, s die Zahl der Träger pro Zentimeter Elektronenbahnlänge 
ist, da o, r und R die Dimension einer Länge haben und da o/r 
und r/R Verhältniszahlen sind, so besitzen die Gl. (43) und (44) die 
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Dimension einer Zahl pro Volumen und pro Sekunde, was mit der 
Definition der Trägerdichte im Sättigungsstrom übereinstimmt. 

Nach unserer in der Einleitung gegebenen Definition stellt das 
Produkt [R-Z-s] die Zahl der Träger dar, die bei Bestrahlung von 
1 cm? Luft pro Sekunde insgesamt gebildet werden. Nun hat 
bekanntlich eine Röntgenstrahlung dann gerade die „Dosis“ von 
„l internat. Röntgen“ oder „Ir“, wenn bei Bestrahlung von 1 cm? 
atmosphärischer Luft von 0°C und 760 mm Hg Druck auf dem 
Wege der Trägerbildung ohne Wandwirkungsausfall und bei Sätti- 
gung eine elektrostatische Einheit der Elektrizitätsmenge entsteht. 
Das Produkt [R-Z.s] ist also proportional der „Dosisleistung“, und 
nach (43) und (44) ist dieser auch die Dichte der Träger in einem 
Punkte des Raumes proportional. Um für den Vergleich der 
Trägerdichte bei Röntgenstrahlen verschiedener Härte ein Maß zu 
haben, wird es zweckmäßig sein, wenn wir unsere Trägerdichte auf 
Röntgenstrahlen gleicher Dosisleistung beziehen. Wir dividieren 
also (43) und (44) durch [R-Z-s] und erhalten 


D; r: dp; r 0 


wobei g, und g, durch (35) und (36) gegeben sind. Die rechten 
Seiten enthalten nun nur noch die Zahlenverhältnisse r/R und g/r. 

Für eine Strahlung bestimmter Härte und damit bestimmter 
Reichweite R sowie bestimmter Intensität oder besser Dosisleistung 
ist R-Z.s konstant. Wählen wir wieder R als Einheit, und drücken 
wir r und o als Bruchteile von R aus, so erhalten wir in Abb. 10 
die Dichtekurven für Strahlenbündel von verschiedenem Radius. 
Man erkennt, daß die Dichte der Träger in der Achse des zylin- 
drischen Strahlenbündels am größten ist und naeh dem Rande des 
Bündels zu abfällt. Am Rande des Bündels hat die Dichtekurve 
offenbar einen Wendepunkt. Außerhalb des Bündels fällt die Dichte 
monoton ab. 

Der Vergleich der Dichtekurven lehrt, daß die Dichte in der 
Mitte des Strahlenbündels nahezu proportional dem Radius r des 
Bündels ist. Außerhalb des Bündels ist die Dichte bei hinreichend 
großem Abstande vom Bündel proportional dem Quadrate des Radius r, 
wie es zu erwarten ist, da die Zahl der Primärelektronen pro- 
portional r? sein muß. Genauer als aus Abb. 10 läßt sich das aus 
Tab. 5 entnehmen, die die berechneten Zahlenwerte für den Fall 


0 -2 selbst wieder gibt; immerhin zeigt aber auch Abb. 10, dab 
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bereits für 9 =2r die Proportionalität der Dichte mit r? in jedem 
Falle erfüllt ist. 

Ist bei gleichem Radius r des Strahlenbündels die Härte der 
Röntgenstrahlen und damit die Reichweite R der Elektronen ver- 
schieden, die Dosisleistung aber dieselbe, so können wir R und o 
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Abb. 10. Absolute Trägerdichte für verschiedener 
als Bruchteile derselben Reichweite R 
a 
Tabelle 5 

Proportionalität von Trägerdichte und Querschnitt des Röntgenstrahlenbündels ur 
in großem Abstande von diesem, berechnet für halbe Reichweite «3 
Dichte, berechnet Dichteverhältnis ) 2 [7 ) 2 4 
1 
0,01371 | 
A > 
1 4,04 10 A 
0,000827 
6,20 6,25 

0,0001332 
50 | 
als Vielfache (oder als Bruchteile) von r darstellen. Wir erhalten Pa 
dann offenbar Kurven, die mit denen der Abb. 9 zusammenfallen, eo 
wenn wir die Bezifferung der Ordinaten in Abb. 9 mit a multi- pad 


plizieren. Da sich also nur der Ordinatenmaßstab ändert, so wurde 
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auf eine neue Darstellung verzichtet. Man erkennt, daß bei gleichem 
Radius r des Röntgenstrahlenbündels und bei gleicher Dosisleistung 
für Strahlen verschiedener Härte die Dichie der Träger in der Achse 
des Strahlenbündels nahezu umgekehrt proportional der Reichweite 
der Elektronen ist. Dieses Ergebnis ist qualitativ einleuchtend, da 
sich bei weichen Röntgenstrahlen die Bildung derselben Trägerzahl 
innerhalb eines engeren Raumes vollzieht als bei harten Strahlen; 
quantitativ ließ es sich indessen kaum voraussehen. 

Man könnte geneigt sein, hieraus den Schluß zu ziehen, daß 
bei einer weichen Strahlung in einer Zylinderkammer schwerer 
Sättigung zu erzielen sei als bei harter Strahlung. Dem ist aber 
entgegenzuhalten, daß sich mit zunehmender Strahlenhärte die 
Trägerbildung von der Achse der Kammer mehr und mehr nach 
der Wandung hin und damit in Gebiete verschiebt, in denen die 
Feldstärke viel geringer ist als in der Umgebung der dem Bündel 
meist benachbarten Elektrode. ; 


IX. Die experimentellen Ergebnisse für die Trägerdichte 

Ein einfach geladener Träger hat die Ladung des Elektrons 
von 4,770.1071° elektrostatischen Einheiten. Hat eine Röntgen- 
strahlung die Dosisleistung von 1 „r“ pro Sekunde, so werden nach 
unserer unter VIII, 2 aufgeführten Definition des „r“ bei Be- 
strahlung von 1 cm? Luft von 0°C und 760 mm Hg Druck in 
jeder Sekunde 1/4,770 - 10*1° = 2,095 . 10° Träger befreit und im 
Sättigungsstrome gemessen. Bei {°C und H mm Hg Druck entsteht 

1 H 


von diesen nur der Bruchteil 


für den sich bei den Messungen ?) (+ 23°C, H = 750 mm Hg) der 
Wert 0,91 ergab. Beträgt die Dosisleistung X „r“ pro Sekunde, 
und hat das Röntgenstrahlenbündel einen Radius von r cm, so ent- 
stehen bei 1 cm Länge des Strahlenbündels 0,91-2,095-10°.r? a-X 
Träger in der Sekunde. Nun ist nach VIII,1 a ee i 

2no de 

der Bruchteil aller durch ein Röntgenstrahlenbündel von 1 cm 
Radius und von 1 cm Länge gebildeten Träger, der im Abstande o 
von der Bündelachse in jedem Kubikzentimeter entsteht. Das Produkt 


(47) N =0,91-2,095-10%. . Pe — 0,952: 109. 7. 
2ng dg e de 
ist demnach die Zahl der Träger, die ein Röntgenstrahlenbündel vom 


Radius r cm und der Dosisleistung X ,,r“ pro Sekunde im Abstande 
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o cm von seiner Achse im Kubikzentimeter während einer Sekunde 
bei 23°C und 750 mm Hg Druck erzeugt. Hierin hat nach (36) 

r 

eg gR 7) 
die Dimension einer Zahl; 0,952-10°® hat die Dimension einer 
Trägerzahl pro Ladung und X hat die Dimension einer Ladung 
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Abstande g von der Achse des Röntgenstrahlenbündels von 1 cm Radius bei 
1 „r“ pro Sekunde (sowie 23°C und 750 mm Hg) (Hauptwolke) 


pro Kubikzentimeter und pro Sekunde, so daß (47) die Dimension = 


einer Trägerzahl pro Kubikzentimeter pro Sekunde hat. * * 

Der Wert von dp,/do läßt sich nun den Differentialkurven der u +f 

Untersuchung I entnehmen. Dividieren wir ihn durch o und multi- 

plizieren wir ihn mit 0,952 - 10%, so erhalten wir 
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d.h. die Zahl der Träger im Kubikzentimeter, die bei 23°C und 
750 mm Hg Druck im Abstande @ von der Achse des Röntgenstrahlen- 
bündels während jeder Sekunde entstehen, wenn dessen Radius 1 cm 
beträgt und wenn die ,,Dosisleistung pro Sekunde 1 ,,r ist. Das y 
ist in den Abb. 11—13 in Abhängigkeit von o dargestellt; aus 4 


(49) N=y-r?-X 


folgt fiir den Biindelradius yon rcm und die Dosisleistung X in 
„r“ pro Sekunde die Trägerdichte N (47). Da hier das N dem r? 


Abb. 12. Zahl der Phototräger pro Kubikzentimeter und pro Sekunde im _ 
Abstande og von der Achse des Röntgenstrahlenbündels von 1 em Radius bei 
1 „r“ pro Sekunde (sowie 23°C und 750 mm Hg) (Schleierwolke) 


proportional gesetzt ist, so folgt aus unserem Befunde von VIII,2, 
daß die Abb. 11—13 nur für 9 =2r Geltung haben. 

Vergleichen wir die Dichte der Hauptwolke der Phototräger 
(Abb. 11) mit der Abb. 7 der Untersuchung I, aus der sie gewonnen 
wurde, so zeigt sich bei 60 und 80 kV eine geringe Verschiebung 
des Maximums auf das Röntgenstrahlenbündel zu. Die Werte für 
o< 1 cm sind nicht so zuverlässig, was auch schon aus der Streuung 
der Punkte hervorgeht. 

Die Dichte im Wolkenschleier der Phototräger gibt Abb. 12. 
Der Ordinatenmaßstab ist hier 100mal so groß gewählt wie in 
Abb. 11. Die systematische Überschneidung der Einzelkurven tritt 
hier wieder sehr deutlich zutage. j 
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Auf die Dichtedarstellung der Wolke der Comptonträger von N 
60—195 kV wurde verzichtet, da die Wolke im Röntgenstrahlen- Si 
bündel selbst ihr Maximum hat und schon in so geringem Abstande 
vom Strahlenbündel abgeklungen ist, daß die Forderung 9 =2r 
kaum erfüllt werden kann. 

Demgegenüber hat die Wolke der Comptonträger von 250 bis 
380 kV eine solche Ausdehnung, daß ihre Dichte in Abb. 13 wieder- 
gegeben werden konnte. Dabei wurde absichtlich nicht weiter als 


A 


ting 

( 

a 
Abb. 13. Zahl der Comptonträger pro Kubikzentimeter und pro Sekunde im : : 
Abstande g von der Achse des Röntgenstrahlenbündels von 1 cm Radius bei a 
1 „r“ pro Sekunde (sowie 23°C und 750 mm Hg) ,o tf 

bis auf 1,5 cm an die Bündelachse herangegangen. Wir erkennen oe 7 
hier eine vorzügliche Übereinstimmung im Charakter der Kurven = 
mit den berechneten Kurven der Abb. 7, und zwar sowohl hinsicht- = 
lich der gegenseitigen Überschneidungen (in Abb. 13 oben durch “4 
Pfeile markiert) als auch bezüglich der Reihenfolge bei der Ein- , 
miindung in die Abszissenachse. 
Daß die Kurven der Abb. 7 in die Abszissenachse unter einem = 
Winkel einmünden, während sie sich beim Versuch (Abb. 13) tan- 
gential an diese heranziehen, erklärt sich zum kleineren Teil durch 7 2 
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die in der Rechnung vernachlässigte Richtungsverteilung der Primär- 
elektronen, vorwiegend aber durch die Inhomogenität der benutzten 
Strahlengemische, die Elektronen sehr verschiedener Reichweiten 
erzeugen. 

Für eine andere Luftdichte als die bei den Versuchen, für die 
auch die Abb. 11—13 gelten, nämlich 23°C und 750 mm Hg Druck, 
erhält man die Trägerdichte leicht durch folgende Überlegung. 
Ändert man weder den Bündelradius noch die Dosisleistung der 
-Réntgenstrahlen, wählt man die Luftdichte aber kmal so groß wie 
beim Versuch, so geht das räumliche Bild der Trägerwolke auf den 
Maßstab 1:k zurück. Dann werden also die Abszissen der Abb. 11 


bis 13 übergehen in 3 


00) 
Die gleiche Zahl von Trägern, die vorher in 1 cm? entstanden, wird 
bei der kfachen Dichte in 1/k? cm? gebildet werden; die Zahl der 
Träger im Kubikzentimeter ist dann also k*mal so groß geworden, 
die Ordinaten der Abb. 11—13 gehen über in 


(51) = y-k’, 
und Gl. (49) nimmt die Form an 


(52) 

Ist die Dichte bei t? C und bei H mm Hg Druck kmal so 
groB wie bei den Versuchen des Verf., für die die Abb. 11—13 
gelten, so ist 


(8) k= 1 + H . H 
~ 1+ 0,00367-¢ 750° + 0,00367 - ¢ 
Zusammenfassung 


1. Ausgehend von der Tatsache, daß die Trägerbildung durch 
Röntgenstrahlen nach Kolonnen erfolgt, wird unter vereinfachenden 
Annahmen die Verteilung der Träger innerhalb und außerhalb eines 
Röntgenstrahlenbündels von kreisrundem Querschnitt in Abhängig- 
keit von dessen Radius r für jeden Raumpunkt im Abstande o von 
der Bündelachse berechnet. 

2. Die „Integralkurven“ und die „Differentialkurven“, die in der 
früheren Untersuchung I des Verf. die Trennung der Trägerwolken 
aus Streuabsorption und Photoabsorption für ein unendlich dünnes 
Röntgenstrahlenbündel auf experimenteller Grundlage ermöglichten, 
werden für ein kreisrundes Röntgenstrahlenbündel vom Radius r 
berechnet. Die früheren Kurven ergeben sich durch den Grenz- 
übergang r —> 0. Die Unterschiede zwischen r=0 und r>0 
werden erörtert. REN 
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3. Die für die Trennung der beiden Trägerwolkenarten grund- 
legenden Überschneidungen von Integralkurven, die für Röntgen- 
strahlen verschiedener Härte auf experimenteller Grundlage gefunden 
wurden, werden durch die Rechnung sichergestellt. 

4. Die Dichte der Trägerwolken zeigt der Rechnung nach auf 
der Achse des Röntgenstrahlenbündels ein Maximum und fällt nach 
dem Rande hin ab; beim Überschreiten des Bündelrandes besteht 
kein Dichtesprung. 

5. Für eine Strahlung gegebener Härte und Intensität ist die 
Trägerdichte in der Achse des Röntgenstrahlenbündels der Näherungs- 
rechnung nach nahezu proportional dem Radius r des Strahlen- 
bündels; außerhalb des Bündels ist die Dichte um so genauer pro- 
portional r?, je größer der Abstand o vom Bündel ist. 

6. Bei gleichem Radius des Strahlenbündels und bei Röntgen- 
strahlen verschiedener Härte, aber gleicher Dosisleistung ist die 
Trägerdichte in der Achse des Bündels der Näherungsrechnung nach 
nahezu umgekehrt proportional der Reichweite R der Primär- 
elektronen. 

7. Im Anschluß an die früheren Ergebnisse '), *) des Verf. 
wird für die untersuchten extrem hart gefilterten Einzelstrahlungen 
zwischen 60 und 380 kV zahlenmäßig angegeben, wieviel Träger 
pro Sekunde und pro Kubikzentimeter in einem beliebigen Raum- 
punkte durch Photoelektronen oder durch Comptonelektronen erzeugt 
werden, wenn der Radius des Röntgenstrahlenbündels und die 
Dosisleistung der Strahlung nach internationalen „r“-Einheiten pro 
Sekunde bekannt sind. 


5) A. Eisl, Ann. d. Phys. [5] 3. S. 277. 1929. 
6) L. Kiepert, Grundriß der Differential- und Integralrechnung, II., 
S. 681 und 682, IX. Aufl., Helwigsche Verlagsbuchhandlung, Hannover 1908. 
7) E. Jahnke u. F. Emde, Funktionentafeln, 2. Aufl., S.150, B. G. Teub- 
ner, Leipzig 1933. 
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4) H. Kulenkampff, Ann. d. Phys. [4] 79. S. 97. 1926. “ 


A 


4 


as... - 
Pax 


n 
n 
ar Br 
ie 
| 
n 
1 ph 
d | 
af 2 
0 
— 
1 
1 
3 
| 
1 
1 
3 
KURSE 
x 


der Ph 


5. Folge. Band 35. 1939 ae 


Zur optischer Wellen, 
insbesondere von Kugelwellen 


Von H.W. Breuninger’) 


(Mit 6 Abbildungen) 


Pflanzt sich eine unbegrenzte ebene Welle von einem Punkte P, 
nach einem Punkte P, fort, wobei die Verbindungsgerade von = 
mit P, auf der W ellenfläche senkrecht stehen möge, und ist (P 


die Phase der Lichtschwingung in P, dann ist 
= 


1 
wobei k = en 4 die Wellenlänge und P,P, die Entfernung der 


beiden Punkte bedeuten. Im Fall einer BE sich von P, nach 
P, fortpflanzenden Welle ist 


w(P,) — w(P,)=k-P,P,— 3. 


ö ist die Phasenanomalie. Bei positivem 3 ist die Phase in P, 
größer im Vergleich zur Phase in P, als bei einer ebenen Welle: 
Positivesöbedeuteteine Phasenbeschleunigung, negatives d entsprechend 
eine Phasenverzögerung. 

1890 hat Gouy?) theoretisch gezeigt, daß eine konvergente 
Kugelwelle beim Durchgang durch ihren Brennpunkt eine Phasen- 
beschleunigung a erleidet. Die experimentelle Prüfung erfolgte mit 
dem Fresnelschen Spiegelversuch unter Verwendung eines Plan- 
und eines Hohlspiegels; die durch Reflexion am Planspiegel entstehende 
Welle divergiert vom Spiegelbild der Lichtquelle aus; durch den 
Hohlspiegel entsteht eine Kugelwelle, die gegen ihren Brennpunkt’) 
konvergiert und dann divergiert. Hinter diesem Brennpunkt war 


1) Mit Ergänzungen versehener Auszug der Jenaer Dissertation „Zur 
Theorie des Verhaltens der Phase konvergenter Kugelwellen, insbesondere in 
der Umgebung des Brennpunktes“ re R. Noske, Borna, 1938). 

2) G. Gouy, Compt. rend. 110. S. 1251. 1890; 11. 8. 33, 910. 1890; Ann. 
chim. phys. (6) 24. S. 145. 1891. 

3) Unter Brennpunkt der Kugelwelle wird hier der Punkt verstanden, 
nach welchem die Kugelwelle konvergiert. Er ist im allgemeinen verschieden 
vom Brennpunkt des Hohlspiegels. In obigem Fall ist er der durch den Hohl- 
spiegel erzeugte Bildpunkt der Lichtquelle. Nur bei deren unendlicher Ent- 
fernung fallen Brennpunkt der Kugelwelle und des Hohlspiegels zusammen. 
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4 5) F. Reiche, Ann. d. Phys. [4] 29. S. 65. 1909. 


H. W. Breuninger. Zur Phasenanomalie optischer Wellen usw. 229 


der mittlere Interferenzstreifen dunkel, wie Gouy zur Bestätigung 
seiner Theorie fand. Später wiesen Joubin!) und Fabry?) diese 
Phasenanomalie 2 mit den Newtonschen Interferenzringen nach. 
Eine andere Versuchsanordnung stammt von Zeeman’), die auch 
Sagnac*) und Reiche’) benutzt haben: Das durch einen Spalt 
eintretende Licht passiert ein polarisierendes Nikol, dann eine durch 
eine kreisförmige Blende begrenzte, plankonvexe Kalkspatlinse und 
schließlich einen Analysator; die Nikols sind parallel oder gekreuzt; 
die optische Achse der Linse steht senkrecht zur Achse der Versuchs- 
anordnung und bildet einen Winkel von 45° mit den Hauptschnitten 
der Nikols. Durch die doppelbrechende Linse entstehen zwei Kugel- 
wellen mit verschiedenen Brennpunkten; diese Wellen besitzen 
gegeneinander eine Phasendifferenz, die durch Drehung einer zwischen 
der Linse und dem Analysator angebrachten planparallelen Kalkspat- 
platte beliebig eingestellt werden kann. Zeeman zeigte, daß wenn 
die Mitte der Interferenzerscheinung der beiden Wellen vor dem 
ersten Brennpunkt dunkel war, sie zwischen den beiden Brennpunkten 
hell und hinter dem zweiten Brennpunkt wieder dunkel wurde. 
Sagnac, der mit kleinerer Blende und größeren Brennweiten 
arbeitete, fand, daß bei kleiner-Blendenöffnung am Brennpunkt kein 
plötzlicher Umschlag der Mitte des Interferenzbildes von hell nach 
dunkel erfolgt, sondern daß ein solcher Umschlag an mehreren 
Stellen zwischen den Brennpunkten stattfindet; die Phasenanomalie 
schwankt in Brennpunktsnähe oszillatorisch; erst bei größerer Blenden- 
öffnung tritt der Helligkeitsunterschied bereits an Punkten kurz vor 
und hinter dem Brennpunkt auf. Die Erklärung der Phasenanomalie 
gab Sagnac auf Grund der beugungstheoretischen Behandlung einer 
durch eine kreisförmige Blende ausgeblendeten Kugelwelle und wies 
auf die Bedeutung hin, die unter Umständen bei den Experimenten 
den Unvollkommenheiten der Versuchsanordnung zukommen kann, 
indem dadurch die Feinheiten verwaschen werden. Etwas genauer 
als Sagnac berechnete Reiche das Verhalten einer ausgeblendeten 
Kugelwelle längs der optischen Symmetrieachse; die seiner Arbeit 
entnommene Abb. 1 stellt sein Ergebnis dar: Weit vor dem Brenn- 
punkt ist die Phasenanomalie ö= 0; dann fängt 0 an, zwischen 


1) P. Joubin, Compt. rend. 115. S. 932. 1892. 
x 2) C.Fabry, Compt.rend.115. S.1063. 1892; Journ. dephyes. (3) 2. S.22. 1893. 
3) P.Zeeman, Phys. Ztschr.1. 8.542. 1900; Arch. Néerland. (2)4. S. 318. 1901. 
4) G.Sagnae, Journ. de phys. (4) 2. $. 721. 1903; Compt. rend. 138. S. 479. 
1904; Compt. rend. 138. S. 619. 1904; Compt. rend. 138. S. 678. 1904; Compt. 
rend. 139. S. 186. 1904; Boltzmann-Festschrift, S. 528. 1904. 
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positiven und negativen Werten zu pendeln; die Amplituden dieser 
Oszillation wachsen gegen den Brennpunkt zu; an diesem ist 


ö=-, und schließlich pendelt ö hinter dem Brennpunkt um gz, 


wobei mit zunehmender Entfernung vom Brennpunkt die Größe der 
Schwankung abnimmt, so daß ö weit hinter dem Brennpunkt gleich 
a wird. Bei der Berechnung hat Reiche vorausgesetzt, daß die 
Brennweite f und die Entfernung p des auf der Achse gelegenen 


A 


es Abb. 1. Aus Ann. d. Phys. [4] 29. S. 82. 1909 


Aufpunktes P von der Blende groß gegen deren Radius R seien: 
f>R, p>R. Für den Fall, daß auf die Blende eine ebene Welle 
auftrifft, und diese durch die dort befindliche Linse in eine nach 
dem Brennpunkt konvergente Kugelwelle verwandelt wird, gelten 
die Formeln: 

(1) 


mit der Nebenbedingung, daß cosd und 


2-2 - Q,%)cos@ 
dasselbe Vorzeichen haben. Dabei ist 


4 
] 
] 
( 
| ! 
i” % 0 if ? 
| | | 
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= ] 
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| 


lle 


H. W. Breuninger. Zur Phasenanomalie optischer Wellen usw. 231 Gr dee 
Entsprechend wird die Intensität 
A? 

I = |Up|? = (f 4 (2 — Q,”) cosQ}, 
wobei up die Lichtschwingung in P in komplexer Schreibweise ist; 


A ist die Amplitude in der Blende. Wegen der Voraussetzungen 
R<f, R<p ist Q,?<2; vernachlässigt man daher in obigen 


Formeln Q,? gegen 2, so wird Pee | ch 

wobei das obere Vorzeichen vor und das untere hinter dem Brenn- = i 
punkt gilt und [x] die größte ganze Zahl kleiner oder gleich z . Br 


bezeichnet. Die Beugungserscheinung, die längs der optischen Achse 7 
auftritt, wenn eine ebene Welle eine Kreisöffnung, in der sich keine = 
Linse befindet, passiert, ergibt sich aus diesen Formeln fiir f = oo. 

Mébius’) hat die Phasenanomalie für Wellen berechnet, deren 
Wellenfläche eine in bezug auf die Ausbreitungsrichtung symmetrische 
Sattel-, Zylinder- oder Rotationsfläche ist. Weiter unten werde ich 
auf seine Ergebnisse noch etwas eingehen und sie ergänzen. 

Debye?) gelang es, das Problem der durch eine Öffnung aus- 
geblendeten konvergenten Kugelwelle dadurch einfach zu gestalten, 
daß er die Blende nach Unendlich verlegte. Dabei ergab sich: Bei 
beliebiger Form der Blendenöffnung und überall in der Öffnung 
gleicher Amplitude ist 


(7) Up= f da, ‘= 

wobei der Aufpunkt P irgendwo im Raum liegen kann; r ist die a 

Entfernung des Aufpunktes vom Brennpunkt, der als Koordinaten- : 
a 


anfang gewählt ist; r ist der Winkel zwischen dem Radiusvektor —_ 
von P und dem Strahl, der in dem Raumwinkelelement dQ liegt. on 

(7) ist eine strenge Lösung von 4up+k?up=0 trotz dr = 
Kirchhoffschen Annahmen, in welchen die Unstetigkeit steckt, daB a 
in der Blende die Lichtschwingung gleich ist wie bei ungestörter ig \ 


Ausbreitung, daß sie aber an allen Stellen hinter dem Schirm ver- . BF 
schwindet; am Blendenrand wird also eine Unstetigkeit vorausgesetzt. = 
Wenn ich nun die Blende so nach Unendlich rücken lasse, daß der u _ 


durch sie ausgeblendete Lichtkegel konstanten Offnungswinkel behält, a we 


1) W. Mobius, Ann. d. Phys. [4] 33. 8. 79. 1910. Su ‘mes 
2) P. Debye, Ann. d. Phys. [4] 30. S. 755. 1909. 
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dann nimmt die Fläche der Blendenöffnung wie r?, der Rand aber 
wie r zu, so daß der Einfluß des Randes im Verhältnis zur Fläche 
verschwindet. Dies erklärt anschaulich, wie trotz der am Rand an- 
genommenen Unstetigkeit (7) eine strenge Lösung ist. 

Für eine kreisförmige Blende folgt aus (7) bei Einführung von 
Polarkoordinaten um die Achse des Lichtkegels vom halben Öffnungs- 


winkel « | 

(8) Up = ik (kr sin sin #)-sin 

Die Koordinaten von P sind dabei r,d,9. #=0 kennzeichnet 

axiale Aufpunkte vor und «+= solche hinter dem Brennpunkt. 

Auf der Achse wird: 


9) 


wobei das obere Vorzeichen für # = 0, das untere fir J = 7 gilt. 
Für nichtaxiale, weit vom Brennpunkt entfernte Aufpunkte hat 
Debye näherungsweise up berechnet. Er teilte dazu das Integrations- 
gebiet in zwei geeignete Bereiche und stellte so up als Summe 
zweier Bestandteile, up; und upy, dar. wp; konnte er genau berechnen; 
für upır fand er einen Näherungswert dadurch, daß er wpy durch 
partielle Integration des betreffenden Integrals in eine Reihe nach 
r-! entwickelte und dann nur das erste Glied berücksichtigte. So 
fand er für weit vom Brennpunkt entfernte nichtaxiale Aufpunkte 


ikr 
(10) Up= — - , wenn 0< ta, did: 
Up= 0 wenn 
-ikr 
(11) Up=— ‚wenn <I < ist. 


Aus dem Vorzeichenunterschied zwischen (10) und (11) folgt, daß 
zwischen weit hinter dem Brennpunkt liegenden nichtaxialen Auf- 
punkten und weit vor dem Brennpunkt liegenden eine Phasen- 
anomalie x besteht. Ebenso haben v. Laue?) und Jentzsch?) aus 
(9) geschlossen, daß zwischen zwei axialen Aufpunkten, von denen 
der eine vor und der andere hinter dem Brennpunkt liegt, die 
Phasenanomalie x auftrete. Picht*%) hat auf diesen Irrtum aufmerksam 
gemacht. Ist die Phasenanomalie des Achsenpunktes vor dem Brenn- 
punkt ö,, diejenige des hinter dem Brennpunkt gelegenen Achsen- 


1) M. v. Laue, Enzykl. d. math. Wiss. 5. Teil 3. S. 439—444. 
2) F. Jentzsch, Handb. d. Physik 18. S. 265— 268. mE : 
3) J. Picht, Ztschr. f. Phys. 65. S. 14. 1930. Kate 
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ber punktes ö,, dann ist die Phasenanomalie zwischen den beiden 

che Punkten ö, — ö,, und diese Differenz kann bei noch so großer Ent- 

an- fernung der Aufpunkte vom Brennpunkt alle Werte zwischen 0 und 
2a annehmen. 

on Unlängst hat sich der Verfasser!) mit diesem Problemkreis 

gs. auseinandergesetzt. Es ergab sich: 


L I. Wenn wie bei der Berechnungsmethode von Reiche die 
‘ Blende im Endlichen liegt, dann schrumpfen die Amplituden der 
q Schwankungen von 6 nicht so rasch mit zunehmender Entfernung 

vom Brennpunkt zusammen, als die Abb. 1 den Anschein erweckt. 
ret Für das Gebiet zwischen Brennpunkt und Blende konnte durch 
kt. graphische Veranschaulichung des in der Kirchhoffschen Formel 
auftretenden Integrals gezeigt werden, daß bei einer ebenen durch 
eine Kreisblende tretenden Welle und ebenso bei einer konvergenten 
Kugelwelle, wenn die Brennweite f> R (Blendenradius) ist, in axialen 


lt. Aufpunkten in der Nähe der Blende ö noch mindestens zwischen 
at > 5 mit sich verändernder Entfernung von der Blende schwankt. 
ni Ein Erlöschen der Oszillationen von ö ergibt sich selbst aus (1) 
n: nicht, solange (1) überhaupt gilt, d. h. solange p>R ist. 
ch II. Hinter dem Brennpunkt finden bei im Endlichen gelegener 
ch Blende nur endlich viele, nämlich Pea Schwankungen von 6 statt. 
Ka Dabei sind die Amplituden der Pendelung von ö hinter dem Brenn- 
. punkt stets so groß, daß ö ein Intervall von mindestens 140° durch- 
läuft. Nur der Mittelwert, um den 0 schwankt, ‘unterscheidet sich 
hinter dem Brennpunkt um z von dem Mittelwert vor dem Brennpunkt. 
Selbst bei großer Brennweite f braucht das Verhalten von 0 sich 
noch lange nicht dem Debyeschen Fall (Verlegung der Blende ins 
Unendliche) zu nähern. 
B III. Nach (1) bis (4) kann 0 als Funktion von g aufgefaßt werden. 
f- Da g vor dem Brennpunkt negativ, hinter ihm positiv ist, ergibt 
1- sich das Verhalten von ö hinter diesem aus demjenigen vor ihm 
1S durch Vorzeichenumkehr bei g, wobei zu beachten ist, daB vor dem 
n Brennpunkt g alle Werte von — oo bis0 annimmt, hinter ihm jedoch 
e nur von Obis 1/f zunimmt; dies scheint Reiche bei Abb. 1 nicht 
n beachtet zu haben. 


1) H. W. Breuninger, Jenaer Dissertation 1938. 
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wobei das obere a für#= 0, das untere für # = 2 gilt. 
Das hierdurch gegebene 


a der Achse wird durch 
21 A A V. Führt man in (1) 

0 bis (6) den Grenziiber- 


| unter den Bedingungen 

f—p| A-f=1 
ae tg« endlich bleiben, dann entstehen Formeln, 


f 
die für kleine «, was den Voraussetzungen R<p, R< f entspricht, 
mit (9) übereinstimmen. 

VI. Es wird gezeigt, daß die Berechnung von up für nicht- 
axiale weit vom Brennpunkt entfernte Aufpunkte, wie sie Debye 
durchgeführt hat, mathematisch nicht einwandfrei ist, indem bei 
der Zerlegung von up in up und upy die Reihenentwicklung von 
Upy nach r—! nicht zugelassen werden darf, weil der Koeffizient 
von r”? einen unendlichen Wert annimmt. Daher wird eine neue 
Berechnung dieses Falles durchgeführt, die darauf beruht, daß 1. für 
die Besselsche Funktion im Integranden von (8) für große Argument- 
werte die Näherungsformel 


und lim = lim 


a 
Jy (0) = +008 (2 4) 


eingesetzt wird, und daß 2. der Integrationsbereich von 0 bis « in Teil- 
strecken eingeteilt wird, z. B. im Falles" (s=kr) in 
a) Obiss-", b)s"bs®—e, c) d) #bis# + ¢, 
e) *#+ ebise, wobei e = 2arcsins~‘ ist. Es wird dann gezeigt, 
daß die Teilintegrale entsprechend den Teilstrecken a), b), e) rascher 
als s~! (aber vermutlich weniger rasch als s~*) für s-> oo ver- 
schwinden, während die Teilintegrale entsprechend c) und d) wie 
s~! nach Null gehen. So ergibt sich für « < 7!) dasselbe wie bei 
Debye. 


1) H.W.Breuninger a.a.0., 8.20 wurde vorausgesetzt cos e>s” '/*. Diese 
Voraussetzung selbst wurde nur auf S. 29 bei dem Beweis gebraucht, daß im 
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VII. Es wird für den Fall einer im Unendlichen liegenden 
Kreisblende, durch die eine konvergente Kugelwelle eintritt, gezeigt, 
daß wenn die Aufpunkte P, und P, in derselben Entfernung vom 
Brennpunkt auf demselben geometrisch-optischen Strahl liegen, aber 


P, vor und P, hinter dem Brennpunkt ist, sich 


(13) Up, = — Up, errr 
. 
ergibt, wobei ö, die Phasenanomalie in P, ee 


Auf Grund des Babinetschen ag läßt sich noch eine 
andere, ähnliche Beziehung ableiten. Es trete durch die im Unend- 
lichen liegende Blende vom halben Öffnungswinkel « eine konver- 
gente Kugelwelle. Die dadurch hervorgerufene Lichtschwingung in 
einem Aufpunkt P sei wieder up Nun falle aus entgegengesetzter 
Richtung ein Lichtkegel vom halben Öffnungswinkel 1 — « ein, der 
durch einen kreisférmigen Schirm aus einer konvergenten Kugel- 
welle ausgeblendet wird. Die dadurch in P erzeugte Lichtschwingung 
sei vp. Dann ist die Summe von u, und vp gleich der Lichtschwingung 
bei allseitig einfallender Kugelwelle. Diese ist unabhängig von # 
und kann daher aus (9) berechnet werden, indem man dorte=n 


. 


Falle und 2%+6> stets |sin(2 +0] >s” ist. Sonst genügte 


überall die Voraussetzung cosa > s~ ‘was für a+e< stets erfüllt ist. 


Ich gebe nun hier einen andern Beweis für |sin (2 # + 9] >s” ‘ls. der nicht 
mehr die Voraussetzung cos«e>s” ‘ls gebraucht, sondern ebenfalls mit 


a +e=— auskommt. Es ist nämlich A: 
2 


und somit 
,@ 


|sin (29 + 8) >sin n—e)=sine= 2cos „sin, >sin 


- 


Es ist daher in den Voraussetzungen (31) auf $.20 der erwähnten Arbeit 


cosa zu ersetzen durch mit s= 2-aresin 8” '“ Diese 


Abänderung bedeutet für die unten folgende Tabelle der u,-Werte eine Ver- 


einfachung der Einteilung der «-Werte. at 
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Sei nun der Aufpunkt P’ auf demselben geometrischen Strahl wie P 
und in gleicher Entfernung vom Brennpunkt, aber auf der entgegen- 
gesetzten Seite, so ist die Lichtschwingung up dann gleich vp, und somit 


an 


(15) Up+ Up = = -sin (kr), 

wenn die Lichtschwingung up in P’ (Poldistanz # = a — #) erzeugt 
wird von einem von vorn') einfallenden konvergenten Lichtkegel 
vom Offnungswinkel « = rn — «, während die Lichtschwingung up in 
P (Poldistanz +) erzeugt wird von einem aus gleicher Richtung ein- 
fallenden Lichtkegel vom Öffnungswinkel «. 


Nachdem die Werte für up für die Fälle 0<u<z und 


9 


vorliegen, läßt sich zunächst mit Hilfe von (13) wp für 
und 2<9 <a angeben und schließlich mit Hilfe von 
(15) auch für - < und OS Es bleibt noch « = 
übrig. Dafür fir *=0 gilt (9. Im Fall = wird 
nach (15) up= up, also 


Up = — (kr). 
r 


Der Fall O0< a= läßt sich analog u 0< er- 


edigen ?): Es sei 
7 
ikr 


Damit u, für große s(= kr) näherungsweise gleich am wird, braucht 


nur gezeigt zu werden, daß die folgenden Größen nicht kleiner als 
sind’): sin é, isin "— Isin(s=" + |sin(2e — 


sin (2+ +, | sin _ a) ‚sin (3 5 + 9) . Fir geniigend groBe s 
ist nun 
1. sine = 2cos “sin => sin+ = s-"; 


2 2 2 
2. sin — > sin (s~'/* — = sin — > sin (0,9-s—") 


1) $= 0 ist die Lichtkegelhalbachse vor, 9 = 1 diejenige hinter dem | 
Brennpunkt. 
2) H. W. Breuninger, a. a. O., S. 26. 


3) H. W. Breuninger, a.a. O. , Bedingungen (33) S. 20, (35) S. 23, (35’) 8.26. 
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bd 1/ 3 1 1 1 


4. Wenn 2% —e=- 5 , dann ist 
|sin (2H — 10015 
weil für < (was angenommen werden darf) 


& < 2,003 m 


dann ist auch 2% — «= — und 


ist; ; 


5. Ist 2%+¢> , dann ist wegen = 

ist 24 — dann ist auch 2% +&> und 

6. SchlieBlich ist u 


sin (3 _ 9) = sin (3 + 9) > sine =>s- 


Nach (13) ergibt sich nun auch up für = + < P<. 


Die folgende Tabelle gibt eine Übersicht über die Werte von up. 
Dazu ist zu bemerken: 

bedeutet, daß die für dieses Fach geltende Wertekombination 
von « und # sinnlos ist, =0, 
_ | bedeutet, daß bei dem zugehörigen «-Wert das #-Intervall in 
mehrere Teilintervalle zerfällt, für welche up je einen besondern 
Wert annimmt, z.B. +eS as z —eund 
hier zerfällt das w#&-Intervall in se’ 
a+eStSa—a—e; der Wert für up ist daher in diesen drei 
Fächern zu suchen. 

Die fett gedruckten Werte gelten fir beliebiges r; die andern 
sind Grenzwerte fiir groBes r. 

Sämtliche verwendeten Abschätzungen') gelten für s = 25-10%, 
d.h. bei 2=5-10-5cm für r=200 cm. Für s = 25-10° wird 


1) Mit einer Ausnahme, nämlich der in Breuninger, a.a.O., S. 29 auf- 
tretenden Abschätzung "> Diese Ungleichung ist aber in der 


erfüllt ist. 


Fußnote 1) S. 234 durch eine ersetzt worden, 109° 
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pin Abhängigkeit von « und # 


n 
0 eikr_ eikrcosa eikr_ 1 eikr _ zikrcosa 2é-sin (kr) 
= 0 < 28-8 
r | r | r r 
| 
<u-e r r = 
=. 0 é-sin (kr) Shand 2é-sin (kr) 
2r r 2r r r 
=ar+e 
en—a—e 
=n-—ate 4 -ikr —ikr | wat 2i-sin (kr) 
=n-ates es 2i-sin (kr) 3i-sin (kr) 
a—8 | | 
0 i-sin (krı | ekr 2i-sin (kr) 
2r r | r 2r r 
<n—8 r r 
| 
0 e-ikrena e-ikr e-ikreosa, 2é-sin (kr) 
=n 
r hr, r r 


welchem der stetige Übergang von (9) zu (10) erfolgt. Dieses 


= 


» = 0,119 = 6%. und & = 6*»°. Will ich nun wp längs eines nicht- 
axialen geometrischen Strahls untersuchen, dessen Poldistanz kleiner 
als 6*«° ist, dann habe ich s so groß zu wählen, daß s—’s» = # wird, 
damit die Näherung (10) einigermaßen erfüllt ist. Je weniger sich # 
von 0 unterscheidet, desto größer muß s sein. Im Endlichen ist up 
stetig; denn zwischen den Winkeln, für welche z. B. die Formel (9) 
und die Näherungsformel (10) gilt, ist ein Intervall auszunehmen, in 


Intervall 


wird desto kleiner, je größer s ist. Im Unendlichen selbst ist up 
allerdings unstetig, was auch zu erwarten ist; denn die ganze Ab- 
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leitung erfolgte auf Grund der Kirchhoffschen Formel mit den 
zugehörigen bekannten Annahmen, wodurch voraussetzungsgemäß im 
Unendlichen eine Unstetigkeit vorbanden ist. (10) gilt, allgemein 
ausgedrückt, für genügend große s in dem durch die beiden Rotations- 


flächen s und F=a—s («= $) begrenzten Raum; nach 


den obigen Beweisen ist o = —; es ist aber ohne weiteres möglich, 


1 
daß für o auch größere Werte als '/, zulässig wären; der Beweis- 
gang müßte entsprechend geändert werden. 

Unter der Voraussetzung kr - sin #< 1 und unter der Annahme 
von kleinem « gilt in der Umgebung der Achse eines konvergenten 
Lichtkegels, der durch Ausblendung einer Kugelwelle durch einen 
im Unendlichen gelegenen Kreisschirm entsteht, folgende Näherungs- 


2.2 gin? 

16) 1 {w, — k?r? gin? yp 

ni 

W — eikrcosd _ gikrcosd cosa ins 


kr - cos? k?r?. cos? 


(sine _ 2tcose 2 )- 
= 


44 


kr-cosd cos? 


Diese Formel gilt auch in der Brennpunktsebene (a= 3) in der 


Umgebung des Brennpunktes. Es läßt sich aus (16) fir #=a<1 
noch eine Näherungsformel für ö ableiten, nämlich 
kro? kr 3 


- vii 
(17) 


~ 


5 


Diese Formel besagt, daB lings einem nichtaxialen Strahl in der 
Umgebung des Brennpunktes in der Entfernung r, vor dem Brenn- 


1) H. W. Breuninger, a. a. O., 8.31; Debye hat eine etwas davon ab- u ah : 
weichende Näherungsformel angegeben, nämlich f 
r- cost 4 
= eikreos® _ ‚ikr cosd cosa 
kr.cos# k?r? cos? & 
kreosdcosa k*r*cos*dcos*a) 
= 


| 
— 
(kr) 
(kr) 
ur & 
(kr) 
= 
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punkt annähernd die Phasenanomalie vorhanden ist, die in einer 
Entfernung r, (>r,) auf der Achse vor dem Brennpunkt herrscht, 
wobei die Differenz r,— r, desto größer ist, je größer der Winkel 
zwischen dem nichtaxialen Strahl und der Achse ist. 

Weicht die Blendenform von einem Kreise ab, so ändert sich 
das Verhalten von u, Nach Rubinowicz') existiert bei einer 
konvergenten Kugelwelle zwischen genügend weit vor und hinter 
dem Brennpunkt gelegenen Aufpunkten stets eine Phasenanomalie x 
bei beliebiger Gestalt der Blende mit Ausnahme einiger spezieller 
Blendenformen, wozu auch die Kreisblende gehört. So tritt im 
Gegensatz zur Kreisblende bei einer rechteckigen Blende längs der 
Achse weit genug hinter dem Brennpunkt die Phasenanomalie auf. 

Die theoretische Ableitung der Phasenanomalie z. B. bei Gouy 
erfolgte durch Einteilung der Wellenfläche in unendlich schmale 
parallele Streifen, wobei keine Abhängigkeit der Länge dieser Streifen 
von ihrer Lage angenommen wurde. Daher gilt die Ableitung bei 
Gouy für eine rechteckige Blende, weswegen er auch die Phasen- 
anomalie z erhielt. 

Möbius hat ö für den Fall untersucht, daß die konvergente 
Welle Wellenflächen verschiedener Gestalt besitze. Für die kon- 
vergente Kugelwelle erhält er zunächst in dem Kapitel „Berechnung 
der Phasenanomalie einer Sattelfläche“ nach einem Rechenverfahren, 
das ähnlich wie bei Gouy im Prinzip eine Rechtecksblende voraus- 
setzt, weit hinter dem Brennpunkt die Phasenanomalie x. In dem 
darauf folgenden Kapitel „Berechnung der Phasenanomalie einiger 
anderer Flächen“ wird für Rotationsflichen ein Verfahren an- 
gewandt, das einer Kreisblende entspricht, ohne daß dies erwähnt 
wird. In der dazu gehörigen Tabelle wird als Phasenanomalie einer 
Kugelwelle hinter dem Brennpunkt der Wert 2 angegeben, obwohl 
sich jetzt wegen der Kreisblende für J ein Wert entsprechend den 
Formeln (1), (6) oder (12) ergeben müßte. Es ist daher angebracht, 
das Rechenverfahren von Möbius durchzunehmen, wobei wir es 
gleichzeitig verallgemeinern. 

Abb. 3 stellt im Schnitt die durch die Blende gehende Wellen- 
fläche dar. Nach Möbius ist die Phasenanomalie ö in unendlicher 
Entfernung von dieser Wellenfläche gegeben durch?) 


do 
(18) = arctg — argi [ e-ir.de, q 


do 


1) A. Rubinowicz, Phys. Rev. 54. S. 931. 1938. 
2) W. Möbius, a.a. O., 8. 109. 
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Es ist dabei über die Blendenöffnung zu integrieren. Ist die Wellen- 
fläche ein Zylinder, dann sei die Blende u Rechteck von der 


(19) do =hy(dt?+ (dd); | 
ist sie eine Rotationsfläche, dann sei die Blende kreisförmig wel 


(20) do = 2n& 


Konkave Wellenflächen unterscheiden sich von konvexen gleicher 
Form nur durch das Vorzeichen von 4; do ist stets positiv; ist 


daher ö, die Phasenanomalie einer alex und Ö, diejenige einer 
konvexen Welle in unendlicher Entfernung von der Blende, dann 
ist nach 18) - 

Dies entspricht (14); denn von einem axialen Aufpunkt vor dem 
Brennpunkt erscheint die Wellenfläche einer konvergenten Kugel- 
welle konvex, von einem Achsenpunkt hinter dem Brennpunkt aber 
konkar. 

Wir beschränken uns wegen (21) auf die Berechnung von Ö,. 
Dabei werden wir davon Gebrauch machen, daß 


<— 


4 
i = 
at 


= 


= 
Mt 


(22) fe e-*. dz, 
Pur 
i] 

integriert über die positive imaginäre Achse, gleich demselben Integral, —_— 4 
erstreckt über die positiv reelle Achse ist, wenn —1 << n<0O ist’. 

1 
1) Vgl. W. Mobius, a. a. O., S. 94. 
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Es ist nämlich das Integral über ein Stück der positiv imaginären 
Achse gleich demjenigen über ein entsprechendes Stück der reellen 
Achse, vermehrt um dasjenige über zwei Viertelkreisbogen (Abb. 4) 
mit den Radien o, und o,. Auf den Bogen ist z = 9 - ei”. ie 


Ist 9 = o,, so betrachte ich lim 9 0, = 00. Der Betrag von (23) ist 
höchstens gleich 


a/2 


0 0 ald 


Die drei letzten Glieder verschwinden für n<0 mit o,-> o. 
Ist n = 0, dann ist 


erstreckt über die positiv imaginäre Ashen, gleich 


Br erstreckt über die reelle 2’-Achse; dies ist verschieden von dem über 
die positiv reelle z-Achse erstreckten Integral (24). Desgleichen 
darf im allgemeinen auch für n > 0 die Integration über die positiv 
imaginäre Achse in (22) nicht ersetzt werden durch diejenige über 
die positiv reelle Achse. 

Betrachten wir zunächst konvexe Zylinderflächen mit der zur 
Mantellinie normalen Schnittkurve 


(25) gq=k4=|k£r. 


Dann ist nach (19) 


Wie Möbius berücksichtigen wir nur die kleinste Potenz von 4 
Dann wird aus (6) | . 


ist das Integral über einen Bogen gleich re 
(23) fe ecosp, e-iesinp, gi(ntly, dq. 

4 
Dies verschwindet aber, wenn 9 = 9, und lim op, = 0 ist fir -—1 <n. — 


242 Annalen der Physik. 5. Folge. Band 35. 1939 ar 4 


A=. 
H 
Dt ie 
4 
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1 
(2 
H 
2 
f 
: 
(: 
I 
P 
if 
i 07) 
— 


en 
für dh O<m<il 
4 
d. h. m= 1 
so 
m>1 
Daher wird nach (18) 
kD 
(27) Ö, fir O<msiI 
n. 
(28) O,=argij/q -dq fir m>1 
? Ich betrachte zunächst m > 1 und setze iq =z; dann wird i 
ikD u 
N. Hier ist vom Nullpunkt aus über ein Stück der positiv imaginären 
z-Achse zu integrieren. Wie oben gezeigt, kann dafür im Grenz- 
fall kD = oo über die positiv reelle Achse integriert werden. Dieser 
Grenzfall ist um so mehr angenähert, je größer m und je größer kD, 
d.h. je größer die Blende ist. In diesem Grenzfall wird somit 
m-1 
Br 6, =arge™ , Ny . 
wobei dessen kleinster positiver Wert zu nehmen ist!), also bat 
2. 
N 1) Der Vektor A 
r > für m > 1 im 3. oder 4. Quadranten; denn sein Imaginärteil ist negativ, . 
4 
0 
positiv ist. Daher muß nach (28) ö, im 1. oder 4. Quadranten liegen. Für a “nl 
m = 2 liegt kein Zweig voni ™ im 4. Quadranten; hier ist die Bedingung, u : 
m—1 
daß der kleinste positive Wert von argi ™ zu nehmen ist, richtig. Für . . 
andere m folgt diese Bedingung daraus aus physikalischen Stetigkeitsgründen. Er 
16* 
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(29) und (30) bedeuten, daß die Phasenanomalie 5 Zylinderwellen, 
die durch (25) gegeben sind, sich dann, wenn der Aufpunkt längs 
der Achse unendlich weit wegläuft, im Fall m > 1 einem bestimmten, 
im wesentlichen nur von der Wellenfläche abhängigen Wert nähert, 

? während im Falle m=1 auch die Blendengröße eine wesentliche 
Rolle spielt. 
En Betrachten wir nun eine Fläche, die durch Rotation der durch 
(25) gegebenen Kurve um die zur $£-Achse senkrechte Gerade ent- 


steht, dann ist — 
k2. — m q Mm +1. dq. 


es Barbi man nur die kleinste Potenz in q, dann wird 


k?.do=2ng" -dq fir 0O<m<l; 43: 


2 m—2 AN 
k do = q dq 
ist, wird < 


> 
Ich untersuche diesen Fall nicht näher. Ist m > 1, dann wird 


kD 
3 0 
Wenn hier weiter m > 2 ist, dann kann dieser Fall wegen — =—* <0 


wie oben der Fall m > 1 der Zylinderwellen behandelt werden, und 
es wird 


n m-?2 
(31) 
Ist m = 2, dann wird 
kD 
v 


wa 


| 
. ar Mobius hat nur die Fille m = 4; 2; 3 und 4 behandelt. un 
m = 1 ist, (33 
Ist 
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und nach (21) 


kD k 
Ist die Wellenfläche eine Kugel vom Radius f, so kann sie dar- 
gestellt werden durch 
2f4- 
Da nur die kleinste Potenz in q, bzw. A berücksichtigt wird, wird 
daraus 2f 4 = $?, woraus folgt 
R? 
D 
Dies in (33) eingesetzt, ergibt 
was auch aus (6) folgt, wenn man p = 00 setzt"). u 


Möbius . außer der Kugel bei den Rotationsflächen noch 
die Fälle m==, 3 und 4 behandelt. Seine Ergebnisse fiir die 


beiden letztern (m > 2) sind richtig, während für m = + dies nicht 
zutrifft, weil dann (31) nicht gilt. Für m = + wird die Phasen- . er 
anomalie in unendlicher Entfernung ebenfalls wesentlich abhängig : 
von der BlendengréBe wie bei der Kugelwelle. 

Denkt man sich die Blende, durch welche die Lichtwelle ein- 
tritt, im Unendlichen und die Kaustik im Endlichen, dann ist in 7 
Punkten weit vor der Kaustik die Phasenanomalie ö = 0,, weil von “" 
diesen Punkten aus die durch die Blende eintretende Wellenfläche 
konvex erscheint, während in Aufpunkten weit hinter der Kaustik 
ö=6, ist. Die Phasenanomalie zwischen den Punkten weit vor 


und denen weit hinter der Kaustik beträgt daher ö, —ö,, was z.B. 


im Falle einer Kreiszylinderwelle nach (29) (m = 2) 7 ergibt. 


> 

Betrachten wir das Verhältnis zwischen (31), gültig für m > 2, 
und (32), gültig für m = 2. Habe ich eine Wellenfläche mit großem m, 
dann berechnet sich nach (31) ö, mit einer gewissen Genauigkeit. 
Deformiere ich die Wellenfläche so, daß m stetig abnimmt, dann 
muß ich, damit sich ö, nach (31) mit konstanter von m unabhängiger 
Genauigkeit ergibt, die Blendenöffnung mit abnehmendem m ver- 
größern, weil ja (31) desto ungenauer wird, je kleiner m ist, und 
desto genauer ö, annähert, je größer die Blende ist. Deformiere ich os 


1) Setzt man in (31), was nur für m > 2 gilt, m = 2, dann wird für die a 5 
Kugelwelle 5, = 0 und 6, = 7; also wäre danach die Phasenanomalie gleich 1 4 
fiir eine konvergente Kugelwelle. 


> 
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so weit, daß m nach 2 geht, dann muß ich die Größe der Blende 
nach Unendlich konvergieren lassen, und für diesen Grenzfall ergibt 
sich nach (31) ö, = 0. (32) hingegen gilt für m = 2 und endliche 
Blendenöffnung. Für eine unendliche Blendengröße ergibt sich nach 
(32) für m = 2 ein unbestimmter Wert für ö,, so daß wenigstens 
kein Widerspruch zwischen den beiden Formeln besteht. 
Vergleichen wir nun noch die theoretischen Berechnungen mit 
den Experimenten. Von Zeeman!) wurde festgestellt, daß die Mitte 
der Interferenzringe der sich überlagernden beiden Kugelwellen mit 
verschiedenen Brennpunkten an jedem Brennpunkt von hell nach 
dunkel oder umgekehrt umschlägt, während sich theoretisch ergab, 
daß bei einer Kreisblende ö nicht nur in der Nähe, sondern auch 
in größerer Entfernung vom Brennpunkt beträchtliche Schwankungen 
ausführt. 

Man?) habe zwei Kugelwellen gleicher Wellenlänge, gleicher 
Ausbreitungsrichtung, mit den Brennweiten f, und f,, mit der 
Phasendifferenz £ in der Blende und den Phasenanomalien 0, und Ö,; 
sie seien entstanden aus einer einzigen linear polarisierten Welle 
durch Passieren einer doppelbrechenden Linse und seien daher 
senkrecht zueinander linear polarisiert. Dann ist die Interferenz- 
erscheinung längs der Achse hinter einem Analysator, dessen 
DurchlaBrichtung den Winkel zwischen den beiden Polarisations- 
richtungen halbiert und auf der Durchlaßrichtung des vor der Linse 
befindlichen Polarisators senkrecht steht, gegeben durch 


2A,’ fı? Qa. 2 fe? 2 Ver 
I = sin? 
wobei I die A, und A, sind 
und entsprechend (4) u 


ist. Denkt man sich an Stelle der WAKE, Linse eine 
gewöhnliche, dann ist A, = A,, fi =f,, &= 0 und nach (34) J = 0, 
da durch die beiden aufeinander senkrechten Nikols das Licht voll- 
ständig ausgelöscht wird. Ist Ap die Entfernung zweier Aufpunkte, 
an denen aufeinander folgende Minima einer Kugelwelle liegen, 
so ist 


Zeeman, Phys. Ztschr. 1. 8. 542. 1900. 
2) Vgl. H. W. Breuninger, a. a. O., S. 36—42. 
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Für 4=5-10-5 cm, p= 33 cm und R=1,4 cm (Zeemanscher 
Versuch) wird 4p = 0,06 cm. Sind nun im Experiment die Brenn- 
punkte nicht ganz scharf, dann sind die Extremalstellen der Inten- 
sität etwas verwaschen, und man wird bei diesem kleinen 4p nur 
einen Mittelwert der Intensität wahrnehmen. Dieser ergibt sich zu 


4 fi? A,’ f,” 
(35) 
| ‚cos (#5 


wobei diese Formel nicht in der nn der Brennpunkte gilt. 
Um hierin im Falle des Zeemanschen Versuches i 
4 


(36) 


berechnen zu können, müßte bei einem bestimmten £ der Radius R 
genauer als auf 0,1°/, bekannt sein. Um daher einen Überblick zu 
bekommen, was theoretisch beim 20 


Versuch von Zeeman zu erwarten I bi 
ist, sind unter der Annahme 
A, = A, = 1 fünf Kurven a, b,c, d, e 150 di 


(Abb. 5) berechnet worden, welche 
für einen Wert des Produktes (36) 
von bzw. 


700 | 


-1,-4,8, +7, +1 | 
gelten. Kurve c gibt die Inten- 

sität, die bei reiner Addition der N 

Intensitäten der beiden Kugel- yy | 

wellen entsteht. a und b verlaufen 0 G AV 5 
vor dem ersten und hinter dem W ht HH 

zweiten Brennpunkt oberhalb vonc, Abb. 5 

zwischen beiden Brennpunkten 

unterhalb von c. Dies besagt, daß bei a und b vor dem ersten 
und hinter dem zweiten Brennpunkt relativ Helligkeit und zwischen 
beiden Brennpunkten Dunkelheit herrscht; das Umgekehrte gilt 
für d und e. Die Ursache liegt darin, daß an den Brennpunkten 
das dritte Glied der rechten Seite von (35) das Vorzeichen umkehrt. 
Der Grund dafür ist folgender: Beide Wellen haben an jedem 
Minimum ihrer Amplituden einen Phasensprung a. Diese Minima 
sind nach (5) bestimmt durch Q = 2n mit ganzem positivem oder 
negativem n. Im Brennpunkt ist n = 0; daher entspricht dieser 
einer Minimalstelle; hier findet aber kein Phasensprung statt, weil 
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die Amplitude von Null verschieden ist. Daher ist die duch. 
schnittliche Phasenanomalie einer Welle hinter ihrem Brennpunkt 
um a größer als davor; da nun im Falle des Zeemanschen Ver- 
suches nur eine durchschnittliche Phasenanomalie beobachtet werden 
kann, hat es den Anschein, als ob am Brennpunkt eine Phasen- 
umkehr stattfinde. 

Anders als bei Zeeman liegen die Versuchsbedingungen bei 
Sagnac. Hier ist f, = 200 cm und f, = 300 cm. Nach (5) liegen daher 


E (+ )| 2 Minima einer jeden Welle zwischen den Brenn 


punkten, wenn R = 0,4 cm 


| und A=5.10°cm ist. Abb.6 
0 gibt unter der Voraussetzung 
| A, =A, = vier nach (34) 

00 berechnete Kurven a, b, c,d 
\ fir diesen Sagnacschen 

40 Wel m Fall, wobei £ bzw. = 120°, 
\ 60°, 30°, — 60° gesetzt ist. 

2 Überall sind 2 Minima, aus- 
genommen a, wo beide zu- 

> >  sammenfallen. Dies stimmt 

200 200cm im großen überein mit der 
/2 Angabe von Sagnac}), daß 


er zwischen den beiden 
Brennpunkten 4 mal einen 
beinahe plötzlichen Wechsel der Intensität von hell nach dunkel 
oder umgekehrt beobachtet habe. 

Zum Schluß fassen wir die Eigenschaften der Phasenanomalie 0, 


sammen ?): 


) Kreisblende 


von der Blende ein Winkelintervall, das selbst vom Abstand von 
der Blende (bzw. vom Abstand vom Brennpunkt) abhängig ist und 
dessen Mitte vor dem Brennpunkt sich um a von der Mitte hinter 
dem Brennpunkt unterscheidet. 

Die Breite dieses Intervalles beträgt in Brennpunktsnähe 2. 
Mit zunehmender Entfernung vom Brennpunkt schrumpft es zusammen. 


7 1) G. Sagnac, Compt. rend. 138. S. 620. 1904. 


i =e 2) Vgl. H. W. Breuninger, a. a. O. S. 42—45. 


die beim Ausblenden einer konvergenten Kugelwelle 


Langs der Achse Fe ö mit veränderlichem Abstand 
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Diese Zusammenschrumpfung geht vor dem Brennpunkt bei 


kt Annäherung an die Blende nur so weit, daß die Intervallbreite in 

T- der Nähe der Blende noch mindestens 2/3 beträgt. Hinter dem 

en Brennpunkt schrumpft das Intervall auch etwas zusammen. Während 

n- aber vor dem Brennpunkt ö praktisch unendlich oft pendelt, ge- 
schieht das bei im Endlichen gelegener Blende nur endlich oft 

ei hinter dem Brennpunkt. 

er Bei Verlegung der Blende ins Unendliche schrumpft das von 6 


durchpendelte Intervall längs der Achse nicht zusammen (vgl. Abb. 2). 
Gleichgültig, ob daher die Blende im Endlichen oder Unend- 
r lichen liege, gibt es keine so große Entfernung E, daß zwischen den 
Achsenpunkten einerseits, die in einer Entfernung größer als E vor 
8 dem Brennpunkt liegen, und den axialen Punkten andererseits, die 
4) in einem E übertreffenden Abstand hinter dem Brennpunkt sind, 
stets annähernd die Phasenanomalie gleich a sei. Daher lassen 
sich die Zeemanschen Experimente nicht mit der obigen Phasen- 
anomalie allein erklären, sondern die scheinbare Phasenanomalie z 
am Brennpunkt kommt nur dadurch zustande, daß 2 Kugelwellen 
mit verschiedenen Brennpunkten sich überlagern. Die durch ihren 
Brennpunkt gehende Kugelwelle läßt hier einen sonst periodisch 
aufgetretenen Phasensprung a einmal ausfallen, so daß sich die 
durchschnittliche Phasenanomalie hinter dem Brennpunkt von dem- 
selben Wert vor dem Brennpunkt um unterscheidet. > 
ö längs nichtaxialen Strahlen geht stetig aus dem Verhalten 
längs der Achse hervor. Denn für einen nichtaxialen Strahl gibt 
es zwar eine Strecke E derart, daß zwischen den Punkten vor dem 
Brennpunkt in einer Entfernung größer als E einerseits und den 
| Punkten andererseits, die hinter dem Brennpunkt in einer Eüber- 
treffenden Entfernung liegen, annähernd die-Phasenanomalie 2 be- " 


steht. E muß aber desto größer sein, je kleiner der Winkel 
zwischen dem Strahl und der Achse ist, und wird unendlich, wenn 
diese beiden zusammenfallen. Das von Ö durchpendelte Intervall 
| schrumpft lings nichtaxialen Strahlen auch zusammen, wenn die 
Blende im Unendlichen liegt, was in diesem Fall lings der Achse 
nicht mehr eintritt. Man darf daher annehmen, daß, ob die Blende 
im Endlichen oder Unendlichen liegt, das von ö durchpendelte 
Intervall längs einem Strahl desto rascher mit der Entfernung vom 
Brennpunkt zusammenschrumpft, je größer der Winkel zwischen Achse 
und Strahl ist. Dies stimmt damit überein, daß in Brennpunktsnähe 
in einer Entfernung r, vom Brennpunkt nach (17) nahezu dasselbe 6 
gilt wie in einer Entfernung r, auf der Achse, wobei r, — 1, desto 
größer ist, je größer der Winkel zwischen Strahl und Achse ist, 
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was vermutlich auch noch in größeren Entfernungen vom Brenn- 
punkt richtig ist. 

Bei im Endlichen gelegener Blende wird bei endlicher Aus- 
dehnung der Lichtquelle die Lichtschwingung in einem axialen Punkte P 
nicht genau gleich der unter Voraussetzung punktförmiger Lichtquelle 
berechneten Größe, sondern gleich dem Durchschnittswert der für 
die Punkte eines in P senkrecht zur Achse stehenden Flächen- 
elementes berechneten Lichtschwingungen. Da nach oben aus- 
gesprochener Vermutung die Phase in achsenbenachbarten Punkten 
gleich der Phase in etwas weiter vom Brennpunkt weg liegenden 
Achsenpunkten ist, kann die Mittelwertbildung über ein Flächen- 
element senkrecht zur Achse ersetzt werden durch die Mittelwert- 
bildung über ein gewisses Längenelement der Achse. Je größer 
das Bild der Lichtquelle in der Brennebene, gemessen durch den 
Winkel, unter dem es vom Blendenmittelpunkt aus erscheint, im 
Verhältnis zur Blendenöffnung ist, desto länger ist das Achsenstück, 
über das gemittelt werden muß, um die Lichtschwingung in P zu 
erhalten. Man sieht so, wie durch Vergrößerung der Lichtquelle 
die axiale Interferenzerscheinung mehr und mehr verschwommen wird. 


b) Rechtecksblende 
Hier konvergiert d auch längs der Achse in großen Ent- 
fernungen vor und hinter dem Brennpunkt nach 0 bzw. m. 
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Messungen von elektrischen Anregungsfunktionen 
an Quecksilber') 


Von Pascal Vetterlein 


(Aus dem Physikalischen Institut der Technischen Hochschule Hannover) 
(Mit 15 Abbildungen) 


Inhalt: $1. Problemstellung. — $2. Grundsätzliches zur Methode. — 
83. Beschreibung der Apparatur. — $4. Auswertung und Fehlerquellen. — 
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81. Problemstellung 


Bei der Untersuchung der Ausbeute unelastischer StéBe zwischen 
Elektronen und Atomen kann man von zwei wesentlich verschiedenen 
Gesichtspunkten ausgehen: Man kann einmal die Geschwindigkeits- _ 
verteilung der Elektronen nach Durchlaufen des Stoßraumes unter- | 
suchen und findet die Ausbeute für die Anregung bestimmter Alom- 
terme, man kann zweitens die Intensität der infolge der Stöße 
emittierten Spektrallinien messen und findet die Ausbeute für de — 
Emission bestimmter Lichtquanten. Die quantitative Beschreibung 2 


it 


der Ausbeute erfolgt dadurch, daß dem betreffenden Ereignis, also 
der Anregung eines Terms — bzw. der Emission eines bestimmten 
Lichtquants — in bekannter Weise ein Wirkungsquerschnitt zugeordnet 
wird. Diese Wirkungsquerschnitte für die Anregung eines Terms 
bzw. für die Anregung einer Linie hängen sehr wesentlich von der 
Energie der stoßenden Elektronen ab und werden — als Funktionen 
dieser Energie betrachtet — als Anregungsfunktionen bezeichnet. 
Man spricht, je nachdem es sich um die Anregung eines Terms 
oder um die Emission einer Linie handelt, von elektrischen bzw. 
optischen Anregungsfunktionen. 

Fraglos geben die elektrischen Anregungsfunktionen unmitte- 
barere Aussagen über die Eigenschaften des Atoms als die optischen 
Anregungsfunktionen. Denn in die letzteren geht im allgemeinen 
infolge der Kaskadensprünge außer der elektrischen Ausbeute für 
alle Terme auch die Verteilung der Übergangswahrscheinlichkeiten 
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in komplizierter Weise ein, so daß nicht ohne weiteres von der 
optischen Anregungsfunktion auf die elektrische Anregungsfunktion 
geschlossen werden kann. Trotzdem sind optische Anregungsfunktionen 
in einem ungleich größeren Umfang gemessen als elektrische An- 
regungsfunktionen, bei denen sich das Material auf ganz wenige 
Messungen beschränkt. Das ist fraglos darauf zurückzuführen, daß 
die Messung der elektrischen Anregungsfunktionen nach den bisher 
benutzten Methoden auf weit größere Schwierigkeiten stößt, als die 
Messung optischer Anregungsfunktionen. Eine Verbesserung der 
Methoden zur Messung elektrischer Anregungsfunktionen schien 
deshalb von wesentlicher Bedeutung. 

Zur Bestimmung des Wirkungsquerschnitts muß die Geschwindig- 
keitsverteilung gemessen werden, mit der die Elektronen den StoB- 
raum verlassen. Für die Messung der Geschwindigkeitsverteilung 
stehen grundsätzlich zwei Möglichkeiten zur Verfügung: Man kann 
erstens den aus dem Stoßraum austretenden Elektronenstrahl durch 
elektrische oder magnetische Felder räumlich in ein Geschwindigkeits- 
spektrum auseinanderziehen, zweitens kann man die Geschwindigkeits- 
verteilung dadurch bestimmen, daß man die austretenden Elektronen 
gegen ein variables Gegenfeld anlaufen läßt. Die Messung im Gegen- 
feld bietet vor der spektralen Zerlegung in verschiedener Hinsicht 
wesentliche Vorteile (größere Einfachheit und Durchsichtigkeit der 
Anordnung, bessere Ausnutzung der aus dem Stoßraum austretenden 
Elektronen). Aus anderen Gründen sind jedoch der Verwendung 
der Gegenspannungsmethode bei der Bestimmung von Wirkungs- 
querschnitten enge Grenzen gesetzt. Denn die Berechnung der 
Wirkungsquerschnitte aus der gemessenen Geschwindigkeitsverteilung 
ist um so sicherer, je größer die freie Weglänge der Elektronen im 
Vergleich zu der im Stoßraum durchlaufenen Strecke ist, je geringer 
also der Anteil der Elektronen wird, die Energieverluste erlitten 
haben. Nun wird aber die Ausbeute bei der Gegenfeldmessung als 
Differenz zweier Ströme gemessen, und dadurch wird die erreichbare 
Meßgenauigkeit um so kleiner, je besser im Stoßraum die Bedingungen 
für die Sicherheit der rechnerischen Auswertung der Messung realisiert 
werden. In der vorliegenden Arbeit wird eine Modifikation der 
Gegenspannungsmethode entwickelt, die diese Schwierigkeit ausschaltet 
und die Empfindlichkeit der Messung wesentlich erhöht, ohne daß 
die oben erwähnten Vorzüge der Gegenspannungsmethode dadurch 
berührt werden. 

Die neue Methode wird an Quecksilber erprobt, weil sich dort 
neben neuen Resultaten Vergleichsmöglichkeiten mit früheren 
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oy §2. Grundsätzliches zur Methode 

aes Die Geschwindigkeitsverteilung der Elektronen, die pro Zeit- 

a einheit den StoBraum verlassen, wird beschrieben, indem der Teil 

J der pro Zeiteinheit aus dem StoBraum austretenden Elektronen, der 

e auf den Geschwindigkeitsbereich 4 V entfällt, gleich N (V) AV gesetzt 

B wird. Befindet sich in dem StoBraum ein Gas, so wird die Funktion N(V) 

a einen Verlauf zeigen, wie ihn in schematischer Form Abb.1 wiedergibt. — u 

” Dabei ist V, die Geschwindigkeit, My u 

ad mit der die Elektronen in den ey u 
Stoßraum eintreten. Es ist ferner | 
der Einfachheit halber angenom- 


men, daß nur eine Art von un- | 
elastischen Zusammenstößen (mit ho 


Energieverlust V,— V, = U,) statt- N: 
finden kann. h, J 

Die Flächeninhalte F, und F, 
unter den beiden Maxima sind 
gleich der Zahl der Elektronen, Abb. 1 7 


die pro Zeiteinheit ohne Energie- 

verlust bzw. mit Energieverlust aus dem Stoßraum austreten. Die 
Form ist bei beiden Maxima nur durch die Geschwindigkeitsverteilung 
bedingt, mit der die Elektronen in den Stoßraum eintreten, und ist 
deshalb bei beiden Maxima die gleiche. Deshalb gilt für das Ver- 


hältnis Er , das zusammen mit der Gasdichte und Stoßraumlänge 
o 1 


den Wirkungs-Querschnitt bestimmt, 
. 


F, hy 


FR+F hth’ 


wo h, und h, die Höhen der beiden Maxima sind. 2 
Bei den üblichen Gegenspannungsmethoden bestimmen wir nun 

aber nicht für jedes V den entsprechenden N(V)-Wert, sondern — 
wir erhalten für jede Gegenspannung V den Summenwert aller — 
Elektronen, deren Geschwindigkeit ausreicht, um gegen das Gegenfeld _ 
a — also aller Elektronen, deren Geschwindigkeit = V ist. 
Far den an der Stelle V gemessenen Stromwert J kann man also 


ta 


1 _ Die Gegenspannungskurve ist somit die Integralkurve der gesuchten 
% Geschwindigkeitsverteilungskurve. Die der oben gezeichneten Ver- 
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teilungskurve entsprechende Gegenspannungskurve würde durch Abb. 2 
wiedergegeben werden. I, und I, — I, sind identisch mit den Inte- 
gralen über die entsprechenden Maxima der Kurve N(V), also mit 

F, und F,. Die Gegenspan- 
| nungsmethode liefert uns also 
auch den zur Berechnung des 
Wirkungsquerschnitts notwendi- 


] gen Quotienten 


vi 


I 
? liefert aber die Ausbeute als 
| | Differenz zweier Stréme. Da- 
| »>/ durch wird die Genauigkeit der 
Messung bei kleinen Ausbeuten 
Abb. 2 außerordentlich begrenzt. 
ee Ich habe deshalb auf Vor- 


schlag von H. Bartels unter Ver- 
wendung eines Gedankens, der in anderem Zusammenhang von 
Dymond!) angewendet war, die Gegenspannungsmethode so modi- 
fiziert, daB N(V) unmittelbar gemessen wird. Kleinheit der Ausbeute 
kann jetzt im Gegensatz zur alten Gegenspannungsmethode grund- 
sitzlich durch Steigerung der Empfindlichkeit ausgeglichen werden. 


Auffanger 
Kathode | 


| 
5 


Abb. 3. Prinzip der Schaltung 


Da überdies N(V) durch eine Wechselspannungsamplitude gemessen 
wird, so ist es weiterhin möglich, die moderne Verstärkungstechnik 
zur Steigerung der Empfindlichkeit einzusetzen. Die Grenzen der 
Genauigkeit sind jetzt lediglich durch Störeffekte gegeben. 

Der neuen Methode liegt folgender Gedankengang zugrunde: 
Der Gegenspannung V in der üblichen Gegenspannungsmethode 
wird eine kleine Wechselspannung AV überlagert (vgl. Abb. 3). 
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1 E. @. Dymond, Proc. of the Royal Soc. Ld. A 167. S. 291. 1925. a 
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Ohne diese Wechselspannung erhalten wir für jeden Wert von V 
den der Gegenspannungscharakteristik entsprechenden Gleichstrom 1. 
Nach Einschalten der Wechselspannung AV schwankt der Auffänger- 
strom I um den Betrag AI, 
der der Differenz AV auf der 
Gegenspannungscharakteristik 
entspricht (vgl. Abb. 4. Macht 
man AV klein genug, so ist die 
Amplitude des Wechselstroms 
AI, der sich dem Gleichstrom I 
überlagert, bis auf einen von V 
unabhängigen Faktor gleich der 
Steigung der Gegenspannungs- 
charakteristik. Variiert man 

V 


also V bei festgehaltenem AV, 

so durchfließt die Primärseite des en 
Transformators T, (vgl. Abb. 3) Abb. 4 
ein Wechselstrom, dessen Am- 

plitude bis auf einen von V inchiehetane Faktor gleich der Ge- 
schwindigkeitsverteilungsfunktion N (V) ist, wenn man von den Fehlern 
absieht, die dadurch entstehen, daß AV endlich bleiben muß, und 
der Grenzübergang 4V —> 0 nicht vollzogen wird. 

Dieser Wechselstrom, durch den N(V) gemessen wird, kann 
so verstärkt werden, daß er mit relativ groben Instrumenten schneller 
Einstelldauer gemessen werden kann, so daß unter Zuhilfenahme 
einer später beschriebenen Potentiometeranordnung eine MeBreihe 
in sehr kurzer Zeit aufzunehmen ist. 


m 


§3. Beschreibung der Apparatur 


Das Vakuumgefäß, in dem die Stoßanordnung untergebracht 
war, bestand, wie Abb. 5 zeigt, aus vier Schliffen, die kreuzförmig 
miteinander verblasen waren. In einen Schliff wurde die eigentliche 
StoBanordnung eingeführt, der gegenüberliegende diente zur Her- 
stellung elektrischer Anschlüsse. Der dritte Schlif verband die 
Anordnungen mit der Pumpe, an den vierten wurde ein kleines 
Vorratsgefäß mit Quecksilber, das mit Wasserkühlung und Thermo- 
meter versehen war, angeschlossen. 

Die eigentliche Stoßapparatur wird in Abb. 6 gezeigt. Die 
einzelnen Teile bestanden aus Messing. Das Messingrohr M, wurde 
in der gezeichneten Weise durch eine wassergekühlte Kittstelle mit 
einem Schliffkern verbunden und in den einen Schliff des oben be- 
schriebenen Vakuumgefäßes eingeführt. M, bildete nun den Träger 
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Abb. 6. StoBanordnung 


unten) und den zur Biindelung der Elektronen dienenden Wehnelt- 
Zylinder M,. Über M, wurde ein Glaszylinder G, geschoben, der 
genau zentrierend das Rohr M, hielt, das den Stoßraum bildete. 
M, wurde zur Kathode hin durch eine Platte mit zentraler Bohrung, 
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der Eintrittsöffnung für die Elektronen, nach der anderen Seite hin 
durch ein Netz aus Al-Draht von 0,06 mm Durchmesser und 1,2 mm 
Maschenweite abgeschlossen. Diesem Netz folgte ein zweites von 
gleichen Abmessungen und schlieBlich der ebenfalls durch das Glas- 
rohr G, zentrierte Auffänger. Dieser bestand aus einem Messing- 
rohr M,, das in gezeichneter Weise durch einen Boden abgeschlossen 
war. Das zweite Netz wurde durch das Glasrohr G, vom Auffänger 
zentriert und gehalten. Die inneren Teile der Messingrohre waren 
zur Vermeidung von Elektronenreflektionen berußt. 

Die Konstruktion der Kathode ist aus Abb. 7 zu ersehen. Die 
Kathode bestand aus einem Platindraht von 0,2 mm Durchmesser. 
Die Emissionsstelle wurde dadurch gebildet, daß der Draht in der 
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Abb. 7. Konstruktion der Kathode 


Mitte verdrillt, die Verdrillungsstelle im Gebläse verschweißt und 
die äußerste Spitze mit einer Spur Bariumoxyd überzogen wurde. 
(Bariumoxyd wurde mit Wasser zu einem dünnen Brei angerieben, 
auf den Platindraht aufgebracht und neben der Flamme eines Bunsen- 
brenners getrocknet.) Der Draht wurde dann in das Steatitrohr S, 
eingeführt, die Oxydstelle durch die Bohrung der Steatitkappe S, 
geführt und das Ganze in der gezeichneten Weise in den Wehnelt- 
zylinder M,, geschoben. Zur Stromzuführung dienten zwei Stecker A, 
und A, an den Enden des Platindrahtes und zwei Buchsen B, und 
B, in der Steatitplatte S,. 

Bei der Konstruktion von Kathode-Stoßraum-Auffänger wurde 
besonderer Wert darauf gelegt, alle Teile leicht auswechseln zu 
können. Die Kathodenkonstruktion hatte daneben noch den Vorteil, 
durch die Form der Steatitröhre im Wehneltzylinder genau zentriert 
zu sein, während andererseits der Wehneltzylinder genau in der Achse 
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des Stoßraums und Auffängers lag. Die Emissionsstelle der Kathode p 

wird nun durch Wärmeleitung geheizt und ist dadurch äquipotential, 2 

durch die Form der Verdrillungsstelle ist sie fast punktförmig. 

Die elektrische Schaltung zeigt Abb. 8. . 

Die Batterie B, lieferte einerseits die Beschleunigungsspannung d 

für die in den Stoßraum eintretenden Elektronen, andererseits zu- te 

sammen mit der vom Potentiometer P, abgoneemnence variierbaren b 
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ir 
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A 

Spannung die Gegenspannung zwischen Stoßraum und Auffänger. 2 
Da die Gegenspannung nur um einen relativ geringen Bruchteil 

(etwa bis 14 Volt) gegenüber der Beschleunigungsspannung variiert A 

wurde, konnten durch die angegebene Schaltung Fehler sehr ge- u 

mindert werden, die dadurch entstehen können, daß sich die Gegen- ii 

spannung und die Beschleunigungsspannung während der Messung ki 
gegeneinander verschieben. Der Wehneltzylinder erhielt durch die 

Batterie B, eine negative Spannung gegenüber der Kathode, die G 

durch das Potentiometer P, geregelt wurde. Durch die Batterie B, ue 

wurde dem Netz zwischen Stoßraum und Auffänger eine positive Vv 

Spannung von etwa 10 Volt gegenüber dem Stoßraum gegeben, um T 
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positive Ionen, die etwa im StoBraum hätten entstehen können, nicht 
zum Auffänger gelangen zu lassen. 

Die Schaltung des Potentiometers P, ist aus Abb. 9 zu er- 
sehen: In den Potentiometerteilen P, und P, sind die Einzelwider- 
stände Präzisionswiderstände von je 1000 Ohm, im Potentiometer- 
teil P, ebensolche Widerstände von je 100 Ohm. Es ist aus der Ab- 


bildung ohne weiteres ersichtlich, daß man die zwischen A und B 


ar Abb. 9. Schaltung Potentiometer 


liegende Spannung schrittweise von Hundertstel zu Hundertstel der 
Gesamtspannung durchlaufen kann, und zwar mit viel größerer 
Genauigkeit, als wenn man die Spannung jeweils mit einem Zeiger- 
instrument einstellen müßte. Außerdem erhöht die Regelung der 
Spannung mit dieser Sonderkonstruktion die Meßgeschwindigkeit 
erheblich, da das Ablesen der Spannung bei jeder neuen Potentio- 
metereinstellung entfällt. 

Zur Erzeugung der überlagerten Wechselspannung AV diente 
ein Röhrengenerator, der durch einen Transformator mit der 
Apparatur gekoppelt wurde. Als zweckmäßigste Frequenz erwiesen 
sich 500 Hz. Die Amplitude der Spannung an der Apparatur be- 
trug einige Zehntelvolt, sie ließ sich durch einen Spannungsteiler 
regeln. 

Der Kondensator C, kompensierte die Kapazität C, zwischen 
Auffänger und Stoßraum, die in der Weise störend wirkte, daß 
auch dann eine Wechselspannung am Transformator T, lag, wenn 
überhaupt kein Elektronenstrom floB, also die Steigung der Stufen- 
kurve gleich Null war. 

Über den Kondensator C, wurde eine kleine Antenne an den 
Gitterkreis der Vorröhre angekoppelt. Durch diese Schaltung 
wurden die Störungsschwingungen kompensiert, die dadurch in den 
Verstärker gelangen, daß die eine Seite der Primärwicklung des 
Transformators T, über die Batterie B, geerdet, die andere Seite 
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jedoch mit den Metallmassen und Leitungen des Potentiometers P,, 
des Röhrengenerators und Auffängers verbunden war. Die Stör- 
_ schwingungen, die durch das Lichtnetz usw. induziert werden, ge- 
_ langen auch über die Antenne an das Gitter der Vorröhre, des 
 Transformators wegen jedoch um 180° phasenverschoben, wodurch 
die Kompensation bei passender Einstellung von C, erreicht wurde. 
Die von außen kommenden Störungen konnten zum größten Teil 
durch diese Maßnahme beseitigt werden. Doch ließ es sich nicht 
ganz vermeiden, daß im Ausgangsinstrument ein kleiner Reststrom 
verblieb, der teilweise aus dem Netzbrummen des Verstärkers, teil- 
weise aus durch das Lichtnetz induzierten Störschwingungen (etwa 
1000 Hz) herrührte. Dieser Reststrom war jedoch völlig konstant, 
m: Lautsprecher leicht zu kontrollieren, er betrug den Maßstab der 
Kurvenbilder gemessen etwa 0,3 Skt. 

Der Wechselstrom im Meßkreis Wy wurde in einem 3stufigen 
Verstärker verstärkt. Mit dem batteriebetriebenen Vorverstärker 
war durch Transformator ein Siemenskraftverstärker gekoppelt, der 
als Endstufe zwei Röhren RE 604 in Gegentakt hielt. Am Aus- 
gangstransformator lag unter Zwischenschaltung einer Kapazität 
von 0,5 uF, die das Netzbrummen herabdrücken sollte, ein Hart- 
mann & Braun-Multavi-Instrument. Außerdem war ein Lautsprecher 
angeschaltet, der eine gute akustische Kontrolle gestattete. 

Die Verstärkungsziffer der Verstärkeranlage war in weiten 
Grenzen konstant und zeigte erst bei sehr großen Werten der durch 
den Primärkreis von T, fließenden Wechselströme ein leichtes Ab- 
sinken. Um die Meßergebnisse sicher auswerten zu können, mußte 
zunächst die Eichkurve aufgenommen werden. Es wurden also durch 
die Primärseite des Transformators T, Wechselströme bekannter 
Amplitude und gleicher Frequenz wie bei den eigentlichen Messungen 
geschickt und zu jeder Stromamplitude der zugehörige Ausschlag 
des Ausgangsinstrumentes abgelesen. Bei der Auswertung der MeB- 
reihen wurde dann so vorgegangen, daß zu jedem Meßwert, d.h. 
Ausschlag des Ausgangsinstrumentes, aus dieser Eichkurve rück- 
warts die betr. Wechselstromamplitude im Meßkreis bestimmt 
und als Wy in willkürlichem Maßstab über der Gegenspannung 
aufgetragen wurde. 

Der herrschende Dampfdruck wurde der Temperatur des Hg 
im Vorratsgefäß entnommen. Die Temperatur wurde mit einem 
Normalthermometer bestimmt und durch Wasserkühlung auf + 0,1° C 
konstant gehalten. Der den abgelesenen Temperaturen entsprechende 
Dampfdruck wurde den Tabellen von Landolt- Börnstein ent- 
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§ 4. Auswertung und Fehlerquellen 
Als typisches Beispiel für die im Quecksilberdampf auf- 
ge- genommenen Geschwindigkeitsverteilungskurven zeigt Abb. 10 zwei 
des Kurven, die bei 60 bzw. 380 Volt gemessen wurden. Als Abszisse 
ırch ist die Gegenspannung V in Volt, als Ordinate der Wechselstrom Wy 
rde. in einem willkürlichen Maß- . 
Teil stab aufgetragen. Das Null- | Wy 
icht maximum bei V = V, ist der 5 6 
rom Ubersichtlichkeit halber in | mw 
-eil- allen Kurvenbildern gegen- | 
twa über dem Verlustmaxima 3s 
ant, um ein Zwanzigstel verklei- 
der nert. Das Nullmaximum ist _ 4 
Nullpunkt für die Bestim- | = . 
. . . 
gen mung der in Voltgeschwin- | er 
ker digkeit gemessenen Energie- y 
der verluste U. Diese Energie- # 80 
us- verluste U sind in einer 
ität neuen Skala mit dem Null- P | oa 
art- punkt in V =V, ein- 
her getragen. | 
Die Kurven in Abb. 10 2 Oo 
ten zeigen, daB sich die ein- 
rch zelnen Verlustmaxima zum i a 
Ab- Teil schon merklich über- 
Bte lagern. Zur genauen Er- ot 
ter querschnitts, der den ein- m 32 3m 3% 378 380 Volt 
8&8 6 4 2 0 
gen zelnen Verlusten zuzuord- * Mebwerte 380V 
lag nen ist, mußte deshalb der 
eß- Kurvenzug mit Hilfe der Abb. 10. Beispiel zweier gemessener > 
h Anfangsgeschwindigkeits- Geschwindigkeitsverteilungskurven 
ck- verteilung, die ja die Form = 
mt aller Maxima gleichmäßig bestimmt, in die einzelnen Verlustmaxima . 
Ing zerlegt werden. Das gelang ohne Mühe, da die Lage der einzelnen + 
Maxima trotz der Störungen in allen Fällen gut erkennbar war, | 
Hg Die Anfangsgeschwindigkeitsverteilung wurde dem Nullmaximum 64 
em entnommen. 
0¢ Das Nullmaximum zeigt bei allen Messungen eine asymmetrische oe 
de Verbreiterung nach wachsendem U. Wird der Hg-Dampf ausgefroren, * 
nt- verschwindet diese Asymmetrie völlig. Abb.11 zeigt das Nullmaximum a 
im Hg-Dampf und nach Ausfrieren. Beide Kurven wurden unmittel- ts @ 


7 
G 


262 Annalen der Physik. 5. Folge. Band 35. 1939 Ee q 
. bar nacheinander aufgenommen. Die Asymmetrie entsteht also 


dadurch, daß sich der völlig symmetrischen, nur durch die Anfangs- 
geschwindigkeitsverteilung bestimmten Form des Nullmaximums ein 
Elektronenanteil überlagert, der kleine Verluste erlitten hat. Die 
Kurven zeigen, daß der Einfluß dieser Verlustelektronen schon 
wesentlich vor der Mitte 
des Nullmaximums ver- 
schwindet. Deshalb läßt 
sich die Verteilung der 
Anfangsgeschwindigkeiten 
| dadurch gewinnen, daß 
man den Abfall des Null- 
maximums nach größerem 
V an der Abszisse des 
Höchstwertes spiegelt. 
| Diese Geschwindigkeits- 
; verteilung zeigt eine Halb- 
wertsbreite von 0,8 Volt, 
dieselbe Halbwertsbreite 
hat auch das Maximum 
| 6,7 Volt, das einzige, das 
| nicht durch benachbarte 
Maxima gestört wird. 
Die Übereinstimmung der 
Halbwertsbreiten darf als 
Beweis für die Richtig- 
keit dieser Zerlegungs- 


— Damf 


Vakuum 


-f methode angesehen wer- 
den. Abb. 12 zeigt eine 

\ Abb. 11 in dieser Weise zerlegte 
Kurve. 


ney Ist Q; der Wirkungsquerschnitt für die Anregung des Terms, 
der zum i-ten Verlustmaximum gehört, so ist 


Dabei ist 1 die Linge ne StoBraumes und p i Druck des 
Hg-Dampfes in mm Hg. Das Integral im Zähler ist, wie der 
Index i andeutet, nur über das i-te Verlustmaximum zu er- 
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strecken, das Integral im Nenner aber über die ganze Verteilungs- 
funktion N (V)). 

Da N(V) nur bis 14 Volt gemessen wurde, ist das Integral 
Vo 
fumav nicht ohne weiteres den gemessenen Kurven zu ent- 
sungen zeigten aber, daB der | 
Bereich zwischen 14 Volt und An 
V, in jedem Fall nur wenig | 
zum Gesamtintegral beiträgt. | ad 
Zudem läßt sich dieser kleine 
Beitrag, den der nicht gemes- 
sene Teil von N(V) zum Integral 1 


V, | 
liefert, aus den | 


Wirkungsquerschnitten Q. fiir 
lonisation berechnen, die den 
jeweiligen StoBspannungen V, 
zugeordnet sind, und deren 
Größe für Hg mehrfach be- 
stimmt wurde und deshalb gut 
bekannt ist. Wir setzen 


2 Volt 


fx arf (yay, 

Abb. 12. Beispiel 
wo die 
ist. Es ist dann 


Vo 


f = 4 


1) Bei dieser Rechnung ist nur das lineare Glied der e-Funktion beriick- 
sichtigt. Hinzunahme des quadratischen Gliedes ändert das Resultat unter 
den gegebenen Bedingungen nur um 1°/,, liegt also innerhalb der Fehler- 
grenzen. 
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Q, wurde den Messungen von Smith!) entnommen. ~ 
Wie bereits oben gesagt, ist die Héhe der einzelnen Maxima 

proportional dem Flächeninhalt unter dem betreffenden Maximum, 

wir kénnen also schreiben 


1 h; 
worin 
h, = Höhe des Verlustmaximums 2 
und 


Sh, = Summe der Höhen aller Maxima unter Einschluß des 
Nullmaximums bedeutet. 


Br Elektronen, die im Stoßraum mit einem Hg-Atom zusammen- 
stoßen, ändern ihre Richtung. Überschreitet diese Richtungsänderung 
einen bestimmten, durch Apparaturabmessung und StoBort gegebenen 
Winkel, so trifft das Elektron auf die Wandung des Stoßraums und 
gelangt nicht auf den Auffänger. Durch diese Streuverluste können 
Fehler entstehen. Denn bei Berechnung der unelastischen Wirkungs- 
querschnitte wird implizite vorausgesetzt, daß alle Elektronen, die in 
den Stoßraum eintreten, auch wirklich nach Durchlaufen des Stoß- 
raums von der Messung im Gegenfeld erfaßt werden. Ist diese Be- 
dingung nicht erfüllt, so wird der unelastische Wirkungsquerschnitt 
zu klein oder groß gefunden, je nachdem die Ausfälle bei den Elek- 
tronen, die Verluste erlitten haben, oder bei den Elektronen, die 
den Stoßraum ohne Energieverlust passiert haben, anteilmäßig größer 
sind. Die Fehler, die durch diese Streuverluste entstehen können, 
werden zwar bei Anregung eines bestimmten Terms um so kleiner, 
je größer die Geschwindigkeit der stoßenden Elektronen wird, weil 
sich ganz allgemein die Richtungsverteilung nach elastischen und 
unelastischen Stößen mit wachsender Geschwindigkeit immer enger 
um die Stoßrichtung herumlegt. Für kleine Geschwindigkeiten sind 
diese Richtungsänderungen aber im Mittel so groß, daß die Meb- 
resultate durch die Streuverluste weitgehend entstellt werden; ordnet 
man also, wie das hier geschehen ist, Stoßraum und Auffänger 
axial hintereinander an, so hat man sich bei der Ausbeutemessung 


PL T. Smith, Phys. Rev. 37. S. 808. 1931. 
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auf Geschwindigkeiten zu beschränken, bei denen sich die durch 
Streuverluste entstehenden Fehler abschätzen lassen und zum Meb- 
wert selbst in einem vernünftigen Verhältnis stehen. Das sind aber 
Geschwindigkeiten, bei denen die kinetische Energie der anregenden 


= Elektronen die Anregungsenergie um ein Mehrfaches übersteigt'). 

ta Die vorliegende Untersuchung beschränkt sich auf Stöße, bei denen 

Dano die Energieverluste höchstens 14 Volt betragen, und bei denen die or j 
na Geschwindigkeit der stoBenden Elektronen zwischen 100 und 300 Volt => 

m, liegt. In diesem Geschwindigkeitsbereich ist die Richtungsverteilung : 


nach elastischen und unelastischen Stößen an Hg von Tate und 
Palmer?) und van Voorhis*) gemessen. Beide Autoren fassen 
allerdings bei der Untersuchung der unelastischen Stöße die Ver- 
luste bis zu 10,2 Volt zusammen. Infolgedessen ist die Richtungs- 
verteilung der unelastischen Stöße, wie sie sich aus den zitierten 
Messungen ergibt, nur eine Mitte- 
es bildung über alle Richtungsvertei- *” 


lungen, die den in der vorliegenden 4+”% be 
Arbeit untersuchten Anregungen | +5 \ 
uf zugeordnet sind. Die Korrektur, 0 \ 
die aus dieser Mittelbildung er- __ \ 
rechnet wird, hängt deshalb auch m? 
nur noch von der Geschwindigkeit 
. des stoBenden Elektrons ab und aie 

ist an alle, einer Kurve entnom- Abb. 13 
2. menen Ausbeuten in gleicher Weise 
‘ anzubringen. - Die Korrekturen, die nach dieser Rechnung an den 
tt gemessenen Werten sasuhringen sind, sind unter gleichzeitiger Be- 
"4 rücksichtigung der durch die Ionisierung bedingten Korrektur in 
Abb. 13 zusammengestellt *). 
a) 1) Es sei ausdrücklich betont, daß diese Beschränkung in keiner Weise 
ny bedingt ist durch den Gedanken, der hier zur Verschärfung der Ausbeute- 
r, messung benutzt wurde. Sie ist vielmehr ausschlieBlich bedingt durch die 
il davon ganz unabhängige Anordnung von Stoßraum und Auffänger, die zu- 
j nächst gewählt wurde, um die Verschärfung der Ausbeutemessung an einem 

möglichst einfachen System zu erproben. Die Beschränkung auf hohe Ge- 
r schwindigkeiten läßt sich durch andere Anordnung von Stoßraum und Auf- 
d fänger grundsätzlich beseitigen. 
k 2) Tate u. Palmer, Phys. Rev. 40. 8. 731. 1932. 
t Be 3) van Voorhis, Phys. Rev. 46. S. 480. 1934. 

4) Der Wechsel des Vorzeichens ist darauf zurückzuführen, daß die Streu- 

j verluste bei elastischen und unelastischen StéBen mit verschiedenen Vorzeichen 
4 eingehen, und bei kleinen Geschwindigkeiten die Wirkung der unelastischen 


Streuung, bei groBen Geschwindigkeiten die Wirkung der elastischen Streuung 
überwiegt. 
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Die mittlere Richtungsverteilung, die der Berechnung der 
Korrektur zugrunde liegt, muß sehr nahe übereinstimmen mit der 
Richtungsverteilung bei Anregung von 2!P, weil Anregung dieses 
Terms bei weitem überwiegt. Dagegen ist die Benutzung dieser 
Richtungsverteilung zur Korrektur der schwachen Ausbeuten mit 
einer Unsicherheit behaftet. Man darf jedoch auf Grund der in 
vielen Einzelheiten geprüften Theorie der Stoßprozesse annehmen, 
daß die Richtungsverteilung bei Anregung der anderen Terme zwar 
ungünstiger, aber nicht so ungünstig wird, wie die Verteilung bei 
elastischer Streuung, d. h. aber, die kleinen Ausbeuten sind höchstens 
um 20°, zu niedrig, Eine größere Sicherheit können aber die 
kleinen Ausbeuten einstweilen sowieso nicht für sich in Anspruch 
nehmen, weil sich bei ihrer Messung schon infolge von Störeffekten 
erhebliche Schwankungen bemerkbar machten. 

Systematische Fehler können ferner dadurch entstehen, daß das 
Spannungsintervall AV notwendigerweise endlich bleiben muß. Zur 
Abschätzung dieser Fehler genügt es, die Kurve J=J(V) im 
Intervall AV durch eine Gerade zu ersetzen. Diese Fehler sind 
rechnerisch leicht zu ermitteln. Anwendung der Taylorschen Ent- 
wicklung bis auf Glieder von mindestens 4. Ordnung in AV ergibt: 


4¥V dv 24 


ir). am”. 


ergibt sich für die Korrektur AN(V), die dem gemessenen Werte 
von N(V) zuzufügen ist, um den richtigen Wert zu erhalten 


dV* 


Wir betrachten zunächst diejenigen Maxima der RER RER 
verteilungskurve, bei denen eine wesentliche Verzerrung durch Über- 
lagerung nicht stattfindet, in unserem Fall also das Maximum der 
Null-Elektronen und das Maximum, das der Anregung von 2!P zu- 
geordnet ist. Da die Form dieser Maxima in gleicher Weise durch 
die Anfangsgeschwindigkeitsverteilung der Elektronen bestimmt ist, 
ist auch der prozentuale Fehler an entsprechenden Stellen überall 
derselbe. Das gilt insbesondere auch für die Höhe h, die allein in 
die Berechnung des Wirkungsquerschnitts eingeht. Ersetzt man zur 
Abschätzung des Fehlers für h den zweiten Differentialquotienten 
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2N 
dV? 
sprechenden Differentialquotienten, der der gemessenen Kurve ent- 
nommen ist, so ergibt sich, daß der Fehler für die Höhe 4°/, sicher 
nicht iibersteigt'). Die Höhen werden also höchstens 4°/, zu klein 
gemessen. Da aber nur die Verhältnisse der Höhen in das Resultat 
eingehen, fällt dieser Fehler völlig heraus, wenn es sich um Maxima 
handelt, die durch Überlagerung nicht gestört werden. Werden 
dagegen die Ausbeutemaxima durch Überlagerung gestört, so wird 
d?N 
‘d V3 
im Verhältnis zur Höhe des Nullmaximums zu groß gemessen. Der 
Fehler kann jedoch, wie aus dem oben gesagten und dem Kurvenbild 
ohne weiteres ersichtlich ist, wenige Prozent nicht überschreiten. Er ist 
also wesentlich kleiner als die Unsicherheit, die den kleinen Ausbeuten 
schon wegen des Streufehlers anhaftet und liegt zudem noch überall 
in entgegengesetzter Richtung. 

Zusammengefaßt läßt sich über die systematischen Fehler 
folgendes sagen: 

Für die Anregung von 2,P sind systematische Fehler höchstens 
im Betrage von 2—3°/, zu erwarten, die kleinen Ausbeuten sind 
wahrscheinlich zu klein ausgefallen, der Fehler kann jedoch auch 
hier 20°/, nicht übersteigen. 

Der Einfluß, den zufällige Fehler auf das Resultat haben können, 
läßt sich abschätzen, wenn man die Streuung der einzelnen Meb- 
punkte um die Kurven betrachtet, die den Gang der Ausbeute- 
querschnitte mit der Voltgeschwindigkeit der anregenden Elektronen 
darstellen. (Abb. 14 im nächsten §.) Für die Anregung von 2!P 
sind diese Fehler sicher nicht größer als + 4°/,, für die kleineren 
Ausbeuten können sie auf + 10°/, steigen. 

Da sich die Anregungsfunktionen der höheren Terme bisher 
nicht einmal größenordnungsmäßig festlegen ließen, darf die Unsicher- 
heit der kleinen Ausbeuten in Kauf genommen werden, zumal die 
Fehlerabschätzung in allen Fällen äußerste Grenzen gibt und der 
Gesamtfehler wahrscheinlich wesentlich kleiner sein wird. 

Es wurden auch eine Reihe von Punkten in dem Bereich 
zwischen 40 und 100 Volt gemessen. Da jedoch die Korrekturen 
für diese Werte durch Extrapolation ermittelt werden mußten, und 
die Streuung der Meßwerte infolge der geringen Stromwerte erheblich 
wurde, kann diesen Messungen kein wesentliches Gewicht gegeben 
werden. 


), der wahren Geschwindigkeitsverteilung durch den ent- 


an der Stelle der Höhe algebraisch größer, es wird also auch h 


1) Bei den Messungen war 4 V etwa = 0,2 Volt. 
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§5. MeBergebnisse 

Zunächst seien die in allen Kurven festgestellten Maxima den 
Termen zugeordnet, deren Ausbeute sie messen. 

Die beiden Maxima bei 5,4 und 4,9 Volt sind der Anregung 
der 2®P-Terme zuzuordnen. Anregung von 2°P, und 2°P, (4,66 
und 4,86 Volt) wird nicht getrennt. 

Das Maximum bei 6,7 Volt entspricht eindeutig der Anregung 
von 2'!P. 

Es folgt ein sehr kleines und stark gestörtes Maximum bei 
7,9 Volt. Die Lage entspricht sehr gut dem 2 !S-Term. (Anregungs- 
spannung 7,89 Volt.) Wäre Anregung von 2°S (Anregungsspannung 
7,69 Volt) wesentlich beteiligt, müßte die Verschiebung — besonders 
bei kleineren Geschwindigkeiten, wo die Ausbeute größer ist — 
merklich geworden sein. 

Die Einordnung der beiden nächsten Maxima würde bei der 
dichten Lage der Terme Schwierigkeiten machen, wenn nicht eine 
Anomalie der !P-Terme die Zuordnung eindeutig vorschreiben würde. 
Bei dem Verlust von 8,8 Volt ist außer der Anregung von 3!P, 3°D 
nicht auszuschließen, wenn man bei der Deutung lediglich von der 
Größe des Energieverlustes ausgeht. Das nächste Maximum liegt 
aber nicht bei 9,5 Volt, wo sich auf ganz engem Bereich 4!D, 4°P, 
4°D und 4°F zusammendrängen, sondern liegt bei 9,8 Volt, wo sich 
infolge der Anomalie der 'P-Terme 4'!P und 5'P dicht aneinander- 
lagern. (Anregungsspannung bei 9,71 und 9,86 Volt.) Es muß sich 
also auch bei dem Maximum, das bei 8,8 Volt erscheint, um die 
Anregung von 3!P handeln. 

Es wird also, wie das auch theoretisch zu erwarten ist, vorzugs- 
weise die Termreihe angeregt, die bei gleicher Multiplizität unmittelbar 
mit dem Grundterm kombiniert. Bei Anregung der Interkombinations- 
linie 119 — 23P orientiert sich die Anregung offenbar nicht nach 
den Auswahlregeln für die inneren Quantenzahlen, sondern nach 
dem Gewicht der einzelnen Tripletterme. 

Das Maximum bei 11,0 Volt entspricht Ionisierungsverlusten. 
Bei den ionisierenden Stößen liegen die Verhältnisse insofern anders 
als bei der Anregung scharf definierter Terme, als hier die Ausbeute 
kontinuierlich über die Voltskala verteilt ist. Das Verhalten der auf 
die Skaleneinheit bezogenen Ausbeute beim Übergang über die 
Ionisierungsspannung ist in quantitativ auswertbarer Form bisher 
nicht gemessen. Das scharfe Heraustreten des Ionisierungsmaximums 
zeigt jedenfalls, daß die auf die Einheit der Voltskala bezogene 
mittlere Ausbeute an anregenden Stößen sich an der Ionisierungs- 
grenze nicht monoton fallend in der ebenfalls auf die Skaleneinheit 
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bezogenen lIonisierungsausbeute fortsetzt. Die Verschiebung des 
Maximums zeigt vielmehr, daß diese auf die Skaleneinheit bezogene 
Ausbeute unmittelbar hinter der Ionisierungsgrenze ansteigt, um dann 
allerdings schon bei etwa 11,5 Volt zu unmerklichen Werten ab- 
zusinken. Es scheinen hier also ähnliche Verhältnisse vorzuliegen, 
wie sie Ditehburn!) beim Übergang von Absorptionsserien ins Serien- 
kontinuum beobachtete. Ditchburn fand, daß die Absorption an 
der Seriengrenze beim Übergang in das Endkontinuum zunächst 
nahezu unstetig zunahm, um erst dann wieder monoton mit wachsender 
Frequenz abzunehmen. Ob allerdings in dem vorliegenden Fall die 


cn? - 
ns 
ip 
> 
4 .. 
3 
2 
yo 8 120 10 200 240 280 320 360 400 
Yo —> Volt 
fiir die Verluste 5,4, 6,7, 9,8, 11,0 Volt. Ihe We: 
=--- Kurven ohne Korrektur 
Kurven nach den in $4 besprochenen Korrekturen si . 


auf die Skaleneinheit bezogene Ausbeute an der Jonisierungsgrenze 
tatsächlich eine Unstetigkeit hat, läßt sich bei der geringen Auflösung 
noch nicht feststellen, ist aber nicht ausgeschlossen. 

Da Messungen mit einer wesentlich verbesserten Apparatur in 
Aussicht genommen sind, wird zunächst auf die Auswertung der 
Maxima bei 4,9 sowie bei 7,9 und 8,8 Volt verzichtet, weil die 
Zerlegung der N(V)-Kurve am Ort dieser Maxima infolge der einst- 
weilen noch geringen Auflösung so unsicher war, daß die Mitteilung 
quantitativer Ergebnisse nicht gerechtfertigt erscheint. Für die 
übrigen Terme, für die die Fehlerdiskussion schon oben gegeben 


1) R. W.Ditehburn, Ztschr. f. Phys. 107. S. 719. 1934. 
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wurde, wird der Anregungsquerschnitt und seine Abhängigkeit von 
der Elektronengeschwindigkeit in Abb. 14 dargestellt. 

Ein quantitativer Vergleich ist nur bei den Messungen von 
Tate und Palmer (a. a. O.) möglich, die durch Integration über 
die von ihnen gemessene Richtungsverteilung den Wirkungsquerschnitt 
für die Gesamtheit der anregenden Stöße bestimmen!) Den von 
Tate und Palmer gemessenen Wirkungsquerschnitt erhalten wir 

durch Addition der hier gemessenen 
um? Einzelquerschnitte. Dabei wurden 
cm . 
0 \ die Wirkungsquerschnitte fiir 4,9 
2 ‘ bzw. 7,9 und 8,8 Volt, die im 
|, 7 ganzen nur eine kleine Korrektur 
N ausmachen, nach den N(Vj-Kur- 
uw ven geschätzt. Die Abweichungen 
oy liegen bis 300 Volt ganz wesent- 
>> lich unter den oben zugelassenen 
Fehlergrenzen, nur für die größten 
Geschwindigkeiten gehen sie etwas 


5 > darüber hinaus. Da bei Tate und 
Abb. 15 Palmer für deren Messungen 
er keine Fehlergrenzen angegeben 


sind, läßt sich nicht ohne weiteres 
entscheiden, worin diese Diskrepanz mit den beiderseitigen Fehler- 
grenzen noch vereinbar ist. 

Zu dem Maximum mit einer Verlustspannung von etwa 1,5 Volt 
ist zu sagen, daß bisher eine physikalische Deutung noch nicht 
gelungen ist. Da das Maximum, wie bereits erwähnt, bei Ausfrieren 
des Hg-Dampfes verschwindet, ist jedoch anzunehmen, daß es sich 
nicht um Störeffekte handelt, die durch die Apparatur bedingt sind. 


Zusammenfassung 


i cs a Die Empfindlichkeit bei der Bestimmung der Ausbeute un- 


elastischer StéBe wird dadurch sehr gesteigert, daB die Gegen- 
spannungscharakteristik durch Uberlagerung einer kleinen Wechsel- 


1) Da die Messungen von Tate und Palmer auch zur Berechnung der 
Korrektur herangezogen wurden, sei ausdriicklich vermerkt: Der Vergleich 
gibt insofern eine völlig unabhängige Prüfung des Resultates, als er zeigt, daß 
die Ausbeute, die in den bei der vorliegenden Messung unmittelbar erfaßten 
Raumkegel hineinfällt, mit dem entsprechenden Ausbeuteanteil, wie er sich 
aus den Messungen von Tate und Palmer berechnet, gut übereinstimmt. 
Zudem sind oberhalb von 200 Volt die Korrekturen so klein, daß von einer 
Abhängigkeit schon deshalb nicht die Rede sein kann. 
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spannung elektrisch differenziert wird, und damit die Geschwindigkeits- 
verteilung des austretenden Elektronenstrahls unmittelbar gemessen 
wird. Es werden elektrische Anregungsfunktionen einzelner Quecksilber- 
terme zwischen 100 und 380 Volt aufgenommen und die Fehlerquellen 
der Methode diskutiert. 


Diese Arbeit wurde im Physikalischen Institut der Technischen 
Hochschule Hannover ausgeführt. Herrn Prof. Dr.H. Bartels spreche 
ich für die Anregung zur vorliegenden Untersuchung und deren 
dauernde Unterstützung meinen herzlichen Dank aus. Der Helmholtz- 
Gesellschaft habe ich für weitgehende Unterstützung zu danken. 


Hannover, Physikalisches Institut der Technischen Hochschule. 


(Eingegangen 15. März 1939) 
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- Quantenmechanische Theorie 


der Réntgenstrahlinterferenzen in Kristallen 


eh: I, Ableitung und allgemeine Diskussion 
der dynamischen Grundgleichungen 


Von G. Moliére 
(Mit 2 Abbildungen) 


In den ursprünglich von P. P. Ewald und C. G. Darwin’) 
entwickelten Formen der dynamischen Theorie der Kristallréntgen- 
interferenzen wurden die einen Kristall aufbauenden Atome oder 
Ionen als punktförmige Streuzentren betrachtet. Demgegenüber beruht 
die später von M. v. Laue?) entwickelte Form der Theorie auf der 
Annahme eines kontinuierlich im Kristall verteilten Streuvermögens 
entsprechend der kontinuierlichen Ladungsverteilung der Elektronen 
nach der Wellenmechanik. Die kontinuierliche Verteilung des Streu- 
vermögens wird dabei in der Weise berücksichtigt, daß der Kristall 
als Dielektrikum mit dreifachperiodischer Dielektrizitätskonstante 
betrachtet wird, in dem sich die Röntgenwellenfelder nach den 
Maxwellschen Gleichungen des Dielektrikums ausbreiten. M. Kohler‘) 
hat eine quantenmechanische Begründung des v. Laueschen Ansatzes 
gegeben. Er berechnete die Stromdichte %, die in einem Kristall 
durch ein elektromagnetisches Wechselfeld induziert wird, indem er 
letzteres als Störungsglied in die Schrödingergleichung der Kristall- 
elektronen einsetzte. Er fand dabei einen Ausdruck für die Strom- 
dichte, welcher in zwei Anteile %” und %® zerfällt, deren erster für 
die Streuung von Röntgenstrahlen genügend hoher Frequenz, denen 
gegenüber die Elektronen praktisch als frei betrachtet werden können, 
allein maßgeblich ist, während der zweite „Resonanz“-Anteil im 
Bereich der anomalen Dispersion wesentlich ist. Bei Berücksichtigung 


1) C. G. Darwin, Phil. Mag. 27. S. 675. 1914; P. P. Ewald, Ann. d. 
Phys. 54. S. 519 u. 557. 1917; P.P. Ewald, Ztschr. f. Kristallogr. 97. S. 1. 1937. 
2) M. v. Laue, Erg. d. exakt. Naturw. Bd. X. S. 133. 1931. Vgl. auch 
E. Lohr, Wiener Ber. 133. S. 157. 1924. 
3) M. Kohler, Berliner Ber. phys.-math. Kl. XIX. $. 1. 1935. Vgl. auch 
F. Galperin, Phys. Ztschr. d. Sowjetunion 11. S. 182. 1937. Galperin führt 
im wesentlichen die gleiche quantenmechanische Rechnung wie M. Kohler 
durch. Seine Schlußfolgerung, man müsse in der v. Laue schen Theorie eine 
dreifachperiodische Dielektrizitätskonstante in Form eines nichthermitischen 
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allein des Anteils %® konnte Kohler den v. Laueschen Ansatz 
einer dreifachperiodischen Jielektrizitätskonstante und damit die 
v. Lauesche Form der dynamischen Grundgleichungen bestätigen. 
In der vorliegenden Arbeit, die an M. Kohlers Arbeit anknüpft, 
wird gezeigt, daß auch der Resonanzanteil %® der induzierten Strom- 
dichte vollständig und in geschlossener Form in die dynamische 
Theorie der Kristallröntgeninterferenzen eingebaut werden kann. 
Dabei ergeben sich die v. Laueschen dynamischen Grundgleichungen 
in einer verallgemeinerten Form [vgl. Gl. (2,13)], in der sie auch für 
den Bereich der anomalen Dispersion gültig sind und dort die 
Ausbreitung — d.h. Brechung, Reflexion usw. und zugleich auch 
die Absorption — der Röntgenstrahlen im Kristall beschreiben. Die 
in den Grundgleichungen auftretenden Koeffizienten können im 
Bereich der anomalen Dispersion nicht mehr als Fourierkoeffizienten 
eines gewöhnlichen Streuvermögens in Form einer einfachen Orts- 
funktion aufgefaßt werden, vielmehr tritt an die Stelle dieses Streu- 
vermögens ein tensorieller Integralkern [vgl. Gl. (1,33)], welcher die 
an einer Stelle r induzierte Polarisation ®(r) mit der Feldstärke E(r') 
an verschiedenen Stellen r‘, über die integriert wird, verknüpft. 
Der Einfluß der Resonanzstreuung bei den Röntgenstrahlinter- 
ferenzen wurde bereits von H. Hönl!) berücksichtigt, jedoch nicht 
im Rahmen einer in sich geschlossenen dynamischen Theorie, sondern 
durch Berechnung eines verallgemeinerten Atomfaktors auf Grund 
der Wallerschen Dispersionsformel für das freie Atom?) Der 
Begriff des Atomfaktors stammt ja aus der elementaren Theorie 
der Röntgenstrahlinterferenzen und ist der dynamischen Theorie 
P.P. Ewalds nur künstlich aufgepfropft. Praktisch vernachlässigt 
Hönls Theorie den Einfluß der Gitterstruktur, der gerade bei der 
Resonanzstreuung merkbar wird und zu den bekannten von R. de 
L. Kronig*) gedeuteten Anomalien des Absorptionskoeffizienten sowie 
zu röntgenoptischer .\nisotropie führt, wie im anschließenden Teil II 
dieser Arbeit gezeigt wird. — An die Stelle des Atomfaktors treten 
in v. Laues bzw. unserer Theorie die in den Grundgleichungen 
auftretenden Koeffizienten; für die letzteren erhalten wir Formeln, 
die denen von H. Hönl bzw. I. Waller sehr ähnlich sind. [Vgl. 
z. B. die Gl. (4,12) und (4,13) dieser Arbeit.] Der unserer (sowie 
bereits M. Kohlers) Theorie zugrunde liegende Ansatz kann als 
Übertragung des kohärenten Streuanteils der Wallerschen Dispersions- 
formel vom einzelnen Atom auf das Elektronensystem des gesamten 


1) H. Hönl, Ann. d. Phys. 18. S. 625. 1933. em 


~ 


2) I. Waller, Ztschr. f. Phys. 51. S. 213. 1928. 
3) R. de L. Kronig, Ztschr. f. Phys. 70. S. 317. 1931; 75. 8.191 u. 468. 1932. 
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Kristalls aufgefaßt werden. Ein wesentlicher Unterschied besteht 
jedoch darin, daß in unserer Theorie» nicht zwischen einfallender 
und gestreuter Welle unterschieden wird; vielmehr wird das gesamte 
elektromagnetische Wechselfeld so bestimmt, daß es zusammen mit 
der von ihm induzierten Polarisation den Maxwellschen Gleichungen 
genügt. Dieses Verfahren ist charakteristisch für v. Laues 
dynamische Theorie. 

Der vorliegende Teil I dieser Arbeit enthält die Aufstellung der 
Grundgleichungen sowie eine eingehende Diskussion der darin auf- 
tretenden Koeffizienten. Der Inhalt der einzelnen Paragraphen ist 
folgender: §1 enthält die zugrundeliegende quantenmechanische 
Störungsrechnung; im §2 werden die verallgemeinerten Grund- 
gleichungen aufgestellt. Nach allgemeinen Betrachtungen über inko- 
härente Streuung und Absorption im §3 dienen die folgenden Para- 
graphen der eingehenden Diskussion der in den Grundgleichungen 
auftretenden Koeffizienten, und zwar werden im $4 bestimmte Aus- 
drücke für dieselben durch Einführung Hartreescher Näherung 
(d.h. Reduktion auf das Einelektronenproblem) hergeleitet, die im 
$5 in Streu- und Absorptionsanteil unterteilt werden. Im §6 wird 
der Tensorcharakter und im letzten §7 die Frequenzabhängigkeit 
der Koeffizienten diskutiert. — Im Teil II dieser Arbeit!) werden 
die Folgerungen für die eigentliche dynamische Theorie der Brechung, 
Reflexion und Absorption von Röntgenstrahlen in Kristallen behandelt. 


§1. Quantenmechanische Berechnung der durch ein 
elektromagnetisches Wechselfeld im Kristall induzierten Polarisation 


Die ungestörte Schrödingergleichung der Kristallelektronen habe 


die Form 

und ihre Lösungen seien is 
wobei »,, = =. und W,, die Energie des m-ten stationären Zustandes 


ist. Durch ein elektromagnetisches Wechselfeld tritt ein kleines 
Zusatzglied H, (t) zur Hamiltonfunktion H, hinzu, so daß die 


gestérte Hamiltonfunktion lautet 


H, (t) sei mit der Frequenz v sili veränderlich, also von der Form 
H, = H+ e2 2zirt H- 
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Die stationären Lösungen des gestörten Problems ergeben sich dann in 
bekannter Weise durch den Ansatz 
(1,5) ys — Um tivmt rs ut +r)t ti (um — 


Entwickelt man gemäß 


(1,6) ut = > b> 
n 


so ergibt Gl. (1,3) fiir die Koeffizienten b* und b-: . 


> 
(1,7) b* = 
wobei 


ae) (W, — Wr) 
(1,9) Yam = Yn — Yn = h 


ist. Speziell für diejenige Lösung ¥Y von (1,3), die im Grenzfall verschwindender 
Störung in den Grundzustand ¥, des ungestörten Problems übergeht, erhält 
man nach (1,5), (1,6) und (1,7) mit m = 0 einen Ausdruck der Form 


2rint 


(1,10) y= "+2 F, e = 


( ) of ) h v) h (rot v) 


Das elektromagnetische Wechselfeld sei durch sein Vektor- 
potential 2 beschrieben, das so normiert sei, daß das zugehörige 
skalare Potential verschwindet. Das Störungsglied H, # in Gl. (1,3) 
hat dann nach M. Kohler!) die Form 


A 


(1,12) a, Fo > (div, (U, #) + (A, grad.) 7]. 

Dabei charakterisiert der Index s die Koordinaten des s-ten Elektrons . 


und W, bedeutet den Wert des Vektorpotentials Y\ am Ort des s-ten 
Elektrons. Man beachte, daß das Vektorpotential in unserem Fall 
nicht divergenzfrei ist. Mit div 21 = 0 würde Gl. (1,12) die von 2 
andern Dispersionsproblemen her bekannte Gestalt annehmen. — 
Das weitere Programm besteht darin, die Strom- und Ladungsdichte 
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\ 
(1,13,2) el) = > | fe war 


zu berechnen’), Hierbei bedeutet e die Ladung und m die Masse 
des Elektrons, grad, den Gradientoperator bezüglich der Koordinaten f 
des k-ten Elektrons und dr,’ das Produkt der Volumenelemente in 
simtlichen Elektronenkoordinaten auBer dem k-ten. Das Vektor- \ 
potential X soll dann nachträglich so bestimmt werden, daß es zu- 
sammen mit der Stromdichte § die Maxwellschen Gleichungen er- 
füllt. Das beschriebene Verfahren ist insofern nicht ganz streng, 
als die Größe W, in (1,12) und (1,13,1) eigentlich den Wert des 
Vektorpotentials am Orte des s-ten Elektrons abzüglich des vom 
s-ten Elektron selbst erzeugten Anteils von % bedeuten sollte. 
W, sollte also in der Ausdrucksweise P. P. Ewalds das „erregende 
Feld“ für das s-te Elektron sein. Der Unterschied zwischen dem 
erregenden Feld und dem Gesamtfeld ist jedoch im Falle von 
Röntgenstrahlen zu vernachlässigen. 


Wir wollen nun die Stromdichte % berechnen, indem wir Gl. (1,10) in 
(1,13,1) einsetzen. Wir zerlegen dazu § in zwei Anteile, 
= 34+ 9” 


indem wir zunächst den X enthaltenden Anteil in (1,13,1) vorwegnehmen: 


Da %, von den Koordinaten des k-ten Elektrons abhängt, über die in (1,14) 
nicht integriert wird, kann es vor das Integralzeichen und sodann wegen 
t= r vor das Summenzeichen gezogen werden, so daß sich mit Benutzung 
von (1,13,2) ergibt: 


eo! 
(1,15) 3° =- Tann. 


Da 4 klein ist, kann in (1,14) in erster Näherung Y durch Y= u, gates 
ersetzt werden, so daß also g in Gl. (1.15) die Ladungsdichte des ungestörten 
Kristalls darstellt. § nach Gl. (1,15) stellt, wie wir sehen werden, denjenigen 
Anteil der Stromdichte % dar, der bei genügend hohen Frequenzen allein mab- 
geblich ist und der, wie bereits M. Kohler gezeigt hat, zum Ansatz einer 
dreifachperiodischen Dielektrizitätskonstante im Sinne der v. Laueschen 


Theorie führt. Der zweite Anteil der Stromdichte in (1,13,1) lautet: x 
(1,16) @) = > If W grad, — W* grad, 


(1,24 


1) Vel M. Kohler a.a.0.; die dortige Gleichung für % enthält ver- 
sehentlich einen Vorzeichenfehler. a 
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Setzt man hier % nach Gl. (1,10) ein und berücksichtigt dabei nur die 
von der Störung herrührenden Glieder, und zwar insofern sie klein von erster 
Ordnung sind, so erhält man unter Benutzung der Strommatrixelemente A 


11) „= inn™ grad, u * —u, grad, u,) 


für den Anteil 3°) der HE 


(1,18) = > {ion + F, o}. 
n 

Schließlich sind noch die Koeffizienten F,, (¢) nach Gl. (1,11) zu berechnen. 

Das Vektorpotential X soll mit der Frequenz » zeitlich oszillieren, also von der 

Form sein 


(1,19) Y= + 


Setzt man dies in (1,12) ein und vergleicht mit (1,4), so erhält man 
0,201) Hr ** D faiv, (a + (grad, 


Snime 


und 
— eh 1% 

(1,20,2) = > { div, 7) + verad, 
8 = 

Für die Matrixelemente Hy, ergibt sich somit: 
ag) {2 div, (%, uo) + 


wobei dr das Produkt der Yuan in sämtlichen Elektronenkoordi- 
naten bedeutet. Das erste Glied des Integranden läßt sich umformen: 

(1,22) u, div, (W,u,) = div, (u, U, uo) — up grad, 

Das von der Divergenz herrührende Oberflächenintegral verschwindet bei 


der Integration über den ganzen Kristall, und wir erhalten mit Benutzung 
von (1,17): 


f 
(1,23,1) H,. (A (t’) )) dt 
In entsprechender Weise ergibt sich: ne fe: 
(1,23,2) (v’)) dr’. 


Dabei bedeutet dr’ das Volumenelement im gewöhnlichen dreidimensionalen 
Raum. Setzt man H,* 9 und H,, nach (1,23,1) und (1,23,2) in den Ausdruck (1,11) 
fiir die Koeffizienten 7 () und diesen sodann in GI. (1,18) ein, so erhält man 
für den Anteil 3 der Stromdichte den Ausdruck 


32 «p= 2c ion) dr’ inn JC dr e~2aivt 


9. hiv. v) 
(1,24) 
Ä + Je erirt ion tt) ion) Qaivt}, 
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Gl. (1,24) ist, wie es sein muß reell, da sich, wie man sieht, jeweils ER 
komplexe Glieder entsprechen. 

Wir gehen nunmehr für alle Feldgrößen zur komplexen Schreib- 
weise über, indem wir an Stelle von (1,19) den komplexen Ansatz 
machen: 

(1,25) Art) = 

sowie entsprechende Ansätze für die Stromdichte %, die elektrische 
Feldstärke € usw., und setzen fest, daB physikalische Bedeutung 
jeweils nur dem Realteil der betreffenden FeldgréBe zukommen soll. 
Diese Schreibweise ermöglicht es, Phasenverschiebungen einfach 
durch komplexe Faktoren zu berücksichtigen, wodurch sich unsere 
Gleichungen wesentlich vereinfachen werden. Gl. (1,15) bleibt beim 
Übergang zur komplexen Schreibweise unverändert. In Gl. (1,24) 
dagegen dürfen wir die W'* enthaltenden Glieder durch ihre konjugiert- 
komplexen ersetzen, da hierbei nur der Imaginärteil von X® ver- 
ändert wird. Unter Berücksichtigung von (1,25) nimmt dabei Gl. (1,24) 
die Gestalt an: 


ion dr’ ion ® (t',2), ig g(t’) dr’ | 


(1,26) ) It 1.23 > v) + h(rno af ») 


Damit haben wir in den Gl. (1,15) und (1,26) handliche Aus- — 


drücke fiir die beiden Anteile der Stromdichte % erhalten, die durch 
ein (durch sein Vektorpotential WA beschriebenes) Wechselfeld im 
Kristall induziert wird. Wir wollen diese beiden Gleichungen noch 


umformen, indem wir an Stelle von % die elektrische Feldstärke 


- 1 ax 2niv 
(1,27) gre- — 
und an Stelle von § die durch die Gleichung _ fon ene 
(1,28) = SF = 


definierte „dielektrische Polarisation“ ® einführen. Die Maxwell- 
schen Gleichungen nehmen bei dieser Substitution die gleiche Form 
an, wie in der makroskopischen Elektrodynamik des Dielektrikums, 
nämlich 


| rot E=— w9=0, 


wobei D=E+4n® gesetzt ist, und an Stelle der beiden An- 
teile %» GI. (1,15) und §” Gl. (1,26) der Stromdichte erhalten wir 


es Du 


(1,32) 
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zwei Anteile PB’ und $ der Polarisation ®B. Für den ersteren ergibt 
sich dabei eine Gleichung der Form = 
eine dreifachperiodische Suszeptibilität darstellt. Bei genügend 
hohen Frequenzen, wo der zweite Anteil %” der Polarisation zu 
vernachlässigen ist, beschreibt Gl. (1,30) zusammen mit (1,31) und 
den Maxwellschen Gl. (1,29) die kohärente Streuung von Röntgen- 
strahlen, und zwar ist dies die klassische Streuung freier Elektronen, 
deren räumliche Verteilung der wellenmechanischen Ladungsdichte 
entspricht. Dies ist das bereits von M. Kohler gefundene Ergebnis’). 
Für den zweiten Anteil der Polarisation ergibt sich nach Gl. (1,26) 


n 


Diese Gleichung läßt sich formal in folgender Weise schreiben: 


(1,33) 4a = Eq 


Dabei ist die Integration über das gesamte Kristallvolumen V zu 
erstrecken. Die Größe A (t,r‘) stellt eine Art von Integralkern mit 
Tensorcharakter dar: 


hivng —”) +») 


wobei die Klammerausdrücke der Form {X; 8} Tensoren (dyadische 
Produkte) mit den Komponenten ®},, = bedeuten. 

Der Anteil % nach Gl. (1,33) ist, ‘wie sich zeigen wird, bei 
Röntgenfrequenzen nahe den Absorptionskanten wesentlich. Wie 
man sieht, ist es in diesem Frequenzbereich nicht mehr möglich, 


1) Die Suszeptibilität ist nach Gl. (1,30) hier so definiert, daß sie 4n ® 
anstatt wie sonst üblich ® mit der Feldstärke € verknüpft. Dadurch wird in 
den Endformeln das Auftreten lästiger Faktoren 4 vermieden. Wir betonen, 
daß in dieser Arbeit das übliche elektrostatische und nicht das Lorentzsche 
Maßsystem benutzt wird; es erscheinen also in den Formeln alle Größen, 
insbesondere auch die elektrische Elementarladung e in den üblichen elektro- 
statischen Einheiten. e bedeutet in dieser Arbeit (außer die Basis der natür- 
lichen Logarithmen) die Ladung des Elektrons, also „minus Elementarladung“. 

2) Vgl. M. Kohler, a.a. O. 
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den Kristall einfach durch eine dielektrische Suszeptibilität zu 
charakterisieren, welche wie in Gl. (1,30) die Polarisation ® an einer 
Steller mit der elektrischen Feldstärke € an der gleichen Stelle r 
verknüpft; vielmehr hängt nach Gl. (1,33) die Polarisation ®(r) an 
einer Stelle r von der elektrischen Feldstärke E(r') an verschiedenen 
Stellen r’, über welche integriert wird, ab. (Strenggenommen kommt 
es dabei auf den Verlauf der Feldstärke € im ganzen Kristall an, 
praktisch jedoch, wie sich zeigen wird, jeweils nur auf einen kleinen 
Bereich von atomarer Größenordnung in der Umgebung der jeweiligen 
Stelle r.) Gl. (1,33) stellt also eine wesentliche Erweiterung des 
v. Laueschen Ansatzes dar, die jedoch, wie im folgenden Paragraphen 
gezeigt wird, nur eine einfache formale Änderung der v. Laueschen 
Grundgleichungen zur Folge hat. 


$2. Ableitung der modifizierten Grundgleichungen 
der v. Laueschen dynamischen Theorie der Röntgenstrahlinterferenzen 


Wir schreiben nach v. Laue?) die Feldvektoren als 3fache 
Fourierreihen (genauer: durch 3fache Fourierreihen modulierte 


— 227i (Kg, ziv 


und entsprechende Ausdrücke für die übrigen Feldgrößen ®,D usw. 
Dabei sind die Ausbreitungsvektoren 8, gemäß 


b, 


miteinander verknüpft, wobei b, die reziproken Gittervektoren sind. 
Die Möglichkeit des Ansatzes (2,1) ergibt sich aus der Tatsache, 
daß die Gl. (1,29), (1,30) und (1,33) linear und homogen in den Feld- 
größen sind, zusammen mit den Periodizitätsbedingungen 


wobei a einen Gittervektor bedeutet. — Weiter entwickeln wir die 

der Gl. 31) in eine 3fache Fourierreihe: 
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aus a 31) zu berechnen sind. Die Integration ist dabei über die 
Gitterzelle A zu erstrecken. Mit den Entwicklungen (2,1) und (2,2) 
folgt aus Gl. (1,30) die Beziehung zwischen den Fourierkomponenten 
der Feldgrößen ® und €: 


(2,4) 4 = -q E, . 
& 


Ferner multiplizieren wir die Fourierreihe (2,1) für B®: 


mit e?7i(%»,t) und mitteln his Ergebnis räumlich. Es würde dabei 
genügen, die Mittelung über eine Gitterzelle zu erstrecken, doch 
wollen wir, um die Endgleichungen symmetrischer zu gestalten, über 
das ganze Kristallvolumen V mitteln. Wegen — 8, = bn—, 
verschwinden dabei alle Glieder außer dem mit g=m, so daB © 


sich ergibt the! 
(t dr To 

Setzt man hier §§(r) nach GI. (1,33) und anschließend fir G(r’) die = 
F (2,1) so ergibt sich 
ba 
Durch Vertauschung der Reihenfolge von Summation und Integration : 
erhält diese Gleichung die Form “% 

26) Ang = de de | 3 
vv 
oder mit Benutzung des Ausdrucks (1,34) für A(t,r‘) . u 


1 {ion (ts jon(t)} 22i(fK - Rar’) , 


ist. Es sei sogleich bemerkt, daß die Koeffizienten A, nicht vom 


Kristallvolumen V abhängen, wie weiter unten gezeigt ‘werden wird. t 
Für die Fourierkomponenten der gesamten Polarisation + 
ergibt sich nach (2,4) und (2,5): 
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wobei zur Abkürzung 


mq 


gesetzt ist. Nach Gl. (2,7) sind die Koeffizienten A, und damit 
auch die C,, im allgemeinen Tensoren. — Es ist zweckmäßig, an 
Stelle der elektrischen Feldstärke & die dielektrische Verschiebung 
D=€+4a$ einzuführen, weil wegen dv®=0 die Fourier- 
komponenten D, senkrecht auf den zugehörigen Ausbreitungs- 
vektoren St, stehen, d.h. die Partialwellen von D rein transversal 
sind. Da bei Réntgenstrahlen D und € sehr nahezu überein- 
stimmen, dürfen wir in (2,8) einfach die ©, durch ®, ersetzen, so 
daß wir erhalten 


[Genau genommen wire in (2,10) C,,, durch eine Reihe Ei = 


Ong= Cng - 


zu ersetzen, doch ist dies praktisch belanglos, da die C,, für 
Röntgenstrahlen von der Größenordnung 10% sind.) 

Schließlich sind noch die Maxwellschen Gleichungen (1,29) in 
Fourierkomponenten umzuschreiben. Eliminiert man aus (1,29) die 
magnetische Feldstärke 9 und ersetzt € durch D— 4a und setzt 
in die entstehende Beziehung zwischen ® und D die Fourier- 
reihen (2,1) ein, so erhält man für die Fourierkomponenten die 


wobei 
v 
(2,12) k=— und Bn Rn 
2 


m 


den Anteil des Vektors ®, bedeutet, der senkrecht auf &,, steht. 
Setzt man ®, nach (2,10) in (2,11) ein, so ergibt sich }) 


1) Eliminiert man nicht die ®, sondern statt dessen die D,, aus (2,10) 
und (2,11), so erhält man eine andere Form der Grundgleichungen, die der- 
jenigen ähnlich ist, dieP.P. Ewald, der nach M. v. Laues Vorgang ebenfalls 
den zum verteilten vollzogen hat, 
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Dies sind die verallgemeinerten Grundgleichungen der v. Laueschen 
dynamischen Theorie der Röntgenstrahlinterferenzen. Ausdrücklich 


j sei betont, daß die Koeffizienten C,,, sich im allgemeinen nicht als 

- Fourierkoeffizienten einer dreifachperiodischen Suszeptibilität auf- 

= fassen lassen. Nur im Fall von Röntgenstrahlen genügend hoher 

8 Frequenz, nämlich bei genügender Entfernung von den K-Kanten 

"i der am Gitteraufbau beteiligten Atome, können die Größen A,, in 

ns Gl. (2,9) vernachlässigt und die ©, somit durch die Fourierkoeffi- 

al zienten %„_, nach Gl. (2,3) ersetzt werden. In diesem Fall nimmt 

si die Gl. (2,13) die ursprünglich von M. v. Laue eingeführte Form der 
Grundgleichungen an ’). 

§ 3. Ergänzungen auf Grund der Diracschen Strahlungstheorie: “7 
inkohärente Streustrahlung und Absorption ws 

Die im § 1 dieser Arbeit benutzte Methode der quantenmecha- 

nischen Störungsrechnung, bei der das Röntgenwellenfeld als Stérungs- 

potential in die Schrédingergleichung der Elektronen eingesetzt wird, 

5 liefert (jedenfalls in der von uns durchgeführten ersten Näherung) 

nur die kohärente Röntgenstreuung. Die inkohärente (Fluoreszenz-) 

: Streustrahlung geänderter Frequenz erhält man mit, wenn man statt 

. dessen die Diracsche Strahlungstheorie benutzt. Die Anwendung 

t der Diracschen Strahlungstheorie auf das Problem der Streuung 

| von Röntgenstrahlen in Kristallen beabsichtigt der Verf. in einer 

; späteren Arbeit, die sich speziell mit der inkohärenten Röntgen- 


streuung befassen soll, ausführlich zu behandeln. 

Wir gehen an dieser Stelle überhaupt darauf ein, weil sich nach 
der Strahlungstheorie zugleich gewisse Ergänzungen unserer im § 1 
erhaltenen Formeln für die kohärente Streuung — nämlich gewisse 
Zusatzglieder in den Resonanznennern — ergeben, auf die es uns \ 
hier ankommt. Das Auftreten inkohärenter Streustrahlung hat 4" 
nämlich Strahlungsdämpfung und infolgedessen Absorption der j 
Primärstrahlung zur Folge. Wir erinnern in diesem Zusammenhang 
an die klassische Dispersionstheorie, nach der im Resonanznenner 
im Ausdruck für die Polarisierbarkeit eines Atoms (d. h. dort: eines 
„elastisch gebundenen Elektrons“) ein imaginires „Dämpfungsglied‘“ 
auftritt als Folge davon, daß dieses Atom Streustrahlung emittieren 
kann. Genau das gleiche ergibt sich nach der Diracschen Strah- 
lungstheorie, und zwar ergibt diese, daß die Resonanznenner 


al 


- 


cal 


1 


iF Ply, 
{ 
- 
= 
_ 
1) Vgl. M. v. Laue, a.a. O. ENDEN a Er 
uf 
X 


in dem Ausdruck A (t,t’) unserer Gl. (1,34) (der ja eine verallgemei- 
nerte Polarisierbarkeit des Kristalls darstellt) durch solche der Form 
h (Yno — v) + ab 

mit positivem b zu ersetzen sind). Die Größe der „Linienbreite“ b 
läßt sich an und für sich nach der Diracschen Strahlungstheorie 
berechnen. Sie ist jedoch sehr klein, und es kommt aus folgendem 
Grunde auf deren genaue Größe gar nicht an: nach einem Satze der 
klassischen Dispersionstheorie absorbiert eine Resonanzstelle aus 
einem Kontinuum eingestrahlter Frequenzen immer die gleiche Energie 
heraus, unabhängig von der Linienbreite ?2). In unserem Falle haben 
wir zwar nur eine einzige eingestrahlte Frequenz », dafür aber ein 
Kontinuum von Resonanzstellen, und das kommt auf das gleiche 
heraus. Wir werden daher nach Ausführung der Summation über n 
in Gl. (1,34) (die wir durch eine Integration ersetzen werden) zur 
Grenze b-» 0 übergehen. Einzig auf die mehr oder minder sprung- 
hafte Änderung der Absorption an den Absorptionskanten kann die 
Größe von b einen Einfluß haben. Die experimentelle Tatsache, 
daß sich die Absorption an den Kanten praktisch völlig sprunghaft 
ändert, zeigt, daß der Grenzfall b—> 0 eine gute Näherung darstellt 
(vgl. 8 5). 

Pr “= § 4. Einführung Hartreescher Näherung 


Um die Koeffizienten A,, nach Gl. (2,7) eingehend diskutieren 
zu können, führen wir Hartreesche Näherung ein, d.h. wir setzen 
die Eigenfunktionen u, der Kristallelektronen als Produkte von Eigen- 
funktionen y, des einzelnen Elektrons im Felde der Kerne und der 
übrigen Elektronen an. Die yw, nehmen wir als zueinander ortho- 
gonal und auf das Kristallvolumen V normiert an, gemäß ers 
(4,1) Yi Wat = 

J k k ; 
Ist N die Anzahl der Kristallelektronen, so sind die Zustände u, 
des gesamten Elektronensystems dadurch gegeben, daB jeweils N 
bestimmte Zustände w, mit je einem Elektron besetzt sind. Speziell 


sei die Eigenfunktion u, des Grundzustandes: A en 
(4,2) Uy = Wi (T1) We (Te) w;(t;) ... Wy (ty). 


Zum Unterschiede wollen wir im folgenden die im Grundzustande 
besetzten , mit lateinischem Index i und die im Grundzustande 


1) P.A.M. Dirac, The Principles of Quantum Mechanics, 2. Aufl. 
Oxford 1935, Kap. XI. 
= 2) Vgl. R.Becker, Theorie der Elektrizität II, $26, 8.11. = 


N 
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wobei hv,,= E,— E, die Differenz der den Zuständen w, und w, 
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unbesetzten yw, mit griechischem Index o bezeichnen. Die angeregten 
Zustände u, gehen aus (4,2) hervor, indem ein oder mehrere w, durch 
vorher unbesetzte w, ersetzt werden. Die Eigenfunktion desjenigen 
Zustandes u,, bei dem ein Zustand w, durch einen Zustand w, er- 
setzt ist, lautet also 
(4,3) = Wr (1) walte) Yw (tw), 
wobei an Stelle des einen Indexes n jetzt die beiden Indizes 4 und o 
treten. Für das Stromdichtematrixelement, das dem Übergang u,—>u,, 
zugeordnet ist, erhält man nach Gl. (1,17) bei Berücksichtigung von 
(4,1) die Gleichung 
(4,4) io = grad y* — y,* grad y,). 
Setzt man nämlich u,, für u, in Gl. (1,17) ein, so verschwinden dort 
in der Summe über k wegen der Orthogonalitat der wy, und y, alle Glieder 
außer dem mit k=i. Aus dem gleichen Grunde verschwinden die- 
jenigen Strommatrixelemente, die Übergängen aus demGrundzustande u, 
zu solchen angeregten Zuständen entsprechen, in denen mehrere Elek- 
tronen angeregt sind. Es kommen also nur Übergänge zu Zuständen 
in Betracht, bei denen jeweils nur ein Elektron angeregt ist}). (An 
diesen Ergebnissen ändert sich übrigens nichts, wenn man die Eigen- 
funktionen nicht einfach in Produktform gemäß (4,2) und (4,3), sondern 
als antisymmetrisierte Produkte in Determinantenform ansetzt, da 
unser Ansatz wegen der Orthogonalität der yw, und yw, und wegen 
der einfachen Besetzung schon an sich dem Pauliprinzip genügt.) 
Nach G].(2,7) ergibt sich nunmehr mit (4,4) für die Koeffizienten A, , 
der Ausdruck: 
+ dr dr’, 
h ”) 


zugehörigen Energiewerte ist. Wegen der 2fachen Spinorientierung 
ist in (4,5) ein Faktor 2 vor die Summe gezogen worden und die 
Summation über i und o ist nunmehr ohne Berücksichtigung der 
Spinorientierung zu verstehen. (Es ist ein Faktor 2 und nicht etwa 
ein Faktor 4 vorzuziehen, da Übergänge zwischen Zuständen mit 


1) Die Annahme orthogonaler y; ist nicht ganz streng gerechtfertigt; bei 
Berücksichtigung der „Fastorthogonalität“ würden, wie bereits H. Hönl be- 
merkt hat, auch Terme auftreten, die der Anregung von mehreren Elektronen 

leich ind, ürf üinfluß klei in. 
zugleich zugeordnet sind, doch dürfte deren Einfluß klein sein oe > 
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entgegengesetzter Spinorientierung nicht vorkommen, d. h. die zu- 
gehörigen Strommatrixelemente verschwinden) Die Summation 
über 2 in Gl. (4,5) ist über die im Grundzustande besetzten und die 
Summation über o über die im Grundzustande unbesetzten Zustände 
zu erstrecken; nur zwischen besetzten und unbesetzten Zuständen 
kommen also Übergänge vor, was hier zwanglos aus dem Ansatz 
für die Eigenfunktionen u, folgt. Die Gesamtheit der im Grund- 
zustande besetzten Zustände w, zerlegen wir in Gruppen entsprechend 
den verschiedenen (K-, L- usw.) Schalen der verschiedenen, das Gitter 
aufbauenden Atome, folgendermaßen: der Index s numeriere die Basis- 
atome (d. h. die Atome der 0-ten Gitterzelle); jedem Basisatom ent- 
spricht ein einfaches Teilgitter, das aus diesem durch Translation 


um die Gittervektoren hervorgeht. Es seien z. B. wi die K-Zustände 


des s-ten Teilgitters. Die Größe A,, nach Gl. (4,5) können wir dann 
folgendermaßen zerlegen: 


(4,6) Ang = (AX + Al 


8 


Die verschiedenen Anteile von A,, sind von der gleichen Form, wie 


das gesamte 4, nach Gl. (4,5) mit dem Unterschiede, daß die Sum- 


mation über ? jeweils nur über die betreffende Gruppe von Eigen- 
funktionen wy, zu erstrecken ist, z. B. bei dem Anteil ye den wir 
im folgenden als Beispiel näher betrachten werden, über die we. 


Wir haben bisher Eigenfunktionen y mit diskreten Indizes i, o numeriert, 
entsprechend der Tatsache, daß bei endlichem Kristallvolumen V die zu- 
gehörigen Energiewerte ein diskretes Spektrum darstellen. Da jedoch das 
Kristallvolumen groß gegenüber atomaren Dimensionen ist, liegen die Eigen- 
werte sehr dicht, und man pflegt daher die Eigenfunktionen y(f) und die zu- 
gehörigen Energieeigenwerte E(f) durch einen kontinuierlich veränderlichen 
Wellenvektor f zu charakterisieren. Die Eigenschaften der y (f) seien hier kurz 
zusammengestellt '): Die yw (f) stellen dreifachperiodisch modulierte ebene Wellen 
dar, sind also von der Form 


(4,7) = > aut 
Zur Normierung wählt man zweckmäßig die Bedingung We? 4 
l 


entsprechend einer mittleren Dichte von einem Elektron pro Volumeneinheit, 
so daß die w(f) vom Kristallvolumen unabhängig sind. Bei großen |f| haben 


1) Vgl. z.B. den Artikel von A. Sommerfeld u. H. Bethe im Handb. 
d. Phys. XXIV/2. Unser Wellenvektor f ist sr om als — we 
vektor bezeichnete. ag 
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die y (f) angenähert die Form ebener Wellen und für die zugehörigen Energie- 
werte gilt angenähert: 

h? 
An gewissen („Bragg-Brillouinschen‘“) Ebenen im f-Raum, die durch die 
Bedingung (f, b,) = sb charakterisiert sind, springt jedoch die Energie E(f), 


und zwar in dem Sinne, daß jeweils ein Energiebereich ausgelassen wird. Die 
dadurch entstehenden „verbotenen“ Energiebänder sind bei großen |f| sehr 
schmal; bei Übergang zu kleineren |f| werden sie breiter, und die „erlaubten“ 
Energiebänder werden entsprechend schmäler und ziehen sich schließlich bei 
den Zuständen, die den innersten Schalen entsprechen, praktisch zu scharfen 
Linien zusammen. Die diesen innersten Schalen zugehörigen Eigenfunktionen 
sind praktisch identisch mit den entsprechenden Eigenfunktionen der freien 
Atome (bzw. einfache Linearkombinationen derselben). 


Da uns hier gerade diese Elektronenzustände der innersten 
Schalen in erster Linie interessieren, werden wir für die im Grund- 
zustande besetzten wi usw. die diskrete Zählung beibehalten und 


nur die im Grundzustande unbesetzten yw, durch die kontinuier- 
lichen if) nach Gl. (4,7) ersetzen. Die Strommatrixelemente (4,4) 
nehmen dann die Form an 


(4,10) % f= (w(f) grad — grad w (f)) . 


Der Anteil = von A... entsteht, indem an Stelle von j,, in Gl. (4,5) 
der Ausdruck (4,10), und zwar genauer [wegen der verschiedenen Nor- 
mierung der w, nach (4,1) und der if) nach (4,8)] Fr je (f) ein- 


gesetzt wird. Zugleich tritt an die Stelle der Summation über o die 
Integration über den im Grundzustande unbesetzten Teil des f-Raumes, 
wobei noch mit der Dichte V (= Kristallvolumen) der Zustände im 
im f-Raum zu multiplizieren ist. — Unter we in (4,10) wollen wir 
die Eigenfunktionen des K-Zustandes des in der i-ten Gitterzelle 
befindlichen Atoms der Sorte s verstehen. Der Index 7 numeriert 
daher jetzt einfach die Gitterzellen. Man erkennt nun leicht, daß 


alle Summanden 27 in ne einander gleich sind; denn es gilt die 
Beziehung = we (t+ a,), wo a, ein Gittervektor ist. Ersetzt 
man im i-ten Summanden von AX : nach Gl. (4,5) zugleich r durch 
t+a, und vr durch r’+a,, so geht dort nach (4,10) und (4,7) ye in 
were und zugleich der Faktor 


woraus das Gesagte folgt. Wir dürfen daher alle Summanden ? durch 
den O-ten und die Summation über © durch Multiplikation mit der 


i (t, a) in sich iiber, 
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Anzahl V/4 der Gitterzellen ersetzen. Wir erhalten mit für den 
K,-Anteil von A, ER den Ausdruck: 


c 
ote 
An = dr, dr Ed - Ex. —h 


E®-Ex thy > 


wobei Ex, den der K-Zustände des s-ten 
bedeutet. Führen wir zur Abkürzung die Vektoren Me Oa 


(4,12) BM ES,,) = fis 3 


ein, so erhalten wir statt (4,11) }): 


4,13) | A E®-E,-hr 
(4, Pr 
Et — Ex, +hy 


in Streu- und Absorptionsanteil 
Ehe wir die Tensoreigenschaften der Koeffizienten A, diskutieren, 
müssen wir auf eine Schwierigkeit bei der Integration bezüglich £ in 
(4,13) eingehen, die darin besteht, daß der Integrand längs der Fläche 

im f-Raum 

(5,1) EB = Ez, +hv= 

singular ist. [Diese Schwierigkeit tritt natürlich nur dann auf, wenn 

die Fläche (5,1) in den Integrationsbereich fällt, d.h. wenn der 

Energiewert E, nicht vollständig in ein „verbotenes“ Energieband 

und auch nicht in den im Grundzustand u, von Elektronen besetzten 

Bereich des f-Raumes fällt] Wir denken uns die Fläche (5,1) von 

zwei benachbarten Flächen 

(5,2) und one 
E(f)= H,—«. 

(E, steht zur Abkürzung für Ex,+ hv) eingeschlossen, und unter- 

teilen den Integrationsbereich von f in zwei Teilbereiche, folgender- 


1) Die Schreibweise (4,13) ist strenggenommen nur dann richtig, wenn 
man voraussetzt, daß die Vektoren $, in (4,12) reell sind. Im II. Teil dieser 
Arbeit werden wir zwar auch komplexe $,„ in Betracht ziehen, doch werden 
deren Imaginärteile stets sehr klein sein, so daß sie in (4, 18) nicht nn. 
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(4,11) 
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maßen: der Integrationsbereich I sei der ganze Integrationsbereich 
(d.h. der im Grundzustand u, unbesetzte Teil des f-Raumes) aus- 
schließlich des Bereiches zwischen den beiden Flächen (5,2) und der 
Integrationsbereich II sei der von den Flächen (5,2) eingeschlossene 
kleine Bereich des f-Raumes. Wir zerlegen nun die Größe A,, in 
zwei Anteile: 

Str. Abs. 
(5,3) A, =AngtAne 


m 


. = 

Ke _ gKuStr. Ku, Abs. 


bzw. den als Beispiel betrachteten K,-Anteil von A 
631) 


= 


Dabei verstehen wir unter "rei diejenige Größe, die entsteht, wenn 
man in (4,13) die Integration bezüglich f über den Integrations- 
bereich I erstreckt und anschließend zur Grenze «—> 0 übergeht. 
Den Anteil a definieren wir (zunächst rein formal) als den 


Anteil von = , der entsteht, wenn in (4,13) die Integration über 
den Integrationsbereich JJ erstreckt wird, ebenfalls mit anschließen- 
dem Grenzübergang e—»> 0. Die Bezeichnung der beiden Anteile (5,3) 
als Streu- (Str.) und Absorptionsanteil (Abs.) wird im II. Teil dieser 
Arbeit gerechtfertigt werden. 

Die Integration über den Integrationsbereich I bereitet keine 
Schwierigkeiten. Bei der Integration über den Integrationsbereich II, 
also bei _ Abe tritt die erwähnte Schwierigkeit wegen des singulär 
werdenden Integranden auf. Diese Schwierigkeit verschwindet nun 
sofort, wenn wir uns an die Ausfiihrungen im § 3 erinnern, wonach 
auf Grund der Strahlungstheorie die Resonanznenner h(v, 9 — ») durch 
imaginäre Zusatzglieder 1b mit positivem b zu ergänzen sind. Ersetzt 
man demgemäß in (4,13) den Nenner 

(5,4,1) E(t) — Ex,— hy 


durch 


so tritt bei der Integration keine Schwierigkeit mehr auf, da die 
Singularitäten ins Imaginire gerückt sind. Die Integration über den 
Integrationsbereich II, d.h. über den Bereich zwischen den beiden 
Flächen (5,2) führen wir nun folgendermaßen aus: wir betrachten 
ein Volumenelement dieses Bereiches in Form eines kleinen Zylin- 
ders mit der Grundfläche do;, dessen Achse senkrecht auf der Fläche 
(5,1) steht (vgl. Abb. 1). Längs der Zylinderachse führen wir eine 
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Koordinate k, ein und zwar so, daß der Wert k = 0 auf der Fläche: 
(5,1) liegt und die positive k-Achse in Richtung wachsender Energie- 


werte E(f) weist; es ist dann 
aE 


wobei grad; den Gradientoperator im f-Raum bedeutet. Entwickelt 
man nun E(f) in eine Taylorreihe um k =0, so lauten bei Be- 
nutzung von (5,5) deren erste Glieder: 


E(t) = Ey + k\gradyE + 


Abb. 1. Das Volumenelement 
im Integrationsbereich II 


Abb. 2. Der Integrationsweg 
in der komplexen k-Ebene 


Bei der Integration über den Integrationsbereich II integrieren wir 
nun zunächst über die Koordinate k. Nach (4,13) und bei Ersetzung 
des Resonanznenners (5,4,1) durch (5,4,2) haben wir somit das Integral 
zu betrachten: 


Die Integrationsgrenzen. k = + 0 entsprechen dabei den Flächen (5,2), 
die ja abgesehen von Gliedern höherer Ordnung gleichen Abstand von 
(5,1) haben. Bei Benutzung der Taylorentwicklung (5,6) geht (5,7,1) 


_ k-|grad, EW|+ib+--- 


Nach $ 3 ist b eine kleine positive Zahl, die singuläre Stelle des 
Intreganden vom (5,7,2) liegt also auf der negativ-imaginären Achse 
nahe dem Ursprung der komplexen k-Ebene. Solange b + 0 ist, kann 
die Integration über k ohne Schwierigkeit längs der reellen k-Achse 


ais 
| 
Er 
dk 
5 
: 3 


a 


2 G. Moliére. Quantenmechanische Theorie usw. I 291 


geführt werden. Gehen wir nun entsprechend den Ausführungen des 
§ 3 zur Grenze b—> 0 über, so rückt die singuläre Stelle vom negativ- _ ER 
imaginären her auf die Stelle k = 0 zu. Der Integrationsweg muß 
daher bei diesem Grenzübergang so deformiert werden, daß er 2 
dieser Stelle in die positiv-imaginäre Halbebene hinein ausweicht 
(vgl. Abb. 2)). Geht man nun außerdem noch zur Grenze e-> 0 
und damit ö->0 über, (d. bh. läßt man die beiden Flächen (5,2) zu- 
sammenrücken), so bleibt von (5,7,2) das halbe negative Residuum 
an der Stelle k = 0 übrig, so daß wir erhalten: 


(5,7,3) 
"9 v 2) kjgrad, 


Fügen wir die fortgelassenen konstanten Faktoren und die über die 
Fläche (5,1) zu erstreckende Integration bezüglich do, wieder hinzu, 
so erhalten wir als vollständigen Ausdruck für den Absorptions- 


(6,8) grad, EM| 


Die Integration bezüglich do; ist dabei über die Fläche (5,1) zu 
erstrecken, sofern diese in einem im Grundzustande u, unbesetzten 
Teil des f-Raumes verläuft; andernfalls verschwindet Ax Abs. Der 
Vollständigkeit halber schreiben wir hier auch den Streuanteil von 
A™ an: 


mq 


| 


_ 2 fe (1,2, ); B"(E 8,)} 


Dabei ist die Integration bezüglich dz, über den Integrationsbereich I, 
) d.h. über den im Grundzustand unbesetzten Teil des f-Raumes zu 
erstrecken unter Auslassung der Fläche (5,1) in dem erläuterten Sinne. 
Die Vektoren 8” (f, &,,) in (5,8) und (5,9) sind durch Gl. (4,12) de- 
finiert. — In entsprechender Weise wie der K,-Anteil lassen sich 
auch die übrigen Anteile der Zerlegung ( (4,6) und somit das ge- 
| samte A, „ gemäß (5,3) in einen N und einen Absorptionsanteil 


zerlegen. 2 


1) Vgl. A. niceties u. G. ‘Schur, Ann. d. Phys. [5] 4. 8.409. 1930. 
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§6. Diskussion der Tensoreigenschaften der Koeffizienten 4,,, 


a) Allgemeine Beziehungen zwischen den Koeffizienten A,, 


Die Koeffizienten _ und = sind im allgemeinen Tensoren, deren 
Symmetrieeigenschaften wir im folgenden untersuchen wollen. Ihre K,-Anteile 
sind nach (5,9) bzw. (5,8) durch Volumen- bzw. Flächenintegrale im f-Raum 
definiert, wobei im Integranden Tensoren der Form {%; 8} auftreten. Ein 
solcher Tensor wird in den adjungierten (d. h. konjugiertkomplex-gespiegelten) 
Tensor überführt gemäß {2; B}+ = {B*; X*, durch Vertauschen des vor- 
deren und hinteren Vektors bei gleichzeitigem Übergang zum konjugiert- 
komplexen, wie man leicht durch Hinschreiben”der Tensorkomponenten sieht. 
Daraus folgt sofort, daß der zu = Str. nach GI. (5,9) adjungierte Tensor aus 
diesem durch Vertauschen von m und g hervorgeht. Das hier am Beispiel des 
K,-Anteils Gezeigte gilt auch für die übrigen Anteile und somit auch für die 
vollständige Größe A°'"-, d. h. es bestehen die Beziehungen 


mq? 


(6,1,1) = Str. 


qm? 

wobei das Zeichen + den Übergang zum adjungierten (d. h, konjugiertkomplex- 
gespiegelten) Tensor bedeutet. (Die „Spiegelung‘ bezieht sich dabei wohl- 
bemerkt auf die räumlichen Tensorkomponenten nach Achsen x, y, z und nicht 
etwa auf die Indizes m und q.) Insbesondere gilt für m = q, daß die Koef- 
fizienten > hermitische Tensoren sind, d. h. solche, die sich aus einem 
reell-symmetrischen und einem imaginär-antisymmetrischen Tensor zusammen- 
setzen lassen. Setzt man entsprechend (2,9) 


tr. Str. 
=Xn -q + 


und beachtet, daß gemäß (2,3) z, = ist, so sieht man, daß auch für die 


| die Beziehung besteht: 


Durch ganz entsprechende Uberlegungen findet man, daB wegen des von der 
Residuenbildung herrührenden imaginären Faktors in Gl. (5,8) für die Absorptions- 


bestehen, d.h. _ und An sind zueinander verkehrt adjungiert, und 


2 sind antihermitische Tensoren. 


speziell die 
Wir werden im folgenden eine Beziehung zwischen den Eigenfunktionen y (f) 


Er gebrauchen: da diese Lösungen einer reellen Wellengleichung sind, muß neben 


em sh y(f) auch das konjugiertkomplexe y (f)* eine Lösung sein, die zum gleichen 


Eigenwert E(f) gehört, und zwar ist offenbar 
(6,3) y(— f) = wif* 


me b) Kristalle mit Symmetriezentrum 


Wir betrachten nun den Sonderfall, daß der Kristall ein Symmetrie- 


ome zentrum besitzt. In diesem Fall gilt neben (6,3) fiir die Eigenfunktionen die 
ee weitere Beziehung: 


anteile der A,,, die Beziehungen 


q 
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Zusammen aus (6,4) und (6,3) folgt, daß für Kristalle mit Symmetriezentrum 
auch die Beziehung 


(6,5) y -D=y (f, r)* 
besteht, woraus fiir die Strommatrixelemente folgt: 
(6,6) id, = 


Wir betrachten nun den Vektor B** (£,8,,) nach Gl. (4,12). Nach Gl. (6,6) 
ergibt sich, daß in (4,12) jeweils die zu r und — r gehörigen Werte des Inte- 
granden zueinander konjugiertkomplex sind, so daß also ya (t, &,,) reell wird. 
Nach (5,9) und (5,8) erhalten wir somit das Ergebnis, das bei Kristallen mit 
Symmetriezentrum die Koeffizienten = (und wie man leicht zeigt auch die 


cer) simtlich reell und die Koeffizienten = rein imagindr sind. 


c) Sätze über die Diagonalelemente und Eigenwerte x 
der Tensoren A„m 


Nach Abschnitt a dieses Paragraphen sind die Koeffizienten _— bzw. 


oS. [Gl. (6,1,1) bzw. (6,1,2) für m = q] hermitische Tensoren und haben als 


mm 
solche reelle Diagonalelemente und Eigenwerte, und die Koeffizienten _ 


[Gl. (6,2) für m = q] sind antihermitische Tensoren und haben als solche rein 
imaginäre Diagonalelemente und Eigenwerte. (Die Worte „Diagonalelement“ 
und „Eigenwert‘“ beziehen sich auf die räumlichen Tensorkomponenten nach 
Koordinaten x, y, z und haben nichts mit den Indizes m und g zu tun.) Auf 
Grund der Gleichungen (5,8) und (5,9) lassen sich diese Aussagen noch bezüglich 
des Vorzeichens vervollständigen, und zwar können wir folgende Sätze aus- 
sprechen: 

a) Die Diagonalelemente und Eigenwerte der Koeffizienten _ sind bei 
genügend hohen Frequenzen » negativ (und streben mit »—»> © gegen Null); 
bei genügend niedrigen Frequenzen sind sie positiv. 

8) Die Diagonalelemente und Eigenwerte der Koeffizienten A 
stets negativ-imaginär. 


y) Die Imaginärteile der Eigenwerte von A,,„ bzw. C,,,, sind stets 


Abs. 
zen sind 


negativ-imaginär. (Das gleiche gilt auch für die „reduzierten“ Tensoren ;;, 
die im zweiten Teil dieser Arbeit eingeführt werden.) 

Zum Beweis dieser Sätze betrachten wir einen Tensor der Form {¥; ®*}, 
aus dem wir mit einem beliebigen (u. U. komplexen) Einheitsvektor e [d. h. 
einem solchen, für den (e,e*) = 1 ist], den Vektor |®;3*}e* und ferner das 
Skalarprodukt (die „Projektion“ des Tensors in Richtung e*) 


(6,7) (e, {B; B*} e*) = |(e, B) |? 


bilden. Dieser offenbar reelle und positive Ausdruck kann als Diagonalelement 
des Tensors {8 ; B*} bezüglich der durch den Einheitsvektor e definierten Achsen- 
richtung aufgefaßt werden. Daraus folgen mit den Gl. (5,8) und (5,9) für m=y 
bereits die Sätze a und f, soweit sie sich auf Diagonalelemente beziehen. Zum : 
Beweis der Aussagen über Tensoreigenwerte fügen wir noch folgende Betrachtung E \ 


an: Es sei A ein beliebiger (u. U. nichthermitischer) Tensor, A einer seiner were 

Eigenwerte und q der zugehörige Eigenvektor. Es gelte also die Gleichung: ur 
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durch skalare Multiplikation mit q* folgt weiter: ER 


(q*, Aq) = A\q\?, 


woraus folgt, daß jeder Eigenwert eines Tensors als ein Diagonalelement 
bezüglich der durch den zugehörigen Eigenvektor definierten Achsenrichtung 
aufgefaßt werden kann. Dies gilt insbesondere auch dann, wenn die zu den 
verschiedenen Eigenwerten gehörigen Eigenvektoren nicht aufeinander senkrecht 
stehen. Daraus folgen die Aussagen der Sätze « und @ über Eigenwerte, sowie 
mit Benutzung der Zerlegung A,, „= + der Satz y. 


d) Ein Spezialfall 


\ Br Wir wollen nochmals GI. (6,3) benutzen und zwar jetzt für den all- 
= 


meinen Fall eines Kristalls mit oder ohne Symmetriezentrum. Wir betrachten 
wieder den Vektor w(t, &,,) nach Gl. (4,12), und zwar wollen wir den Spezial- 
fall behandeln, daß der Faktor erriAm,t im Integranden von (4,12) in dem 
wesentlich in Betracht kommenden Integrationsbereich (d. h. in dem Bereich, 
in dem die K-Eigenfunktion der s-ten Atomsorte wesentlich von Null verschieden 
ist) — = 1 gesetzt werden kann. Dies ist dann der Fall, wenn 


dn = 78.1 © | genügend groß ist gegenüber dem Bohrschen Radius der s-ten 


Atomsorte. Die Vektoren Be (t, 8.) sind in diesem Falle näherungsweise 
unabhängig von &,,, also 


&,,) = (t, 0). 
Es gilt daher für sie nach (6,3) die Beziehung ug 
(6,8) = BX" (-1,8,), 


woraus weiter folgt, daß in (5,9) und (5,8) jeweils die zu f und —f gehörigen 
Werte des Integranden zueinander konjugiertkomplex sind. Es ergibt sich 
also, daß die AS sofern sie der gemachten Voraussetzung genügen, sämtlich 
einander gleich und reell sind, und daß die AÄ"* unter der gleichen Bedingung 
einander gleich und rein imaginär sind. 

Nehmen wir, indem wir unseren Spezialfall noch weitergehend spezialisieren, 
an, daß die rg verschwinden oder vernachlässigbar klein sind; dann können 
alle A,,, durch den reellen symmetrischen Tensor ny ersetzt werden und 


die Größen C,,, nach Gl. (2,9) nehmen die Form an 


(6,9,1) 400 


oder abgekürzt, da jetzt C,, „ nur von der Differenz m — q abhängt: 


(6,9,2) One = 

In das Gleichungssystem (2,13) eingesetzt ergibt dies die Grundgleichungen 
in einer Form, in der sie auf Grund anderer Voraussetzungen von L. Posener') 
aufgestellt und diskutiert wurden. 


pe 
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— 


DL. Posener, Ann. d. Phys. [5] 19. 8. 849. 1934. 
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§7. Frequenzabhängigkeit der Koeffizienten 4A,,, 


Es bleibt noch die Abhängigkeit der Koeffizienten A) und 4," 
ent von der Röntgenfrequenz » zu diskutieren. Vorweg bemerken wir, 


daß bei großen die Vektoren B“(f, &,) der Gl.(4,12) verschwinden, 
= weil dort die Eigenfunktionen w(f) und daher nach Gl. (4,10) auch 
wie die Strommatrixelemente m) im Integranden von (4,12) räumlich 
rasch oszillieren. 
Wir betrachten nun zunächst den Anteil AN Abs. der nach 
Gl. (5,8) durch ein Integral über die Fläche (5,1): 
= Ef) = Ex +h» 
fe im f-Raum dargestellt wird, sofern diese in einen im Grundzustand 
ich unbesetzten Teil des f-Raumes fällt. Bei Frequenzen » unterhalb 
a der zu den Atomen der Sorte s gehörigen K-Kante »;, verschwindet 
enn der Absorptionsanteil A” ‚bs weil dort die Fläche (5,1) in Bereiche 
ten fällt, die im Grundzustand von Elektronen besetzt sind. Dicht 
pise oberhalb von vx, nimmt re Abe. seine größten Werte an und fällt 


bei weiter wachsender Frequenz » ab, wobei wegen der an den 
Bragg-Brillouinschen Ebenen im f-Raum auftretenden Energie- 
sprünge unregelmäßige Schwankungen auftreten‘). Bei großen Fre- 


quenzen » schließlich verschwindet A,” Abs weil dort die Fläche (5,1) 


in Bereiche großer |fj fällt, wo nach der obigen Bemerkung die 
Vektoren ®#(f,®,) und daher der Integrand von (5,8) verschwinden. 


Wir betrachten nun den Streuanteil A* gay nach Gl. (5,9). 


Danach sieht man, daß bei kleinen Frequenzen »—> 0 _ Str. einem 


endlichen positiven Grenzwert zustrebt. (Die Bezeichnung „positiv“ 
en, bezieht sich dabei auf die Diagonalelemente bzw. Eigenwerte des 


Tensors Au 86, Abschnitte). Bei Annäherung der Frequenz» 
in 


an die K-Kante vx, steigt ty “*- an und erreicht etwa bei » = vx, 
seinen größten Wert. Bei Überschreiten von vx, wechseln Teile 
| des Integranden von (5,9) das Vorzeichen, und regi fallt daher 

| ab und wird schließlich negativ. Bei weiter wachsender Frequenz » 
geht schließlich _ Str. nach Überschreiten eines Minimalwertes, ohne 


en nochmals das Vorzeichen zu wechseln, gegen Null, weil bei großen » 
r') der Teil des Integranden, bei dem der Resonanznenner klein ist, 


ts i 1) Vgl. R. de L. Kronig, a. a. O. ; 
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in Bereiche großer |f| fällt, wo die Vektoren ®“ verschwinden, 
während der Bereich des Integranden bei kleinen |f|, wo 8% merklich 
von Null verschieden ist, wegen des dort großen Resonanznenners 
nichts beiträgt. 

Bei großen Frequenzen » verschwindet also sowohl der Streu- 
wie der Absorptionsanteil von A,,,; dort sind also die Tensoren C,,, 
ir u G1.(2,9) durch die skalaren Größen y,,_, zu ersetzen, wodurch 
die Grundgleichungen (2,13) ihre urspriinglich von M. v. Lene auf- 
gestellte Form annehmen 


1) Vgl. M. v. Laue, a. a. O. 


Berlin-Dahlem, Kaiser-Wilhelm- Institut für 
Chemie und Elektrochemie. 
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